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Abstract 


Landform  morphometry  is  related  to  causative 
factors  of  climate,  mantle  characteristics,  vegeta¬ 
tion  density,  and  lithology.  Drainage  basins  ana¬ 
lyzed  are  in  mature  fluvial  development,  free  from 
obvious  structural  influence,  and  in  climates  rang¬ 
ing  from  arid  to  humid.  Over  80  basins  in  Arizona, 
Colorado,  New  Mexico,  and  Utah  were  inspected  in 
the  field;  22  were  subjected  to  detailed  field  inves¬ 
tigations. 

Morphometric  properties  of  all  basins  were 
measured  from  recent  1:  24,000  topographic  maps, 
or  special  drainage  maps,  and  show  great  ranges: 
valley-side  slope  has  a  four-fold  range,  drainage 
density  and  basin  relief  about  100-fold,  and  basin 
area  5,000-fold. 

Climate  was  measured  by  the  Thornthwaite 
precipitation- effectiveness  (P-E)  index,  and  the  in¬ 
tensity  of  average  five-year,  one-hour  rains.  Lith¬ 
ologies  distinguished  are  shale,  clastic  rocks, 
schist,  limestone,  granitic  rocks,  and  acid  volcanic 
rocks,  assumed  to  indicate  roughly  the  proportion 
of  clay  in  residual  soils.  Mantle  and  cover  charac¬ 
teristics  are  summarized  by  relative  infiltration 
capacity,  wet  and  dry  soil  strength,  percentage  of 
bare  area,  and  surface  roughness. 

Validity  of  statistical  inferences  is  affected  by 
sampling  procedures,  which  were  mainly  non- 
random.  Analysis-of-variance  tests  show  that 
shale  and  schist  produce  high  drainage  densities, 
whereas  granitic,  clastic,  and  acid  volcanic  rocks 
produce  low  densities,  with  no  significant  differ¬ 
ences  among  these  three.  Shale  and  schist  have 
low  infiltration  capacity,  whereas  clastic  and  gran¬ 
itic  rocks  produce  high  values.  Valley-side  slopes 
are  greatest  in  shale  and  limestone,  and  lowest  in 
schist,  on  the  average,  subject  to  differences  in 
erosional  history. 

Drainage  density  averaged  highest  in  basins 
having  gullied  slopes,  intermediate  in  unmodified 
basins,  and  lowest  in  basins  with  trenched  main 
channels.  Differences  were  not  judged  significant, 
even  with  adjustment  for  infiltration,  though  low 
power  renders  the  test  inconclusive. 


Using  punched  cards,  variables  were  analyzed 
in  pairs  for  significant  correlations  and  regres¬ 
sions.  Texture  of  topography  decreases  exponenti¬ 
ally  with  valley-side  slope.  Relative  relief  and 
ruggedness  number  corrected  for  order  increase 
exponentially  with  slopes.  Texture  questionably 
varies  inversely  with  relative  relief. 

High  correlations  exist  among  channel  fre¬ 
quency,  drainage  density,  and  frequency  and  den¬ 
sity  of  first-order  channels.  It  is  inferred  that 
average  channel-segment  length  of  any  order  var¬ 
ies  inversely  with  the  total  channel  length,  for 
areas  of  the  same  size  but  differing  drainage  den¬ 
sities.  Average  first-order  channel  length  is  a  con¬ 
stant  proportion  of  average  channel  length  of  all 
orders,  and  is  inversely  proportional  to  total  chan¬ 
nel  length. 

Multiple  correlation  analyses  show  that  valley- 
side  slopes  and  drainage  density  are  related  to 
climate,  and  to  properties  of  mantle  and  vegetal 
cover,  as  measured  in  the  field.  Slopes  increase 
with  infiltration  capacity  and  P-E  index,  but  vary 
inversely  with  wet  soil  strength  and  runoff  inten¬ 
sity-frequency.  Greater  runoff  produces  lower 
valley-side  slopes.  Climate  enters  further  by  the 
increase  of  infiltration  capacity  with  P-E  index. 
Greater  infiltration  permits  interflow,  reducing 
the  ratio  of  slope  runoff  to  channel  flow.  Channels 
are  deepened  with  less  opportunity  for  erosion  of 
slopes,  producing  steep  slopes. 

Drainage  density  varies  with  percent  bare 
area  and  runoff  frequency-intensity,  but  inversely 
with  P-E  index  and  infiltration  capacity,  fitting 
theoretical  observations  of  Gilbert,  Davis,  and 
Johnson,  and  confirming  Horton’s  infiltration  the¬ 
ory  of  erosion.  Primary  controls  of  topographic 
texture  are  lithology  and  climate,  which  act 
through  the  agency  of  secondary  surface  proper¬ 
ties.  Accidental  changes  in  these  secondary  con¬ 
trols,  unrelated  to  climatic  changes,  will  quickly 
alter  the  drainage  density  to  another  value  than 
that  expected  from  consideration  of  climate  and 
lithology  alone. 
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Introduction  and  Acknowledgements 


Climatic  influence  on  topographic  form  proper¬ 
ties  of  fluvially  dissected  landscapes  has  long  been 
recognized  by  geologists  and  geographers.  Al¬ 
though  considerable  literature  exists  on  this  sub¬ 
ject,  three  vital  procedures  are  lacking  in  early 
studies:  (1)  Isolation  of  climatic  elements  and  ef¬ 
fects.  Climate  was  usually  referred  to  as  “arid” 
or  “humid,”  in  regard  to  moisture;  “hot,”  “tem¬ 
perate,”  or  “cold,”  in  regard  to  temperature.  (2) 
Investigation  of  relative  importance  of  processes. 
Knowledge  of  process  was  largely  taken  for  granted 
emphasis  being  placed  on  historical  interpretation 
and  classification.  (3)  Expression  of  laws  in  quan¬ 
titative  terms.  Research  was  largely  of  a  qualita¬ 
tive  nature  because  it  was  thought  that  features  of 
the  landscape  are  so  variable  that  mathematical 
treatment  would  be  impracticable.  Attempts  at 
quantitative  application  of  principles  were  resticted 
to  simple  cases. 

Recently,  however,  a  number  of  geomorpholo¬ 
gists  have  recognized  the  availability  of  mathema¬ 
tical  tools  capable  of  coping  with  the  variability, 
or  chance  element,  in  landforms.  Their  approach 
has  been  chiefly  to  examine  in  detail  the  variations 
in  properties  of  areas  similar  in  major  character¬ 
istics  and  history.  If  it  can  be  shown  that  these 
variations  are  too  great  to  be  of  random  nature, 
then  they  must  be  the  result  of  some  basic  differ¬ 
ences  in  rate  or  intensity  of  operation  of  process. 

In  these  recent  studies,  climate  has  been  for  the 
most  part  of  secondary  interest. 

The  main  purpose  of  this  study  is  to  provide  a 
foundation  for  further  research  into  the  effects  of 
availability  of  moisture  to  vegetal  cover,  and  of 
rainstorm  intensity  and  frequency  on  processes  of 
landform  evolution.  Consideration  is  given  to  pro¬ 
perties  of  the  surface  mantle  as  they  are  affected 
by  climate,  and  in  turn  as  they  affect  processes. 
Because  action  of  the  processes  is  too  slow  to  be 
observed  during  a  short  time,  relative  rates  are 
inferred  from  measurements  of  the  landforms 
themselves. 

The  factors  influencing  topographic  texture  and 
average  slope  can  be  considered  as  fundamental  in 
control  of  the  fluvial  landscape.  The  controls  most 
often  mentioned  in  descriptive  treatments  of  the 
subject  are  “climate,”  “rock  permeability,”  and 
“rock  strength.”  Structural  properties  of  the  rock 
such  as  jointing,  initial  or  available  relief,  stage 
in  the  erosional  cycle,  and  condition  of  the  surface 
soil  are  frequently  mentioned  as  well  (Thornbury, 


1954,  p.  128).  Factors  influencing  the  steepness 
and  mode  of  retreat  of  slopes  in  fluvially- controlled 
landscapes  have  received  much  attention  in  testing 
Walther  Penck’s  ideas  of  slope  development. 

Gilbert  (1877,  p.  109-118)  analyzed  the  processes 
of  slope  and  channel  equilibrium  as  effected  by  run¬ 
ning  water.  He  reasoned  that  fluvially  eroded  areas 
in  arid  regions  are  “degraded”  at  a  rate  limited  by 
the  slowness  of  the  production  of  fine  waste  mate¬ 
rial;  removal  by  transport  is  rapid  once  waste  has 
formed.  Thus,  differences  in  rock  texture  are  accen¬ 
tuated  in  arid  regions.  Greater  availability  of  water 
hastens  the  rate  of  production  of  waste  material, 
and  likewise  increases  the  rapidity  of  transport. 

But  the  greater  vegetal  cover  in  humid  regions 
slows  transportation  and  is  generally  effective  in 
causing  an  accumulation  of  waste.  Thus,  differences 
in  rock  texture  are  obscured  by  accumulated  waste 
in  humid  regions.  A  climate  intermediate  between 
arid  and  humid,  Gilbert  concluded,  must  produce 
the  most  rapid  and  effective  erosion. 

In  badlands  areas,  declivities  decrease  as  the 
amount  of  water  flowing  over  them  increases,  and 
slopes  are  greater  “in  proportion  as  they  are  near 
divides,  unless  they  are  very  near  divides”  (Gil¬ 
bert,  1877,  p.  117).  The  convexity  of  hilltops  evi¬ 
dently  puzzled  Gilbert  at  the  time  he  wrote  the  re¬ 
port  on  the  Henry  Mountains,  and  his  solution  to  the 
problem,  as  indicated  later,  was  the  agency  of  creep 
(Gilbert,  1909,  p.  344-350).  Little  has  been  added  to 
the  qualitative  knowledge  of  slope  equilibrium  since 
Gilbert’s  publications;  quantitative  study  of  these 
processes  described  has  only  just  begun. 

Walther  Penck  (1953,  p.  119-120,  and  130-133) 
evidently  believed  that  the  same  landforms  can  ap¬ 
pear  under  all  climates,  though  developing  at  dif¬ 
ferent  rates.  Climatic  differences  may  change  the 
relationship  of  mass-wasting  to  transportation,  but 
they  do  not  produce  any  fundamental  change  in  the 
denudational  forms.  Slope  gradients  depend  only  on 
the  intensity  of  erosion;  so  if  all  slopes  in  a  region 
are  of  the  same  steepness,  it  is  an  expression  of  the 
fact  that  erosion  has  acted  upon  them  with  the  same 
intensity.  It  follows  from  this,  according  to  Penck, 
that  gradients  are  determined  by  endogenetic  fac¬ 
tors,  and  less  so  by  rock  hardness.  The  form  of 
slope  profiles  reflects  changes  in  endogenetic  rates. 
Old  highlands,  separated  from  a  former  base-level 
of  erosion  by  a  topographic  break,  invariably  dis¬ 
play  concave  slopes.  Convex  and  straight  slopes  are 
associated  with  zones  of  intense  recent  erosion. 
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Where  the  mean  slope  gradient  corresponds  to  the 
maximum  gradient  for  the  rock  concerned,  convex 
slopes  are  replaced  by  straight  profiles;  these  find 
their  “purest  expression  in  the  steep  relief  of 
mountains  which  escaped  Pleistocene  glaciation, 
whether  these  are  covered  with  tropical  virgin  for¬ 
est  or  belong  to  arid  regions”  (Penck,  1953,  p.  132). 
Penck’s  views  of  relation  of  slope  form  to  endogene- 
tic  factors  have  largely  been  discredited  in  univer¬ 
sal  application.  It  is  possible  that  qualifications 
may  have  originally  been  present  in  his  book  that 
have  been  lost  in  translation,  so  great  is  the  mass 
of  verbiage. 

The  first  truly  quantitative  approach  to  the  sub¬ 
ject  of  the  control  of  topographic  texture  was  by 
hydrologists.  Robert  E.  Horton  (1945,  p.  283,  285) 
defined  the  terms  drainage  density,  the  total  length 
of  channels  divided  by  the  area  drained,  and  stream 
frequency,  the  number  of  channel  segments  divided 
by  the  area  drained,  showing  that  together  they  de¬ 
termine  the  fineness  of  topography.  He  concluded, 
perhaps  erroneously,  that  areas  with  low  drainage 
density  are  poorly  drained,  and  areas  with  high 
drainage  density  are  well  drained.  Horton  em¬ 
phasized  the  importance  of  surface  characteristics 
in  determining  drainage  density,  particularly  in¬ 
filtration  capacity,  and  resistance  of  the  surface 
to  erosion  (1933,  p.  260,  and  1945,  p.  284). 

Forest  and  range  conservationists,  working  un¬ 
der  the  U.S.  Department  of  Agriculture,  have  ob¬ 
tained  much-needed  quantitative  information  on 
infiltration  and  erosion  rates  related  to  land  use, 
type  of  soil  and  cover,  and  distribution  of  rain  dur¬ 
ing  storms  (Dortignac,  1951;  Love  and  Johnson, 

1952;  Love  and  Dunford,  1952;  Love,  1953;  Anderson, 
1954).  Of  special  importance  to  geomorphologists 
are  studies  of  rainfall- runoff  proportions,  in  cases 
of  known  accelerated  erosion,  in  both  forested  and 
grazed  watersheds.  Schumm  (1955)  has  demon¬ 
strated  that  sediment  loss  from  basins  can  be  re¬ 
lated  to  form  properties  within  the  basin,  and  to 
infiltration  capacity  of  the  soil. 

A.  N.  Strahler  has  emphasized  the  use  of  statis¬ 
tical  methods  in  quantitative  investigations,  rather 
than  deductive  analysis  (1950,  1954a).  Recently  he 
has  used  the  Pi  Theorem  and  dimensional  analysis 
in  combining  factors  that  rationally  can  be  said  to 
contribute  to  the  main  causes  determining  drainage 
density  (Strahler,  1956b).  Dimensionless  numbers 
were  derived  that  have  the  property  of  independence 
of  size  of  corresponding  units  of  the  drainage  basins 
considered.  Thus,  if  drainage  density  can  rationally 
be  said  to  depend  upon  (1)  the  average  rate  of  flow 


stated  in  volume  per  unit  of  area  of  basin,  (2)  an 
erosion  proportionality  factor,  (3)  the  height  of  the 
basin  above  its  mouth,  (4)  density  of  surface  ma¬ 
terial,  (5)  viscosity  of  water,  and  (6)  the  accelera¬ 
tion  of  gravity,  then  these  quantities  can  be  arranged 
into  four  dimensionless  numbers.  These  are  (1)  rug¬ 
gedness  number  (product  of  drainage  density  and 
relief  of  basin),  (2)  Froude  number,  (3)  Reynolds 
number,  and  (4)  Horton  number  (product  of  erosion 
proportionality  factor  and  the  average  volume  rate 
of  flow  per  unit  area).  Full  evaluation  of  the  rela¬ 
tions  among  all  these  quantities  would  take  many 
years  of  research.  The  present  study  treats  only 
portions  of  these  relations  that  are  inherent  in  ex¬ 
pression  of  relief,  strength  of  surface  material, 
average  rate  of  runoff,  measures  of  roughness,  and 
extent  of  cover. 

Generalizing  from  these  discussions,  there 
seems  to  be  agreement  on  factors  influencing  slope 
development  and  equilibrium,  and  drainage  density. 
Considerable  disagreement  exists  as  to  which  fac¬ 
tors  are  dominant  in  effect  and  which  are  secondary. 
Where  topographic  features  and  related  quantities 
are  measurable,  it  is  generally  not  known  whether 
these  quantities  are  related  by  simple  proportion, 
by  power  functions,  by  exponential  functions,  or  by 
other  transcendental  functions.  The  writer  has  en¬ 
deavored  to  advance  the  understanding  of  these  fun¬ 
damental  questions  by  using  (1)  form  measurements 
from  topographic  maps,  (2)  simple  field  measure¬ 
ments  of  soil  properties  and  characteristics  of 
vegetal  cover,  (3)  analysis  of  climate  elements,  and 
(4)  combination  of  numerical  data  by  statistical 
analysis.  Correlation  analysis,  in  particular,  allows 
separation  of  major,  minor,  and  definitely  unre¬ 
lated  effects. 

This  investigation  formed  part  of  a  quantitative 
study  of  erosional  landforms  sponsored  through  Co¬ 
lumbia  University  by  the  Geography  Branch  of  the 
Office  of  Naval  Research  as  Project  Number  NR 
389-042,  under  Contract  N60NR  271,  Task  Order  30. 

The  writer  is  indebted  to  Professor  A.  N. 

Strahler  for  guidance  in  planning  and  carrying  out 
field  and  map  research.  Professor  Howard  Levene 
of  the  Department  of  Mathematical  Statistics  of 
Columbia  University  advised  on  a  method  of  analy¬ 
sis  of  variance  and  critically  read  the  text.  Mr. 

Gary  Gordon  of  Mt.  Vernon,  Iowa,  assisted  in  the 
field  during  two  months  of  the  summer  of  1955. 

Sara  Melton  assisted  during  two  field  seasons,  and 
helped  in  preparing  the  manuscript.  Styron  Douthit 
and  Steven  Haines  of  the  University  of  Oklahoma 
aided  in  preparation  of  the  illustrations. 


Methods  of  Map  Analysis 

Geologists  familiar  with  quantitative  landform 
analysis  will  understand  most  of  the  terms  used 
here,  although  some  are  original.  To  avoid  misun¬ 
derstanding,  all  terms  and  symbols  are  defined  in 
Table  1,  which  gives  units  and  dimensions  of  all 
quantities.  Map  data  is  included  from  3  fifth-order 
basins,  22  fourth-order  basins,  and  83  third-order 
basins.  Field  measurements  are  given  for  1  fifth- 
order  basin,  19  fourth- order  basins,  and  9  third- 
order  basins.  These  are  divided  into  58  third-order 
basins  for  analysis.  The  remaining  basins  were  in¬ 
spected  briefly  to  determine  the  general  condition 
of  the  vegetal  cover,  rock  type,  and  condition  of 
channels. 


MAP  ACCURACY 

The  amount  of  usable,  reliable  information  con¬ 
tained  on  topographic  maps  increases  with  the  ac¬ 
curacy  of  the  map.  Because  maps  differ  in  accuracy 
depending  on  scale  and  contour  interval,  and  on  the 
method  by  which  they  were  made,  it  is  important  in 
comparative  studies  that  data  be  obtained  from 
maps  of  comparable  relative  accuracy.  Otherwise, 
comparisons  between  areas  differing  in  some  re¬ 
spect  will  be  in  error. 

Only  tho  se  maps  published  by  the  U.  S.  Geologi¬ 
cal  Survey  at  the  scale  1:24,000,  and  compiled 
from  aerial  photographs,  have  been  used  in  this 
study,  with  three  exceptions:  (1)  maps  of  the  Neder¬ 
land,  Black  Hawk,  and  Morrison  quadrangles  in 
Colorado,  were  published  at  the  scale  1:24,000,  but 
were  made  by  plane-table  surveys  in  1942  and  1938, 
before  photogrammetric  methods  were  in  use  by 
the  U.  S.  Geological  Survey.  The  drainage  nets  of 
the  basins  studied  in  these  areas  were  drawn  from 
aerial  photographs  obtained  from  the  U.S.  Forest 
Service,  by  using  a  KEK  Stereoscopic  Plotter.* 
Considerable  revision  in  the  drainage  nets  of  basins 
in  those  areas  was  necessary.  (2)  A  small  badlands 
area  in  the  Chinle  shale,  north  of  Cameron,  Ari¬ 
zona,  is  not  included  in  any  quadrangle  map,  and 
could  not  be  adequately  shown  at  the  scale  1:24,000. 
Two  small  drainage  basins  in  the  badlands  were 
mapped  at  1  inch  to  50  feet  by  E.  D.  Koons,  S.  A. 
Schumm,  and  A.  N.  Strahler,  in  the  summer  of 
1953.  When  examined  by  the  writer,  two  years  later, 
two  small  channels  had  become  obliterated  and 

*Manufactured  by  Kail  and  Associates,  Denver,  Colorado. 
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were  removed  from  the  drainage  net;  otherwise, 
the  map  was  extremely  accurate.  (3)  An  area  of 
badlands  in  the  Saguaro  National  Monument,  south¬ 
east  of  Tucson,  Arizona,  is  covered  only  by  the 
1904  Tucson  Quadrangle,  at  the  scale  1:125,000.  In 
order  to  use  the  area  in  the  study,  a  fifth-order 
basin  was  mapped  at  1  inch  to  100  feet  with  plane 
table  and  alidade  by  the  writer  and  assistant,  in 
December,  1955.  The  drainage  nets  mapped  in 
these  three  areas  are  believed  comparable  in  ac¬ 
curacy  to  the  drainage  nets  shown  on  photogram¬ 
metric  maps  of  other  areas. 

In  all  basins  in  which  field  measurements  were 
taken,  the  drainage  net  was  drawn  on  the  map  to 
correspond  with  the  existing  channels.*  It  was 
found  that  in  most  maps  a  high  percentage  of  the 
first-order  channels  is  represented  by  the  small¬ 
est  cusps  in  contour  lines.  In  a  few  areas  (e.g. 

Spring  Creek  basin,  Springville  Quadrangle,  Utah), 
the  smallest  cusps  represent  only  irregularities  in 
the  surface,  not  persistent  drainage  lines.  In  other 
areas  (e.g.  Quartz  Hill  basin,  Paxton  Springs 
Quadrangle,  N.  M.),  the  smallest  cusps  in  the  con¬ 
tour  lines  represent  second  or  third-order  channels. 
In  many  basins  shallow  gullies  established  on  the 
valley  sides  are  not  shown  by  the  contour  lines.  In 
accordance  with  these  considerations,  about  95% 
of  the  basins  studied  needed  at  least  minor  correc¬ 
tions  of  the  drainage  net  as  shown  on  the  contour 
map.  (See  Appendix  IV  for  photo- copies  of  portions 
of  maps,  showing  drainage  nets  of  all  basins  studied.) 

METHODS  OF  OBTAINING  QUANTITATIVE  DATA 

FROM  MAPS 

Many  of  the  basins  studied  are  rather  small,  as 
shown  on  the  original  maps.  In  measuring  features 

♦Horton  (1945)  used  stream  with  the  same  meaning  that 
channel  has  here  and  in  previous  reports  by  A.  N. 

Strahler  (1950,  1952a,  1953),  V.  C.  Miller  (1953),  K.  G. 
Smith  (1953),  S.  A.  Schumm  (1954)  and  D.  R.  Coates 
(1955).  The  term  channel  means  a  permanent,  clearly 
defined  trench  or  trough  clearly  showing  evidence  of 
scour  by  channel  flow  and  bounded  by  valley  sides  slop¬ 
ing  toward  the  channel  axis .  Excluded  by  this  definition 
are  ephemeral  or  seasonal  rills  cut  in  otherwise  smooth 
slopes  which  do  not  contribute  runoff  to  the  rill  and 
which  have  no  clearly  defined  watershed  and  drainage 
divide.  Channel  is  preferred  since  emphasis  is  placed 
on  the  depression  in  which  water  flows,  rather  than  a 
stream  of  water  itself.  Most  channels  studied  here  would 
be  classed  as  ephemeral  or  intermittent  streams. 
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in  these  basins,  the  percentage  of  error  can  be  re¬ 
duced  by  enlarging  the  maps  photographically; 
measurements  on  the  prints  can  then  be  made  with 
greater  relative  precision.  The  maps  of  all  basins 
studied  in  detail  in  the  field  were  enlarged  two  to 
eight  times,  depending  upon  original  size.  Prints 
were  made  by  a  Photostat  machine,  or  by  making 
full-sized  negatives  on  large  sheets  of  film  and 
printing  directly  on  blue-line  paper.  The  prints 
were  air  dried.  Distortion  of  the  image  was  present 
in  almost  all  cases,  but  by  marking  several  one- 
inch  intervals  with  various  orientations  on  the  or¬ 
iginal  and  measuring  the  corresponding  interval  on 
the  enlarged  print,  it  was  established  that  differ¬ 
ential  shrinkage  caused  slight  variation  in  scale. 

In  two  perpendicular  directions,  parallel  to  the 
edges  of  the  print,  the  greatest  difference  in  scale 
was  1.6%. 

Measurements  of  length  were  made  from  the 
enlarged  maps  by  means  of  a  Dietzgen  Map  Meas¬ 
ure,  No.  1719B,  scaled  to  1/32  inch,  and  on  which 
measurements  can  be  estimated  to  the  nearest  1/64 
inch.  Hence,  on  an  enlarged  map  with  scale  of  about 
10  inches  per  mile,  measured  lengths  can  justifi¬ 
ably  be  recorded  to  the  nearest  0.01  mile. 

Area  measurements  were  made  with  a  Keuffel 
and  Esser  Compensating  Polar  Planimeter,  Model 
4242.  Highly  consistent  results  were  obtained  in 
repeated  measurements  of  the  same  outline.  Re¬ 
versing  the  position  of  the  pole  with  respect  to  the 
tracer  head  eliminated  a  small  systematic  error. 
Deviations  greater  than  0.02  square  inches  are  un¬ 
usual  in  individual  area  measurements,  hence,  on 
a  map  with  scale  100  square  inches  per  square 
mile,  it  is  justifiable  to  record  measurements  to 
the  nearest  0.001  square  mile. 


PROPERTIES  MEASURED  FROM  MAPS 

Of  all  the  information  that  can  be  obtained  from 
maps  by  measurement,  two  general  categories  may 
be  distinguished:  (a)  quantities  measured  (or 
counted)  directly,  with  no  more  computation  than 
multiplication  by  a  scale  factor,  and  (b)  quantities 
that  must  be  computed  from  direct  measurements 
because  they  cannot  in  themselves  be  measured. 

The  first  category  includes  length  and  area  meas¬ 
urements,  numbers  of  items  obtained  by  enumera¬ 
tion,  and  certain  angular  measurements.  Lower¬ 
case  Roman  letters  are  used  as  notation  for 
properties  of  length,  area,  or  number;  lower-case 
Greek  letters  for  angular  measurements.  To  de¬ 
velop  a  satisfactory  notation  for  the  numerous 
morphometric  properties,  full  use  of  upper  and 
lower-case  letters,  and  of  subscripts,  is  necessary. 
Subscript  notation  must  be  systematic  and  consist¬ 
ent.  Subscripts  are  appended  to  letter  symbols  to 
distinguish  measurements  of  the  same  kind  that 


are  made  on  features  differing  in  some  way,  such 
as  n1?  number  of  first-order  channels,  n2,  num¬ 
ber  of  second-order  channels,  and  so  forth.  The 
notation  inherited  by  geomorphologists  and  now 
employed  by  hydrologists  is  remarkably  poor  in 
the  use  of  subscripts.  The  principal  objection  is 
that  the  same  symbol,  for  instance  F,  is  used  for 
such  different  quantities  as  force  and  frequency, 
with  the  difference  indicated  by  a  subscript.  The 
writer  has  followed  a  conservative  practice  in  us¬ 
ing  subscripts  only  to  denote  the  order  of  the 
feature  upon  which  the  measurement  is  taken, 
such  as  n^,  number  of  first-order  channels,  and 
ll ,  total  length  of  first-order  channels. 

Channel  order  states  the  relative  magnitude  of 
a  segment  of  stream  channel  in  a  hierarchy  de¬ 
termined  by  arrangement  of  tributaries  with  re¬ 
spect  to  the  main  trunk.  The  system  of  ordering 
channel  segments  proposed  by  Horton  (1945)  and 
modified  by  Strahler  (1952b,  footnote,  p.  1120)  is 
followed  here.  The  smallest  permanent  channels, 
those  having  no  tributaries,  are  designated  first- 
order  channels.  At  the  junction  of  two  first-order 
channels,  a  second-order  channel  begins;  at  the 
junction  of  two  second-order  channels,  a  third- 
order  channel  begins,  and  so  forth.  A  channel  of 
any  given  order  can  directly  receive  tributaries 
of  all  lesser  orders.  The  trunk  stream,  or  master 
stream,  through  which  all  discharge  passes  has 
the  highest  order.  The  area  drained  by  a  channel 
of  order  u,  down  to  its  junction  with  another  chan¬ 
nel  of  order  equal  to  or  greater  than  u,  is  desig¬ 
nated  a  basin  of  order  u. 

Because  of  the  way  in  which  channel  ordering 
is  defined,  it  is  possible  to  state  the  minimum 
possible  number  of  channels  of  each  order  in  a 
basin  of  a  given  order.  Thus,  for  a  second-order 
channel,  at  least  two  first-order  channels  must 
be  present;  for  a  third-order  channel,  at  least 
two  second- order  channels  and  four  first-order 
channels  are  needed,  for  a  fourth-order  channel, 
two  third-order  channels,  four  second-order  chan¬ 
nels,  and  eight  first-order  channels  are  needed. 

In  general,  in  a  basin  of  order  uQ  the  minimum 
number  of  channel  segments  of  order  u,  is 

nu  =  2  u°‘u  (1) 

The  minimum  total  number  of  segments  of  all 
orders  is 

£nu=2]2Uo"1  ® 

U  =1  U=1 

Basin  area,  or  planimetric  area,  a,  is  the  area 
in  square  miles  of  the  outline  of  the  watershed  of 
a  channel  as  projected  onto  the  horizontal  map 
plane.  This  area  is  commonly  used  to  calculate 
the  total  amount  of  water  entering  the  system 
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Figure  1.  Inflection  angle  as  measured  on  various  forms  of  first-order  channels. 


during  a  rain.  Basin  area  is  a  fundamental  dimen¬ 
sional  property  in  the  majority  of  the  quantitative 
expressions  of  drainage-basin  characteristics. 

Basin  perimeter  is  the  length  along  the  major 
divide  as  projected  onto  the  horizontal  plane  of  the 
map.  The  symbol  used  here  for  the  perimeter  is  p; 
the  unit  of  measure  is  miles.  The  perimeter  is  a 
linear  measure  of  the  size  of  the  basin,  and  is 
largely  independent  of  the  texture  of  the  topography. 
It  is  used  in  conjunction  with  basin  area  to  give  a 
measure  of  the  departure  of  the  basin  from  a  true 
circle  (see  basin  circularity),  and  in  conjunction 
with  the  relief  to  give  a  measure  of  general  steep¬ 
ness  of  the  basin  (see  relative  relief). 

Basin  relief  is  the  vertical  distance  between  the 
mouth  of  the  basin  and  the  highest  point  in  the  ba¬ 
sin,  which  in  the  basins  in  this  study  is  always 
located  on  the  major  divide.  The  symbol  for  relief 
is  r;  the  unit  of  measurement  is  feet.  In  order  to 
find  the  relief,  it  is  usually  necessary  to  inter¬ 
polate  between  contour  lines  at  the  mouth  and  high¬ 
est  point.  This  is  done  to  the  nearest  one-fourth 
contour  interval,  since  it  is  desirable  to  avoid  the 
rigorous  but  clumsy  expression  of  r  given  with  a 
possible  range  of  two  contour  intervals. 

Total  length  of  channel  segments  of  each  order 
within  the  basin  were  measured  and  recorded  un¬ 
corrected  for  gradient,  in  units  of  miles.  The  sym¬ 
bol  used  for  the  total  length  of  channel  segments  or 


order  u  is  lu.*  The  mean  length  of  channel  seg¬ 
ments  of  order  u  is  lii,  where  n„  is  the  number  of 

nu 

segments  of  order  u,  and  is  represented  by  the 
symbol  lu.  The  total  length  of  channels  of  all  or¬ 
ders  in  a  basin  is  S  lu. 

Inflection  angle  of  contour  lines  is  the  angle 
which  a  contour  line  makes  with  itself  where  it  de¬ 
picts  a  channel  and  is  given  the  symbol  ip.  The  in¬ 
flection  angle  has  apparently  not  previously  been 
studied.  In  this  study,  the  minimum  inflection  angle 
of  each  first-order  stream  was  measured  from  a 
map  by  use  of  protractor  and  straightedge.  Figure 
1  shows  several  common  forms  of  first-order  chan¬ 
nels,  as  delineated  by  contour  lines,  and  the  inflec¬ 
tion  angle  for  each  is  illustrated. 

It  is  obvious  that  geometrically  the  gradient  of 
a  channel  can  vary  from  zero  up  to  the  angle  of 
slope  of  the  adjacent  valley  sides,  but  no  steeper. 
Increasing  difference  between  the  channel  gradient 
and  the  slope  of  the  adjacent  valley  sides  produces 
a  more  acute  inflection  angle.  Thus,  the  inflection 

♦Confusion  may  arise  because  of  the  identity  of  form  of 
the  letter  1  and  numeral  1.  In  this  notation  all  subscripts 
refer  to  order  number,  hence  a  subscript  1  always  means 
the  numeral  one.  The  letter  1  is  always  followed  by  a 
subscript  denoting  channel  order.  Therefore,  1  appearing 
without  a  subscript  always  means  one. 
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angles  of  first-order  channels  are  usually  greater 
than  of  higher  order  channels.  In  Figure  1  a,  the 
channel  gradient  is  nearly  the  same  as  the  slope  of 
the  valley  side,  and  gives  a  broadly  obtuse  inflec¬ 
tion  angle.  In  Figure  1  b,  the  gradient  is  some¬ 
what  less  than  the  slope  of  the  adjacent  valley 
sides,  and  the  inflection  angle  is  acute;  in  Figure 
1  c,  the  channel  gradient  is  much  less  than  the 
adjacent  valley-side  slopes,  and  the  inflection  an¬ 
gle  is  even  more  acute.  As  valleys  are  deepened 
the  angle  of  inflection,  as  shown  on  a  map  of  the 
channel,  may  be  expected  to  become  more  acute. 
Thus,  the  evolution  of  a  first-order  basin,  initiated 
by  a  small  gully  having  as  its  gradient  the  slope  of 
the  surface  upon  which  it  formed,  to  a  mature 
basin  form  in  which  channel  gradient  is  consider¬ 
ably  lower  than  the  slope  of  the  adjacent  valley 
sides,  might  be  traced  by  changes  in  the  inflection 
angle. 

New  tributaries  to  existing  channels  tend  to 
develop  perpendicular  to  the  contour  line,  so  as  to 
have  the  steepest  possible  slope.  The  projection  of 
the  new  tributary  on  a  map  makes  an  angle  with 
the  main  channel  equal  to  the  complement  of  the 
angle  which  the  contour  line  makes  with  the  main 
channel.  Main  channels  which  have  highly  acute 
inflection  angles  will  tend  to  develop  tributaries 
entering  at  nearly  right  angles.  The  relation  be¬ 
tween  angle  of  entry  (axil  angle*),  valley-side 
slope,  and  channel  gradient  was  found  by  Horton  to 
be 


where  £  is  the  axil  angle,  y  the  channel  gradient, 
and  0  the  valley-side  slope  (Horton,  1945,  p.  349). 
This  equation  can  be  derived  from  equation  19  in 
Appendix  I,  by  setting 


The  relations  among  channel  gradient,  valley- 
side  slope,  dihedral  angle  between  valley-sides, 
and  the  inflection  angle  can  be  expressed  in  geo¬ 
metric  laws,  and  the  derivation  of  these  laws  is 
given  in  Appendix  I. 

The  measurements  of  inflection  angles  of  first- 
order  channels  taken  in  this  study  are  biased  by 
the  habit  or  convention  of  indicating  small  streams 
adopted  by  the  individual  cartographer.  It  is  there¬ 
fore  advisable  in  future  studies  to  take  inflection- 
angle  measurements  on  small  channels  in  the  field. 
Reliable  results  from  map  measurements  might  be 
obtained  for  second  or  third- order  channels,  but 

*Axil  angle  has  been  used  by  botanists  for  150  years  to 
denote  the  upper  angle  between  the  main  stem  and  a 
branch.  Its  use  in  this  sense  is  analogous. 


no  more  than  a  general  indication  of  the  inflection 
angles  can  be  obtained  from  map  measurement  of 
first-order  channels. 


DERIVED  QUANTITIES 

Derived  quantities  are  those  that  must  be  ob¬ 
tained  by  calculation  from  the  properties  measured 
or  counted  directly  from  maps.  The  accuracy  of  the 
derived  quantities  depends  directly  on  the  accuracy 
of  measured  quantities,  and  their  frequency  dis¬ 
tributions  depend  on  the  distributions  of  all  the 
measured  items  involved  in  their  calculation.  A 
greater  amount  of  freedom  is  involved  in  their  es¬ 
timation.  For  these  reasons  it  is  desirable  to  dis¬ 
tinguish  between  measured  and  derived  quantities. 
Upper-case  Roman  letters  are  used  here  as  sym¬ 
bols  for  derived  quantities,  except  that  derived 
angular  quantities  are  represented  by  lower-case 
Greek  letters. 

Drainage  density,  D,  is  the  ratio  of  the  total 
length  of  channels  of  all  orders  in  a  basin  to  the 

y*  1 

area  of  the  basin,  - and  is  therefore  the  num- 

ber  of  miles  of  channel  maintained  by  a  square 
mile  of  drainage  area  (Horton,  1945,  p.  283).  Drain¬ 
age  density  is  a  numerical  expression  of  the  fine¬ 
ness  of  texture  of  the  basin,  but  because  a  particu¬ 
lar  value  is  usually  nearly  constant  over  an  area  of 
homogeneous  lithology,  climate,  and  stage  of  de¬ 
velopment,  it  is  independent  of  the  basin  size. 
Drainage  density  has  the  dimension  of  inverse  of 
length,  and  is  a  useful  indicator  of  scale  of  size 
of  planimetric  units  in  a  fluvially  controlled  land¬ 
scape. 

Density  of  first-order  channel  segments,  Di,  is 
the  total  length  of  first-order  channel  segments  in 
a  basin  of  some  larger  order,  divided  by  the  area 

of  the  entire  basin,  — 1  .  It  is  useful  in  determin- 

a 

ing  whether  a  drainage  density  higher  in  one  area 
than  in  another  is  due  to  a  greater  length  from 
headward  growth  of  the  first-order  channels. 

Channel  frequency,  F,  is  the  ratio  of  the  total 
number  of  channel  segments  of  all  orders  in  a 

y  n 

basin  to  the  area  of  the  basin,  - —  (Horton,  1945, 

cl 

p.  285).  It  is  a  supplementary  measure  of  the  fine¬ 
ness  of  texture  of  the  topography,  but  may  be  in¬ 
dependent  of  drainage  density  because  it  depends 
only  on  number  of  channel  segments,  not  on  the 
lengths. 

Frequency  of  first-order  channels,  Fi,  is  the 
total  number  of  first-order  channels  divided  by  the 

area  of  the  entire  basin,  ■  n  *— *  This  measure  is 

a 

useful  in  determining  whether  a  drainage  density 
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higher  in  one  area  than  in  another  has  resulted 
from  growth  of  new  channels,  or  simply  a  lengthen¬ 
ing  of  existing  channels. 

Ruggedness  number,  H,  is  the  product  of  drain- 
age  density  and  relief  of  basin,  both  in  the  same 

units:  H  =  or  H  =  (Strahler,  1954c). 

It  is  a  dimensionless  number,  since  D  has  the  di¬ 
mension  [l  1],  and  r  has  the  dimension  [L].  Areas 
of  low  absolute  relief,  but  high  drainage  density 
are  as  rugged,  when  their  small  size  is  considered, 
as  areas  with  great  relief  but  low  drainage  density. 
Areas  that  are  geometrically  similar  in  planime- 
tric  and  relief  features  will  have  the  same  value 
of  the  ruggedness  number. 

Relative  relief,  R,  which  is  a  measure  of  the 
general  steepness  of  the  basin  from  summit  to 
mouth,  can  be  measured  in  a  variety  of  ways. 
Schumm  (1954,  p.  22)  defined  relief  ratio  as  “the 
ratio  between  the  total  relief  of  a  basin.  .  .and 
the  longest  dimension  of  the  basin  parallel  to  the 
principal  drainage  line.  This  relief  ratio  is  a  di¬ 
mensionless  height-length  ratio  and  is  equal  to  the 
tangent  of  the  angle  formed  by  two  planes  inter¬ 
secting  at  the  mouth  of  the  basin;  one  representing 
the  horizontal,  the  other  passing  through  the  high¬ 
est  point  of  the  basin.”  In  oddly  shaped  basins,  or 
basins  with  the  highest  point  at  some  place  other 
than  at  the  position  farthest  removed  from  the 
mouth,  there  is  a  question  as  to  what  should  be 
taken  as  the  length  measurement,  and  the  interpre¬ 
tation  of  the  relief  ratio  as  the  tangent  of  the  angle 
between  two  planes  is  not  appropriate. 

For  this  study  relative  relief,  R,  is  defined  as 
the  ratio  of  the  basin  relief,  expressed  in  units  of 

miles,  to  the  length  of  perimeter:  R  =  The 

advantage  of  this  measure  over  relief  ratio  is  that 
both  relief  and  perimeter  had  been  measured  for 
other  purposes;  the  position  of  the  highest  point  of 
the  basin  with  respect  to  the  mouth  does  not  affect 
the  measure.  A  rational  relation  between  the  meas¬ 
ure  of  relief  used  here  and  Schumm’s  relief  ratio 
can  be  developed  by  assuming  various  shapes  for 
the  basin. 

Ratio  of  channel  lengths  to  basin  perimeter,  P, 
is  a  dimensionless  measure  of  fineness  of  topo¬ 
graphic  texture  relative  to  basin  size.  Two  basins 
of  different  sizes,  but  with  geometrically  similar 
planimetric  properties,  will  have  different  values 
of  drainage  density,  but  identical  values  of  the 
channel-perimeter  ratio. 

Basin  circularity,  C,  measures  the  degree  to 
which  a  basin  approximates  a  circle  in  planimetric 
outline.  Circularity  is  defined  as  the  area  of  the 
basin  divided  by  the  area  of  the  circle  with  equal 
perimeter  (Miller,  1953,  p.  8).  The  circle  can  be 


defined  uniquely  as  that  plane  figure  in  which  the 
greatest  area  is  bound  by  the  shortest  possible 
perimeter.  Therefore,  if  the  area  of  any  irregular 
figure,  such  as  a  drainage-basin  outline,  be  divided 
by  the  area  of  the  circle  having  the  same  perimeter, 
the  quotient  will  be  a  number  between  zero  and  one. 
The  measure  is  ambiguous,  though,  as  it  does  not 
discriminate  between  departures  from  perfect  cir¬ 
cularity  due  to  irregularities  in  gross  outline,  or 
to  crenulations  in  the  perimeter. 

Bifurcation  ratio,  BU:u  +  l  =  — ^ — ,  is  the  ratio 
-  nu+i 

between  the  number  of  channel  segments  of  a  given 

order  and  the  number  of  channel  segments  of  the 

next  higher  order.  According  to  Horton’s  law  of 

stream  numbers,  values  of  B  tend  to  be  the  same 

for  all  values  of  u  within  a  given  basin  (Horton,  1945, 

p.  286).  The  value  of  B  characteristic  of  all  orders 

within  the  basin  is  the  slope  of  the  regression  line 

when  numbers  of  channel  segments  are  plotted  on 

semi-log  paper  against  order  number  (Strahler, 

1954b,  p.  345). 

Length  ratio,  L„.u-i  =  is  the  ratio  of  the 

lu-1 

mean  lengths  of  channels  of  a  given  order,  to  the 
mean  length  of  channels  of  the  next  lower  order. 
According  to  Horton’s  law  of  stream  lengths,  the 
mean  channel  length  for  each  order  increases  as 
an  exponential  function  of  channel  order  (Horton, 

1945,  p.  286).  Whether  the  values  of  Lu;u.i  for  a 
basin  are  constant  within  a  basin,  or  change  sys¬ 
tematically  with  order  is  a  question  that  must  be 
answered  by  statistical  analysis. 

Maximum  valley-side  slope,  0,  is  the  slope  in 
degrees  of  the  valley  side  at  its  steepest  point  along 
an  orthogonal  profile.  This  definition  assumes  that 
short,  oversteepened  sections  of  the  valley  side 
near  channels,  that  are  the  obvious  result  of  recent 
rejuvenation,  are  not  included  in  the  profile.  It  is 
not  likely  that  such  sections  of  slope  would  be  in¬ 
cluded  in  map  measurements,  but  it  is  necessary 
to  exclude  them  from  the  definition,  as  valley-side 
slopes  are  also  measured  in  the  field  (see  maxi¬ 
mum  valley-side  slopes  under  Measured  Properties). 
The  measure  is  obtained  from  a  map  by  computing 
the  horizontal  distance  between  two  contour  lines 
that  are  separated  vertically  by  say,  100  feet,  and 
converting  the  tangent  to  the  angle  in  degrees.  De¬ 
pending  on  the  variability  of  slopes  within  the  basin, 
15  to  75  slope  readings  are  needed  to  give  an  ade¬ 
quate  estimation  of  the  mean  value. 

The  importance  of  the  valley-side  slopes  in  stud¬ 
ies  of  erosion  lies  in  the  fact  that  the  energy  avail¬ 
able  for  erosion  in  foot-pounds  per  pound  of  water 
per  linear  foot  is  equal  to  the  sine  of  the  slope. 

Dihedral  angle  between  valley  sides,  £,  meas¬ 
ured  in  degrees,  is  the  angle  between  the  valley 
sides  measured  perpendicular  to  the  stream.  It 
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cannot  be  measured  easily  either  in  the  field  or 
from  the  map,  but  it  can  be  calculated.  Figure  2 
is  an  illustration  of  the  dihedral  angle,  inflection 
angle,  gradient,  and  valley-side  slope.  Appendix  I 
gives  the  geometrical  relation  among  the  quantities 
Q,&,  C,  and  y. 

Gradient  of  stream  channel,  y,  expressed  in  de¬ 
grees,  can  be  obtained  from  the  map  by  measuring 
the  distance  along  the  channel  between  two  or  more 
contour  lines.  This  quantity  is  of  great  importance 
in  the  equilibrium  of  the  basin,  since  it  is  related 
to  rates  of  transport,  runoff,  and  celerity  of  flood 
crests. 


Figure  2.  Relations  among  valley-side  slope,  inflection 
angle,  channel  gradient,  and  dihedral  angle  between 
valley  sides  for  an  idealized  channel. 


Problems  of 


METHODS  OF  FIELD  OBSERVATION 

Field  investigations  were  carried  out  (1)  to  ob¬ 
serve  the  accuracy  of  the  map,  particularly  in  re¬ 
gard  to  the  extent  of  the  drainage  net,  and  (2)  to 
obtain  measurements  of  characteristics  of  the  soil 
and  vegetal  cover  that  are  believed  to  influence 
morphometric  properties  of  the  basin. 

At  representative  points  within  each  basin 
several  stations  were  established  at  which  a  ser¬ 
ies  of  measurements  and  observations  was  made. 

In  large  basins  about  10  stations  per  square  mile 
were  established;  in  small  basins  the  density  in¬ 
creased  to  200  stations  per  square  mile.  To  esti¬ 
mate  the  number  of  measurements  required  to  de¬ 
tect  a  difference  of  a  specified  amount  in  some 
property  of  interest,  among  a  number  of  areas, 
a  prior  knowledge  of  the  standard  deviation  of  the 
property  is  necessary.  Estimates  of  the  standard 
deviations  of  the  properties  studied  in  the  field, 
excepting  maximum  valley-side  slopes,  were  not 
available  at  the  beginning  of  the  study.  Without 
such  estimates  it  is  impossible  to  judge  in  ad¬ 
vance  the  necessary  number  of  readings  within 
each  basin,  or  the  number  of  basins,  needed  to  de¬ 
tect  differences  in  means  of  the  properties  studied. 
It  was  assumed,  however,  that  for  a  given  property, 
greater  variability  would  exist  among  the  means 
for  all  basins  than  among  readings  within  the 
basins;  analysis  of  the  data  proves  this  assumption 
to  be  correct  (Table  6).  Therefore,  the  emphasis 
was  placed  on  gathering  data  from  a  large  number 
of  basins,  and  fewer  measurements  were  taken 
within  each  basin. 

A  major  problem  in  the  field  study  of  each  area 
was  to  decide  which  linear  depressions  should  be 
included  as  channels,  and  which  should  be  omitted. 
The  crux  of  the  problem  is  in  the  identification  of 
the  smallest  “finger-tip,”  or  first-order  channels, 
because  the  channels  of  higher  order  present  no 
difficulty  once  all  first-order  channels  are  drawn 
on  the  map.  First-order  channels  can  be  defined 
as  the  smallest  permanent  depressions  found  in  a 
connected  drainage  system,  whose  form  is  main¬ 
tained  indefinitely  against  agents  tending  to  obliter¬ 
ate  it,  by  the  action  of  water  flowing  in  a  concen¬ 
trated  stream.  However,  any  definition  depending 
on  the  knowledge  of  the  behavior  of  the  depression 
through  several  years,  or  seasons,  is  unusable  in 
studies  in  which  observations  are  made  in  each 
basin  during  a  single,  relatively  brief  period.  For 


Field  Study 

the  purposes  of  field  identification,  true,  functional, 
erosional,  first-order  channels  are  defined  as  linear 
depressions  having  the  following  characteristics: 

(1)  Water  entering  the  depressions  is  forced  to 
flow  in  a  concentrated  stream,  and  under  no 
conditions  can  become  sheetflow. 

(2)  At  its  lower  end  the  depression  (a)  joins 
another  similar  depression  of  equal  or 
greater  size,  (b)  disperses  onto  an  alluvial 
fan  or  cone,  or  (c)  ends  at  a  sinkhole  or 
other  localized  area  of  recognizably  great 
permeability. 

(3)  Erosion  and/or  creep  have  produced  (a) 
graded  valley  sides  that  lead  down  to  the 
depression,  and  have  orientations  noticeably 
different  from  the  hillside  adjoining  the 
area  drained  by  the  depression.  The  gradient 
of  the  depression  in  this  case  is  lower  than 
the  slope  of  the  adjacent  surfaces  (Gilbert, 
1877,  p.  115).  Or  (b)  the  depression  is  steep¬ 
sided  with  rectangular  or  trapezoidal  cross- 
section,  its  depth  being  on  the  order  of  ten 
times  the  maximum  micro-relief  on  the 
older  unmodified  surfaces  near  the  depres¬ 
sion. 

This  definition  excludes  steep-gradient,  grassy 
swales  through  which  water  might,  after  a  mild 
rain,  pass  in  the  form  of  sheetflow.  However,  per¬ 
haps  such  channels  need  not  be  excluded  from  con¬ 
sideration  as  higher  order  channels,  uncommon 
though  they  may  be. 


OBSERVED  WATERSHED  PROPERTIES 

Vegetation:  The  typical  vegetation  assemblages 
with  their  elevations  and  exposures  were  noted  in 
each  basin  (Table  3).  Identifications  of  species 
were  made  from  Southwestern  Trees  (Little,  1950). 
Certain  species  of  trees  and  bushes  provide  a 
basis  from  which  inferences  about  the  climate  of 
the  basin  can  be  drawn  (Little,  1950;  Costello,  19?). 
Estimation  of  the  precipitation-effectiveness  index 
for  each  basin  was  obtained  from  a  regression  line 
for  precipitation-effectiveness  against  elevation 
and  confirmed  by  field  data.  For  instance,  the 
saguaro  cactus,  creosote  bush,  and  paloverde  of  the 
desert  or  Lower  Sonoran  Life  Zone  indicate  a  pre¬ 
cipitation-effectiveness  index  of  less  than  15;  the 
shrub  live  oak  and  cliffrose  of  the  chaparral  or 
Upper  Sonoran  Life  Zone  indicates  a  P-E  index  of 
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35  to  40,  and  the  Engelmann  spruce  and  alpine  fir  of 
the  Hudsonian  Life  Zone  indicate  a  P-E  index  of  90 
to  120  (see  Climatic  Data). 

Lithology:  Basins  were  selected  for  study  partly 
on  the  basis  of  similarity  of  lithology  throughout.  In 
a  few  instances,  basins  were  found  to  have  areas  of 
rock  differing  from  that  of  the  major  part  of  the 
basin.  In  addition,  in  four  basins  the  rocks  lie  in 
near-horizontal,  alternating  weak  and  strong  layers, 
consisting  of  sandstone,  limestone,  and  shale,  or 
lava  and  pumice.  Descriptions  of  the  lithology  pre¬ 
vailing  in  each  basin  are  given  in  the  general  de¬ 
scription  of  basins.  If  consistent  structural  planes 
are  present  throughout  the  area,  their  strike  and 
dip  are  noted  also. 

Percentage  of  basin  lying  within  each  type  of 
cover:  Each  basin  studied  in  the  field  was  divided 
into  portions  covered  by  forest  and  bushes,  grass¬ 
land,  and  bare  soil  or  rock.  The  percent  of  total 
area  lying  within  each  type  of  cover  was  then  esti¬ 
mated,  either  from  aerial  photographs,  or  from  a 
vantage  point  providing  a  general  view  of  the  entire 
basin.  It  was  found  that  on  those  maps  with  green 
overprint  representing  vegetative  cover,  there  was 
often  little  correspondence  with  the  actual  cover. 
Data  are  given  in  the  general  description  of  each 
basin. 

Range  evaluation:  Many  of  the  basins  studied  in 
the  field  were  being  grazed  at  the  time  the  field 
work  was  done,  or  had  recently  been  grazed.  The 
extent  of  grazing  was  estimated  from  observation 
of  the  general  length  of  grass  stubble.  A  simple 
count  of  animal  density  is  inadequate;  for  instance, 
near  Tucson,  three  animals  per  square  mile  result 
in  considerably  over-grazed  range,  whereas  in  the 
Colorado  Front  Range,  many  more  animals  can  be 
supported  without  serious  over-grazing.  Data  for 
each  basin  are  given  in  the  general  description  of 
basins. 

Presence  of  gullies:  It  is  assumed  here  that  gul¬ 
lies  usually  signify  readjustment  of  the  drainage 
pattern  to  a  changed  environment.  Where  numerous 
gullies  were  present,  excluding  those  few  caused 
by  some  factor  not  affecting  the  entire  basin,  note 
was  made  of  whether  they  were  primarily  on  the 
valley  slopes  or  in  the  channels.  In  a  single  case 
extensive  gullying  was  found  in  both  situations.  It 
seems  likely  that  channel  trenching  and  gullying  on 
the  upland  slopes  may  be  due  to  opposite  causes. 

The  surface  characteristics  of  the  mantle  of  the 
basins  may  have  changed  greatly  with  the  advent  of 
human  activity.  Grazing,  logging,  burning,  and  min¬ 
ing  activities  leave  very  apparent  effects.  Grazing 
can  drastically  reduce  the  amount  of  cover,  and 
even  moderate  grazing  lowers  the  infiltration  capac¬ 
ity  of  the  soil.  Watershed  management  studies  in  the 
Colorado  Front  Range,  on  granitic  soils  similar  to 
the  majority  of  those  studied  here,  and  in  ponderosa 
pine  forest,  show  that  heavy  grazing  more  than 


doubles  the  amount  of  erosion  per  unit  area,  for 
each  inch  of  runoff  (Love,  1953,  p.  214).  It  could 
therefore  be  contended  that  measurement  of  surface 
characteristics,  including  vegetal  cover,  as  they 
are  now  will  not  reflect  the  condition  of  the  basin 
when  it  was  developed,  and  are  therefore  valueless. 
All  evidence  available  to  the  writer  indicates  that 
increase  of  drainage  density,  brought  about  by 
lowered  resistance  to  erosion,  concentration  of 
runoff,  decreased  infiltration  rate,  disruption  of 
plant  cover,  etc.,  occur  very  soon  after  the  change 
of  conditions.  New  gullies  on  slopes  below  a 
new  road,  or  in  areas  of  recently-burned  forest 
are  convincing  evidence  of  the  rapidity  of  the 
change.  Gullies  in  slopes  are  included  in  the  drain¬ 
age  net  and  it  is  therefore  considered  likely  that 
measurements  of  surface  properties  do  bear  rela¬ 
tion  to  the  extent  of  channels.  Part  of  the  difficulty 
is  in  the  assumption  inherent  in  the  objection.  Ba¬ 
sins  were  not  developed  then  left  unchanging  and 
static  for  the  next  surface  geologist  to  study,  but 
must  reflect  by  continual  development  the  changes 
in  climate  and  environment  within  the  basin. 

Perhaps  there  is  more  reason  to  believe  that 
degree  of  slope  will  not  be  influenced  immediately 
by  changes  in  surface  properties,  as  time  is  re¬ 
quired  to  accomplish  the  regrading  of  a  slope. 
Erosion  of  a  large  amount  of  material  is  necessary 
to  change  the  average  basin  slope  even  a  degree 
or  two.  Since  rates  of  creep  are  not  well-known, 
supporting  evidence  for  this  view  is  needed. 


GENERAL  PLAN  OF  FIELD  MEASUREMENTS 

At  each  station  established  in  a  basin,  a  series 
of  measurements  was  repeated  in  as  nearly  the 
same  manner  as  possible  from  basin  to  basin.  All 
measurements  should  have  been  taken  at  randomly 
located  points,  but  unfortunately,  it  was  impossible 
to  do  this  consistently.  In  cobbly  and  bouldery  soils 
it  was  often  difficult  to  locate  places  on  which 
several  of  the  measurements  could  be  made.  In 
other  areas,  such  as  the  Chinle-shale  badlands, 
where  full  randomization  would  be  relatively  easy, 
it  was  less  important  due  to  the  high  degree  of 
homogeneity  of  the  area. 


MEASURED  PROPERTIES 

Infiltration  capacity:  The  present  study  did  not 
include  a  general  investigation  into  the  factors  in¬ 
fluencing  the  infiltration  capacity  of  the  soil.  The 
infiltration  capacity  measurements  were  taken 
solely  with  the  intention  of  obtaining  comparative 
data,  not  necessarily  the  actual  rate  of  infiltration 
during  a  rain. 

The  infiltrometer  used  in  this  study  was  designed 
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by  the  writer  to  be  portable,  light  in  weight,  and 
conservative  in  use  of  water.  These  requirements 
were  important  because  in  many  areas  it  was  nec¬ 
essary  to  carry  the  apparatus  and  water  over 
rugged  terrain. 

The  principle  of  measurement  is  simple:  water 
from  a  tank  is  drawn  through  a  metering  device, 
an  adjustable  valve,  a  sprinkler,  and  distributed  by 
inspection  evenly  over  a  plat  of  soil  one  square  foot 
in  area  at  the  maximum  rate  of  infiltration.  Dis¬ 
charge  is  measured  by  the  meter  in  ounces  per 
minute,  and  is  converted  to  inches  (depth)  per  hour 
over  an  area  of  one  square  foot.  The  apparatus  is 
mounted  on  a  standard  surveying  tripod,  and  weighs 
45  pounds  when  full  of  water  (Figure  3).  The  reser¬ 
voir  tank  holds  2  1/2  gallons,  and  will  supply 
enough  water  for  two  to  five  tests,  depending  on  the 
length  of  run,  and  the  infiltration  capacity  of  the 
soil. 


Figure  3.  Infiltrometer  assembly  diagram. 


The  one  square-foot  area  is  bounded  by  a  stain¬ 
less  steel  ring  with  a  sharpened  lower  edge,  that 
can  be  driven  into  the  ground.  Doubtless  a  double¬ 


ring  arrangement,  with  water  filling  an  outside  ring 
to  lessen  the  amount  of  lateral  flow,  would  be  de¬ 
sirable  were  accuracy  the  main  consideration.  The 
added  weight  of  another  ring  and  the  necessarily 
greater  quantity  of  water,  would  reduce  the  distance 
which  the  instrument  could  be  carried  during  a  day’s 
work.  The  infiltrometer  used  in  this  study  is  defin¬ 
itely  a  reconnaissance  instrument,  and  cannot  com¬ 
pare  in  accuracy  with  more  elaborate  equipment. 

The  tank  containing  the  water  is  made  of  stain¬ 
less  steel,  and  is  built  on  the  principle  of  the  Mari- 
otte  bottle  (Westphal,  1952,  Zweiter  Teil,  p.  15). 

The  tank  can  be  closed  air  tight  by  a  rubber  stopper; 
the  only  air  entering  the  tank  comes  through  a 
copper  tube  extending  down  almost  to  the  bottom  of 
the  tank.  Thus,  the  only  free  surface  is  at  the  bot¬ 
tom  of  the  tube,  and  once  the  air  has  started  bubbling 
out  of  the  tube  into  the  tank,  the  pressure-head  is 
measured  from  the  end  of  the  tube.  Therefore, 
throughout  the  course  of  the  test,  though  the  level 
of  the  water  decreases,  the  hydrostatic  head  re¬ 
mains  the  same.  This  is  necessary  because  the 
measuring  instrument  is  accurate  under  constant 
head.  Once  a  definite  rate  of  flow  has  been  estab¬ 
lished,  it  is  not  necessary  to  make  repeated  changes 
in  the  valve  setting,  as  would  be  the  case  were  the 
head  continuously  decreasing. 

The  measuring  device  is  a  Flowrator,*  which 
is  basically  a  tapered  glass  tube  with  a  metal 
float  fitting  closely  within  the  tube.  For  protec¬ 
tion,  a  pressed  steel  cover  surrounds  the  tube 
and  supports  the  pipe  connections.  Water  flows  up¬ 
ward  through  the  tube;  as  the  discharge  increases 
the  frictional  force  of  the  water  on  the  float  in¬ 
creases,  and  an  upward  acceleration  is  imparted. 
Higher  in  the  tube,  the  diameter  of  the  opening 
and  the  area  of  discharge  are  greater;  the  vel¬ 
ocity  of  water  flowing  around  the  float  is  less,  and 
hence  the  frictional  force  is  less.  The  float  there¬ 
fore  seeks  a  level  at  which  its  weight  is  exactly 
balanced  by  the  upward  force  of  the  moving  water, 
and  this  position  is  a  measure  of  the  discharge. 

The  maximum  discharge  possible  is  11.2  ounces 
per  minute  or  8.4  inches  per  hour.  Figure  4  gives 
the  conversion  graph  from  percent  of  total  flow  to 
inches  per  hour  over  one  square  foot.  After  leav¬ 
ing  the  Flowrator,  the  water  passes  down  through 
a  sprinkler  and  is  distributed  evenly  over  the  area 
within  the  ring. 

The  surface  of  the  test  area  is  prepared  by  re¬ 
moving  loose  leaves,  twigs,  accumulated  pine  duff 
and  needles,  leaving  only  grass  on  the  largely  min¬ 
eral  soil.  At  the  beginning  of  a  test,  the  discharge 
is  set  at  75%  of  the  maximum,  8.4  ounces  per  min¬ 
ute,  or  6.3  inches  per  hour.  If  the  soil  is  dry,  it 
will  take  15  to  180  seconds  to  saturate  the  surface 

♦Manufactured  by  Fisher  and  Porter  Company,  Hatboro, 
Pennsylvania. 
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Figure  4.  Conversion  graph:  percent  of  total  flow  (100%  = 
11.2  ounces  per  minute)  to  inches  per  hour  over  an  area 
of  one  square  foot. 


and  to  form  excess.*  The  time  at  which  the  excess 
forms  is  recorded,  and  the  discharge  is  decreased. 
Some  practice  is  required  to  adjust  the  discharge 
to  match  the  infiltration  rate,  and  at  best  it  takes 
one  minute  to  determine  if  the  machine  is  adjusted 
correctly.  On  most  poorly-developed,  gravelly 
soils,  during  the  first  15  minutes  of  the  test,  the 
rate  of  infiltration  decreases  rather  rapidly. 
Thereafter,  the  rate  is  fairly  constant,  though  it 
may  decrease  slowly  if  considerable  clay  is 
present  in  the  soil.  It  is  usually  possible  to  adjust 
the  instrument  to  within  1%  of  flow  after  20  to  25 
minutes  have  elapsed,  giving  an  adjustment  error 
on  0.084  inches  per  hour,  and  the  rate  of  infiltra¬ 
tion  usually  changes  little  thereafter.  The  “instan¬ 
taneous”  rates  for  3,  5,  10,  15,  20,  25,  30  and 
sometimes  40  and  50  minutes  were  recorded.  The 
final,  stable  rate  (usually  the  minimum  rate)  was 
taken  as  the  infiltration  capacity  (Horton,  1942, 
p.  480).  The  values  of  infiltration  capacity  obtained 
with  this  instrument  are  probably  higher  than  the 
actual  value  prevailing  during  a  hard,  general  rain, 
since  no  allowance  is  made  for  the  lateral  flow  of 
water  when  it  passes  below  the  edge  of  the  ring. 


*L.  D.  Love  (personal  communication,  1955)  and  others 
believe  that  because  of  wetting  difficulties,  most  reliable 
results  may  be  obtained  on  pre-wetted  plats  that  have 
drained  12  to  24  hours.  The  writer  encountered  wetting 
difficulties  commonly  in  soils  with  a  deep  pine-needle 
deposit.  In  a  few  cases  a  small  amount  of  detergent  was 
added  to  the  water,  and  this  usually  shortened  the  time 
needed  to  saturate  the  surface.  Too  few  trials  were  made 
to  determine  whether  this  altered  the  final  infiltration 
capacity. 


A  greater  source  of  error,  however,  may  be  the 
disturbance  of  the  soil  structure  caused  by  pound¬ 
ing  in  the  ring,  thereby  opening  up  fissures  not 
otherwise  present.  The  values  obtained  are  believed 
to  be  consistent  within  themselves,  since  values  of 
infiltration  capacity  vary  greatly  between  areas  of 
different  rock  type,  and  less  so  within  areas  of 
similar  type. 

Soil  strength,  s,  is  measured  under  instantan¬ 
eous  impact  of  a  12  pound  iron  shot,  dropped  a 
standardized  distance,  to  impart  60  foot  pounds  of 
energy  normal  to  the  ground  surface.  This  distance 
is  5  feet  for  horizontal  surfaces,  and  5  feet  divided 
by  the  cosine  of  the  slope  on  sloping  surfaces.  Sat¬ 
isfactory  results  can  be  obtained  on  slopes  up  to 
35  degrees.  The  volume  of  the  imprint  of  the  shot 
divided  into  60  foot-pounds  gives  a  measure  of 
strength  in  pounds  per  square  foot. 


60  ft. -lbs. 
x  ft.3 


DU  ,,  2 

=  —  lbs.  per  ft. 

x 


The  volume  of  the  imprint  can  be  obtained  from  a 
graph  relating  diameter  of  the  base  of  a  spherical 
segment  to  volume.  The  diameter  of  the  imprint  is 
easily  measured.  Appendix  III  gives  the  derivation 
of  the  relation  between  the  diameter  of  imprint  and 
volume  of  imprint.  Figure  5  illustrates  the  method; 
Figure  6  is  a  graph  relating  diameter  of  imprint 
directly  to  strength  of  soil. 


Figure  5.  Principle  of  measuring  soil  strength  by  drop¬ 
ping  a  spherical  weight  a  standardized  distance  and 
measuring  the  volume  of  imprint. 
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Figure  6.  Conversion  of  diameter  of  imprint  of  12-lb. 
shot,  diameter  =  0.378  ft.,  to  soil  strength,  whendropped 
the  distance  necessary  to  impart  60  ft. -lbs.  normal  to 
the  soil  surface. 


To  standardize  the  moisture  conditions  under 
which  the  strength  tests  were  made,  two  or  three 
repetitions  of  the  test  were  made  on  the  soaked 
plat  of  ground  used  in  the  infiltration  test  to  obtain 
a  value  of  the  strength  of  the  soil  approximating 
its  condition  after  a  hard  rain.  In  addition,  three 
or  four  tests  were  made  in  the  surrounding  area 
on  the  soil  in  its  natural  condition.  In  most  areas 
the  soil  was  extremely  dry;  only  in  a  few  cases 
was  it  damp. 

It  would  have  been  desirable  to  randomize  all 
measurements  of  strength.  However,  the  measure¬ 
ment  of  wet  soil  strength  was  necessarily  limited 
to  the  area  selected  for  infiltration  measurements, 
which  were  in  turn  limited  by  the  stoniness  of  the 
ground.  Randomization  over  dry  surfaces  would 
produce  tests  on  cobbles,  logs,  etc.,  and  the  re¬ 
sults  would  not  be  indicative  of  soil  strength.  The 
measurements  were  probably  not  systematically 
biased,  as  there  was  no  way  of  estimating  the 
strength  of  the  point  directly  under  the  shot  prior 
to  making  the  test. 

Size-distribution  of  soil  was  described  at  each 
station  within  the  infiltrometer  ring  by  the  visual 
and  manual  methods  developed  by  the  Department 
of  Soil  Mechanics  at  Columbia  University  (Bur- 
mister,  1952,  p.  48).  The  average  soil  descrip¬ 
tion  for  test  plats  within  each  basin  is  given  in 
Table  4. 

Percent  bare  area  was  obtained  by  measuring 
the  exposed,  bare,  mineral  soil  with  a  steel  meas¬ 
uring  tape  extended  for  50  feet  over  the  ground  in 
a  randomly  oriented  direction.  The  ground  at  each 
foot  mark  on  the  tape  was  inspected,  and  if  bare  of 
all  cover,  such  as  grass,  fallen  leaves,  twigs, 
limbs,  cobbles,  or  larger  rock  fragments,  it  was 


given  a  count  of  one.  For  each  foot  mark  lying 
above  a  section  of  ground  that  was  covered  in  any 
way,  the  count  was  zero.  Upon  counting  through  50 
feet,  the  number  of  foot  marks  lying  above  bare 
ground  was  multiplied  by  two  to  give  an  expression 
in  percent  of  total  area. 

The  direction  of  traverse  and  point  of  origin 
were  randomized  by  throwing  a  geologic  hammer  in 
the  air  with  a  spinning  motion,  and  extending  the 
tape  in  the  direction  in  which  the  pick  pointed  upon 
landing.  Two  such  readings  were  made  at  each  sta¬ 
tion.  The  average  of  the  results  for  each  basin  are 
given  in  Table  3. 

Roughness,  or  number  and  size  of  rock  fragments 
lying  loose  on  the  surface,  or  partially  buried,  was 
found  by  counting  and  measuring  all  fragments  1/2 
inch  or  greater  in  mid-diameter  within  a  circle  with 
a  radius  of  one  foot.  The  center  of  the  circle  was 
located  by  selecting  a  number  of  feet  between  9  and 
100  from  a  table  of  random  numbers;  an  orientation 
was  found  by  throwing  a  geologic  hammer  into  the 
air  in  the  same  manner  as  described  above  for  per¬ 
cent  bare  area.  Fragments  were  measured  to  the 
nearest  1/16  inch  by  small  calipers.  The  distance 
which  a  buried  fragment  extended  above  the  ground 
was  recorded.  One  circle  was  measured  at  each 
station.  The  roughness  number,  M,  is  the  average 
total  length  of  diameters  of  pebbles  in  each  circle 
for  the  entire  basin.  M  is  then 

N  N 

fimi  J2 
i  =  1 

N  N 

where  fi  is  the  number  of  rock  fragments  in  the  i’th 
circle  of  one  foot  radius,  im  is  the  average  mid¬ 
diameter  of  rock  fragments  in  the  i’th  circle,  x  i  • 
is  the  diameter  of  the  j’th  pebble  in  the  i’th  circle, 
and  N  is  the  number  of  circles  measured  in  each 
basin.  The  roughness  number  increases  rapidly 
with  increasing  number  of  loose  fragments,  and 
also  with  increasing  average  size  of  rock  fragments. 
Among  the  basins  studied,  those  having  a  large 
number  of  smaller  fragments  indicate  greater 
roughness  than  basins  with  a  small  number  of  lar¬ 
ger  fragments. 

Maximum  valley-side  slopes,  0,  were  measured 
in  the  field  with  an  Abney  hand  level,  calibrated  to 
read  to  1/2°.  Within  each  basin,  about  40  measure¬ 
ments  were  usually  taken,  which  would  provide  suffi¬ 
cient  probability  of  detecting  differences  of  1  1/2°  or 
2°  in  the  means  of  two  basins  when  the  standard  devi¬ 
ation  is  about  4°.  The  mean  slopes  of  all  basins  for 
which  data  were  taken  are  given  in  Table  2. 

In  measuring  slopes,  the  writer  walked  along  the 
valley  side  near  the  upper  limit  of  the  steepest 
portion  of  the  slope,  with  an  assistant  near  the  foot 
of  the  slope.  At  fairly  regular  intervals,  slope  and 
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orientation  readings  were  recorded.  In  basins  in 
which  trenching  or  rejuvenation  of  the  channel  was 
apparent,  the  short,  steep  slope  near  the  channel 
was  ignored,  and  the  long  slope  on  the  main  valley 
side  was  measured. 

In  large  basins,  those  of  one  square  mile  or 
greater  area,  a  system  of  randomized  sampling 


that  will  still  allow  rapid  measurement  is  needed. 
In  such  cases  a  tendency  to  ignore  slopes  of  lower 
angle,  and  take  too  many  readings  on  steep 
slopes  would  bias  the  mean.  In  smaller  basins, 
where  several  readings  are  taken  on  each  valley 
side,  bias  is  less  likely  because  a  greater  percent 
of  the  area  is  sampled. 


Climatic  Data 


.  .  .  little  essential  advance  can  be  made  in  the  investigation  of  rainfall  by  continuation  of  readings  at  arbitrar¬ 
ily  selected  localities  which  happen  to  be  centers  of  population.  Intensive  study  of  rainfall  for  a  relatively  short 
period  at  carefully  selected  critical  localities  would  advance  our  knowledge  far  more  than  the  continuation  for 
many  years  of  routine  readings  at  sporadically  located  stations.  (Turnage  and  Mallery,  1941,  from  Introduction, 
by  Forest  Shreve) . 


Since  1941  there  has  been  no  advance  in  system¬ 
atic  collection  of  weather  records,  especially  tem¬ 
perature  data,  in  mountainous  regions.  In  compiling 
climatic  data  for  this  study,  considerable  reliance 
was  placed  on  records  (some  of  only  a  few  years 
duration)  obtained  from  high  mining  camps  in  the 
Colorado  Front  Range  (1860  to  1915),  from  ranches 
in  New  Mexico  and  Arizona,  and  from  Forest  Service 
ranger  stations  and  fire  lookouts.  Perforce,  much 
of  the  data  is  from  the  years  before  1930,  and 
obtained  from  Bulletin  W,  Climatic  Summary  of  the 
United  States  (U.  S.  Weather  Bureau,  Department 
of  Agriculture,  1933).  Although  the  Weather  Bureau 
under  the  Department  of  Commerce  has  placed  em¬ 
phasis  on  improving  techniques  of  forecasting, 
there  is  a  great  need  for  research  into  the  effects 
of  topography,  exposure,  and  elevation  on  precipi¬ 
tation  amounts  and  intensities,  temperature,  and 
snow  accumulation,  in  connection  with  hydrologic 
studies  for  flood  control,  erosion  problems,  forest 
and  range  management,  and  other  conservation 
programs.  This  need  is  being  met  partially  by 
utilization  of  weather  records  of  Forest  Service 
Experimental  Stations,  and  by  regional  studies 
such  as  those  prepared  by  the  Hydrologic  Services 
Division,  U.  S.  Weather  Bureau  (U.  S.  Weather 
Bureau  Technical  Paper  No.  24,  1954,  Parts  I  and 
II).  But  these  summaries  would  be  greatly  im¬ 
proved  were  adequate  mountain  weather  data  avail¬ 
able. 

The  data  secured  consist  of  average  monthly 
precipitation,  monthly  average  of  mean  daily  tem¬ 
perature,  average  number  of  days  per  year  with 
.01  inches  or  more  precipitation,  and  elevation  of 
the  station.  In  addition,  the  depth  of  precipitation 
of  the  average  five-year,  one-hour  storm  for  the 
general  region  was  obtained  from  Technical  Paper 
No.  24,  Part  II  (U.  S.  Weather  Bureau,  1954).  Data 
are  presented  by  state,  section,  and  station  in 
Table  5. 

The  Thornthwaite  Precipitation- Effectiveness 
Index,  or  P-E  index,  is  a  measure  of  the  availabil¬ 
ity  of  moisture  to  vegetation,  and  depends  on  the 
amount  and  distribution  of  precipitation  and  evapor¬ 
ation;  the  latter  in  turn  depends  on  variations  of 


temperature  and  wind  movement  (Thornthwaite, 
1931).  The  exact  rational  formula  is 

i  =  ioi:_p 

12  E  ’ 

where  P  is  the  average  precipitation  for  each 
month,  and  E  is  the  average  evaporation  for  each 
month.  The  index,  I,  is  therefore  a  dimensionless 
ratio.  An  empirical  formula  that  is  an  approxima¬ 
tion  to  I,  when  no  evaporation  data  are  available, 
is 


I  =  115  L 

12 
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where  T  is  the  monthly  average  daily  temperature. 

Thornthwaite  divided  the  range  of  values  of  I  as 
follows: 


Characteristic  vegetation  P-E  Index 

Rain  forest  >  128 

Forest  64  -  127 

Grassland  32  -  63 

Steppe  16  -  31 

Desert  ^  15 


The  P-E  index  for  each  station  for  which  sufficient 
data  are  available,  or  for  which  temperatures  can 
be  estimated  by  interpolation  from  nearby  stations, 
is  provided  in  Table  5. 

In  the  areas  of  the  West  for  which  data  were 
obtained,  the  precipitation  increases  with  eleva¬ 
tion,  above  a  lower  limit,  and  temperature  decreases 
(Turnage  and  Mallery,  1941).  Within  a  limited  area, 
the  P-E  indices,  when  plotted  against  elevation, 
show  a  significant  linear  increase.  The  points  on 
the  graphs  deviate  from  linearity  in  Colorado  and 
Utah  for  elevations  greater  than  10,000  feet  and 
8,000  feet  respectively.  The  graphs  obtained  were 
then  used  to  estimate  the  mean  P-E  index  for  each 
basin.  Ignored  are  the  differences  of  precipitation 
effectiveness  on  north  and  south-facing  hillsides, 
the  effect  of  nearby  mountains,  and  the  situation  of 
the  station.  The  value  of  I  is,  therefore,  only  ap¬ 
proximate.  However,  the  results  agree  well  with 
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known  vegetation  in  the  basins.  Figure  7  shows 
graphs  of  P-E  index  against  elevation.* 


P-E  Index 


Figure  7  (a) 

P~E  Index 


Figure  7  (b) 

*Lines  or  curves  fitted  by  eye  only. 


P-E  Index 


Figure  7  (c) 

Figure  7.  Variation  of  precipitation— effectiveness  index 
with  elevation; 

(a)  in  Park  Park  Plateau  and  Front  Range,  Colo,  and 
Wasatch  Mtns.,  Utah, 

(b)  in  southeastern  Arizona  and  southwestern  New  Mex¬ 
ico, 

(c)  in  north-central  and  west-central  New  Mexico. 

Intensity- duration-frequency  relations  are  ex¬ 
pressed  by  the  amount  of  rain  received  from  a 
storm  of  one  hour  duration  expected  to  occur  once 
in  five  years.  The  information  was  obtained  from 
Technical  Paper  24,  Part  II  (U.  S.  Weather  Bureau, 
1954,  p.  6-8),  by  multiplying  the  two-year  one-hour 
amounts  as  given  on  isopleth  maps,  by  1.4  which 
is  the  ratio  of  the  amounts  for  the  two  return  per¬ 
iods  in  the  western  United  States. 

Season  of  maximum  rainfall  is  secured  directly 
from  the  tabulated  amounts  of  precipitation  for 
each  month. 

Ratio  of  mean  January  precipitation  to  1/12 
mean  annual  precipitation  is  obtained  directly  from 
the  values  for  mean  annual  precipitation,  and  mean 
January  precipitation.  The  measure  is  important 
because  of  the  greater  effectiveness  of  rain  for 
erosion  when  the  vegetal  cover  is  least  extensive 
as  it  is  in  this  winter  month.  The  calculated  value 
was  found  to  be  within  the  range  of  the  value  given 
by  the  map  in  Climatic  Atlas  of  the  United  States 
(Visher,  1954,  p.  204). 


General  Description  of  Basins 


This  section  treats  general  aspects  of  each 
basin.  Obtained  largely  from  map  and  field  inspec¬ 
tion,  the  description  is  provided  for  reference 
only;  no  effort  is  made  to  discuss  the  contribution 
of  the  information  on  each  basin  in  regard  to  the 
conclusions.  Location,  physiographic  province  and 
section,  type  of  surroundings,  grazing  history  and 
range  evaluation,  lithology  and  structure,  presence 
of  gullying,  type  and  extent  of  cover,  and  other 
features  are  best  explained  verbally,  since  most 
of  the  information  is  non-numerical.  Climatologic 
data  are  included  in  the  form  of  estimated  precipi¬ 
tation-effectiveness  index,  and  the  distribution  of 
rainfall  by  months.  The  basins  are  listed  alpha¬ 
betically  by  state  and  quadrangle  in  which  they  are 
located.  Maps  of  all  basins  are  given  in  Appendix 
IV. 


BASINS  IN  ARIZONA 

Chinle  Badlands,  Cameron,  including  two  major 
basins  in  this  locality,  are  covered  by  a  special  map 
made  in  August,  1953,  by  E.  D.  Koons,  S.  A. 

Schumm,  and  A.  N.  Strahler,  at  the  scale  one  inch 
to  50  feet  (1:600).  The  topography  is  “micro-relief,” 
that  is,  an  area  of  fully-developed,  mature  topo¬ 
graphy,  relatively  rugged,  but  composed  of  ex¬ 
tremely  small  sized  units.  The  area  is  located  at 
the  edge  of  the  Navajo  section  of  the  Colorado 
Plateau,  in  the  Painted  Desert.  The  lithology  is  the 
blue-gray  Chinle  shale,  in  thickly  platy  laminae  of 
silt  and  clay,  and  fine  sand,  containing  nodules  of 
gypsum  and  thin  plates  of  iron- stained  chalcedony. 
Beds  are  essentially  flat-lying.  No  vegetation  what¬ 
soever  is  present  anywhere  in  these  basins.  The 
surface  layer  of  soft,  unstable  sand  and  clay,  one 
to  four  inches  thick,  with  numerous  drying  cracks, 
has  been  produced  from  the  original  shale  by  swell¬ 
ing  of  the  clay  during  rains,  and  partial  desiccation. 
Any  disturbance  of  the  surface  produces  a  cascade 
of  the  crust  material,  exposing  a  spot  of  bare,  un¬ 
altered  shale,  and  burying  the  foot  of  the  slope 
more  deeply.  (See  also,  Schumm,  1954,  p.  80). 

The  nearest  weather  station  at  Tuba  City,  about 
18  miles  north,  is  probably  in  a  meteorologically 
similar  area.  Of  the  yearly  rainfall,  46%  occurs  in 
the  months  July  through  October.  The  region  is 
definitely  desertic,  with  a  P-E  index  of  12.8. 

Finley  and  Adams  Canyon,  and  Sycamore  Canyon 
(Harshaw  Quadrangle):  Four  fourth-order  basins 


are  included  in  this  section.  Field  work  was  not  de¬ 
tailed;  only  infiltration  measurements  were  taken 
in  these  areas.  The  basins  lie  in  the  Patagonia 
Mountains,  in  the  Basin  and  Range  province,  Mexi¬ 
can  Highland  section.  Underlying  the  basins  is  a 
granite  complex  that  is  deeply  weathered  in  places, 
with  many  residual  boulders  and  dikes  standing  in 
relief.  No  gullying  is  present.  Vegetation  varies 
from  oak  chaparral  on  south-facing  slopes,  to  oak 
woodland  and  pinon-juniper  on  the  north-facing  and 
higher  slopes.  Grass  is  plentiful  in  open  areas  and 
is  ungrazed.  About  85%  of  the  area  is  in  forest  and 
bushes,  12%  to  14%  is  grassland,  and  the  remaining 
is  bare  soil  or  rock.  The  P-E  index  is  estimated  at 
40  to  65,  depending  on  elevation  and  exposure.  The 
nearest  source  of  weather  records  is  Nogales,  14 
miles  southwest,  where  68%  of  the  yearly  precipi¬ 
tation  is  in  July  through  October. 

Tributary  Basins  to  Harshaw  Creek  on  Red 
Mountain  (Mt.  Hughes  Quadrangle),  located  in  Secs. 
16,  17,  20,  21,  Twp.  22  S.,  Ra.  16  E.,  are  one  fourth- 
order  basin  and  two  third-order  basins.  Red  Moun¬ 
tain  is  very  near  the  boundary  between  the  Sonoran 
Desert  and  Mexican  Highland  sections  of  the  Basin 
and  Range  province,  and  stands  prominently  2,300 
feet  above  the  Sonoita  Creek  valley.  The  rock  is 
wholly  igneous;  a  thick  series  of  mineralized 
felsite  overlies  prophyritic  basalt.  Cliffs  form 
easily  in  both  series,  and  slide  rock  slopes  are 
common.  Several  mine  drifts  are  present,  but  none 
are  active. 

The  main  cover  in  the  basins  is  oak  chaparral 
and  brush,  grading  downward  into  semi-desert 
shrub  and  grassland,  and  upward  into  pinon-juniper 
woodland,  and  near  the  top  of  Red  Mountain  on  north¬ 
facing  slopes,  larger  Arizona  white  oak.  An  esti¬ 
mated  50%  of  the  area  is  in  brush  and  chaparral, 

30%  is  grassland,  and  20%  bare  rock.  The  grassland 
area  is  grazed  lightly,  and  the  range  appeared  to  be 
in  fair  to  good  condition. 

The  P-E  index  is  estimated  to  be  35  to  45,  de¬ 
pending  on  elevation  and  orientation.  The  nearest 
weather  station  is  at  Patagonia,  three  miles  distant. 
There,  63%  of  the  yearly  precipitation  occurs  in  the 
months  July  through  October. 

Whipple  Basin  (Prescott  Quadrangle)  is  a  fourth- 
order  basin  in  the  northwest  quarter  of  Sec.  26, 

Twp.  14  N.,  Ra.  2  W.,  at  the  north  boundary  of  the 
Whipple  Military  Reservation,  in  the  Mexican  High¬ 
land  section,  Basin  and  Range  province.  The  basin 
lies  in  an  extensive  structural  valley,  three  miles 
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from  the  north  end  of  the  Bradshaw  Mountains.  The 
immediate  relief  in  the  area  is  low,  not  more  than 
250  feet.  The  rock  is  biotite  granite,  not  greatly 
weathered  where  exposed,  and  overlain  by  thin 
well- indurated  bouldery  alluvium.  The  cover  is  es¬ 
timated  to  be  67%  grassland,  30%  trees  and  bushes, 
and  3%  bare  rock.  Grazing  is  moderate,  and  the 
range  is  in  fair  condition.  Few  new  gullies  have 
formed  on  the  hillsides,  and  the  channels  show 
rejuvenation,  but  are  not  trenched. 

Prescott,  about  two  miles  to  the  southwest, 
probably  has  very  nearly  the  same  climate.  The 
P-E  index  is  40,  and  the  months  of  greatest  pre¬ 
cipitation  are  July,  August,  September,  December, 
January,  and  February,  with  74%  of  the  annual 
precipitation  occurring  in  those  six  months. 

Hog  Canyon  Basin  (Sonoita  Quadrangle)  is  a 
fourth-order  basin  tributary  to  Hog  Canyon,  and  ex¬ 
tends  across  Sec.  20,  Twp.  20  S.,  Ra.  16  E.  The 
area  is  in  a  dissected  alluvial  fan  on  the  east  flank 
of  the  Santa  Rita  mountains,  in  the  Mexican  High¬ 
land  section,  Basin  and  Range  province.  The  im¬ 
mediate  relief  is  500  feet,  but  the  mountains  rise 
sharply  to  the  west  to  a  height  of  4,600  feet  above 
the  basin,  doubtless  causing  increased  precipita¬ 
tion  from  maritime  tropical  air  masses  from  the 
south  and  southeast.  South-facing  slopes  are  al¬ 
most  entirely  grassland,  while  the  north-facing 
slopes  have  only  small  grassland  areas.  About  55% 
of  the  basin  is  grassland,  45%  is  oak  woodland, 
and  less  than  1%  is  bare.  Grazing  is  light;  the 
range  is  fair  to  good,  with  much  long  bunch  grass 
remaining.  Gullies  are  present  in  some  small 
channels  and  in  many  places  extend  into  fresh 
slopes. 

The  nearest  weather  records  are  from  the 
Crown  C  ranch,  one  mile  east,  where  precipitation 
data  are  available.  The  P-E  index  is  estimated  to 
be  35  to  40,  depending  on  exposure.  Greatest  pre¬ 
cipitation  is  in  the  months  July  through  October, 
when  62%  of  the  annual  rainfall  occurs. 

Basins  in  Saguaro  National  Monument  (Tucson 
Quadrangle,  1904):  Two  fourth-order  basins  near 

the  west  boundary  of  the  Monument  and  one-half 
mile  north  of  the  headquarters,  were  mapped  by 
the  writer  and  assistant  in  December,  1955,  on  a 
scale  100  feet  to  the  inch  (1:1,200).  The  area  is  in 
the  Sonoran  Desert  section,  Basin  and  Range  prov¬ 
ince.  The  basins  are  in  an  area  of  low  relief,  a 
dissected  pediment,  near  the  foot  of  the  Tanque 
Verde  Mountains.  Rock  type  is  mainly  chlorite  and 
sericite  schist,  with  injected  fine-grained  granite. 
Attitude  of  the  planes  of  schistosity  is  variable, 
but  the  strike  averages  N  10°  E,  dip  25°  to  45°  W. 
About  60%  to  70%  of  the  area  is  bare  soil  and  rock, 
the  remainder  is  covered  by  saguaro  “forest.”  The 
prickly  pear  is  over-grazed,  with  three  animals 
(horses)  per  square  mile,  and  the  range  is  very 
poor.  As  the  result  of  over-grazing  and  disturbance 


of  the  surface  by  the  animals’  hooves,  numerous 
gullies  have  become  established  on  the  hillsides, 
and  a  few  small  channels  deepened.  Main  channels 
are  broad  and  filled  with  sand  and  gravel. 

The  climate  data  for  this  area  were  taken  from 
Tucson  Magnetic  Observatory  and  the  old  Vail  and 
Pantano  records.  The  P-E  index  is  estimated  at 
15.  Of  the  annual  rainfall,  52%  is  in  the  months 
July  through  October. 


BASINS  IN  COLORADO 

Dory  Hill  Basin  and  Mesa  Gulch  (Black  Hawk 
Quadrangle)  are  two  fourth-order  basins  in  Secs. 

32,  33,  Twp.  2  S.,  Ra.  72  W.,  and  Secs.  4,  5,  8,  9, 
Twp.  3  S.,  Ra.  72  W.,  in  the  Colorado  Front  Range, 
Southern  Rocky  Mountains  province.  The  basins 
lie  in  a  dissected  pediment  or  other  old  erosion 
surface  of  moderate  relief,  sloping  eastward.  The 
immediate  relief  is  1,500  feet.  Higher  mountains 
near  the  Continental  Divide,  10  to  15  miles  to  the 
west,  reach  heights  of  4,500  feet  above  the  basin 
level. 

The  rock  underlying  the  basins  is  granite-gneiss 
and  schist,  with  a  probable  general  northward  dip. 
About  60%  of  the  area  is  grassland  and  is  grazed, 
35%  to  37%  is  in  forest,  the  remainder  is  bare.  The 
estimated  P-E  index  is  76  to  80,  according  to  the 
nearest  records  taken  from  Frances  weather  sta¬ 
tion  at  nearly  the  same  elevation.  The  months 
March  through  July  receive  54%  of  the  mean  an¬ 
nual  precipitation,  though  rainfall  is  distributed 
rather  evenly  throughout  the  year. 

Cabrini  Gulch  (Morrison  Quadrangle)  is  a  small 
third- order  basin  below  the  Mother  Cabrini  Shrine, 
off  Highway  40,  in  the  southeast  corner  of  Sec.  16, 
Twp.  4  S.,  Ra.  70  W.,  and  is  tributary  to  Mt.  Ver¬ 
non  Canyon.  It  is  about  a  mile  west  of  the  moun¬ 
tain  front,  in  the  Colorado  Front  Range,  Southern 
Rocky  Mountains  province.  The  rock  underlying 
the  basin  is  hard,  black,  biotite-hornblende  schist, 
with  infrequent veinlets  of  feldspathic  material. 
About  95%  of  the  basin  is  grassland,  the  remainder 
is  covered  by  small  trees  and  bushes.  Grazing  is 
moderate  to  heavy,  and  the  range  is  fair.  The  main 
channel  of  the  basin  is  gullied  to  a  depth  of  20  feet 
near  the  mouth  and  less  than  one  foot  near  the  head. 
A  single  gully  in  the  slopes  was  caused  by  concen¬ 
tration  of  runoff  from  the  road  which  passes 
across  the  top  of  the  basin.  The  P-E  index  is  es¬ 
timated  at  46  on  the  basis  of  elevation  alone,  al¬ 
though  this  figure  may  be  high,  considering  that  a 
large  percentage  of  the  basin  faces  south.  The 
months  of  greatest  rainfall  are  April  through 
August,  when  61%  of  the  yearly  precipitation  occurs. 

Green  Mountain  Basin  (Morrison  Quadrangle)  is 
a  small  third-order  basin  on  the  southeast  end  of 
Green  Mountain,  in  the  center  of  Sec.  19,  Twp.  4S., 
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Ra.  69  W.,  near  the  western  boundary  of  the  Colo¬ 
rado  Piedmont  section,  Great  Plains  province, 
about  two  miles  east  of  the  mountain  front.  Under¬ 
lying  the  basin  is  a  coarse,  bouldery  conglomerate 
(the  Green  Mountain  conglomerate)  of  alluvial  ori¬ 
gin,  consisting  of  igneous  and  metamorphic  rock 
fragments.  Below  the  conglomerate  is  a  partially 
consolidated,  greenish-brown  crossbedded  sand. 

The  mountain  is  a  remnant  of  Tertiary  alluvium, 
protected  by  the  resistant  conglomerate.  An  ex¬ 
posure  of  granite  cobbles  reveals  partial  disinte¬ 
gration  to  a  depth  of  two  and  one  half  feet.  Almost 
all  cobbles  have  a  caliche  coating  on  the  under 
sides. 

About  95%  of  the  area  is  grassland,  3%  to  4%  has 
small  trees  and  brush,  and  the  remainder  is  bare. 
The  basin  is  moderately  grazed,  and  the  range  is 
fair  to  good.  The  main  channel,  about  8  feet  deep, 
is  trenched  and  has  a  trapezoidal  cross-section. 

The  smaller  channels,  higher  in  the  basin,  are 
gullied  less  deeply  and  have  V-shaped  cross-sec¬ 
tions. 

The  P-E  index  is  estimated  at  38  to  44.  The  pre¬ 
cipitation  is  assumed  to  be  the  same  as  that  of 
Morrison,  3  miles  southwest,  where  61%  of  the  an¬ 
nual  precipitation  occurs  in  the  months  of  April 
through  August. 

Delonde  Gulch  (Nederland  Quadrangle)  is  a 
fourth-order  basin  in  Secs.  3  and  4,  Twp.  1  S., 

Ra.  73  W.,  in  the  Colorado  Front  Range.  The  im¬ 
mediate  area  has  about  1,800  feet  of  relief,  while 
six  or  seven  miles  west  the  peaks  of  the  Conti¬ 
nental  Divide  rise  to  heights  of  4,000  feet  above 
the  basin.  The  rock  in  the  area  is  granite-gneiss 
and  schist,  disordered,  but  with  a  general  north¬ 
ward  dip.  Considerable  mining  was  done  immed¬ 
iately  west  of  the  basin  until  the  late  1930’s, 
though  no  mines  were  seen  in  the  basin  itself.  More 
recently  logging  has  become  the  major  activity, 
but  its  eventual  extent  cannot  be  predicted.  Forest 
covers  nearly  100%  of  the  basin.  The  estimated 
P-E  index  is  84,  and  the  precipitation  is  heaviest 
in  the  months  of  March  through  August,  according 
to  the  nearest  weather  records,  at  Gold  Hill. 

Hicks  Gulch  (Nederland  Quadrangle)  is  two 
miles  south  of  Delonde  Gulch,  in  Secs.  15,  16,  17, 
Twp.  1  S.,  Ra.  73  W.,  at  about  the  same  elevation, 
and  is  very  similar  to  it,  except  as  described  be¬ 
low.  Many  inactive  mine  drifts  are  present  in 
Hicks  Gulch.  To  provide  timber  for  mines,  most  of 
the  large  limber  pines  in  the  upper  part  of  the 
basin  were  logged  35  to  50  years  ago.  The  basin  is 
grazed  lightly,  and  the  range  is  good  to  excellent 
in  the  grasslands.  Aerial  photographs  show  about 
75%  of  the  basin  covered  by  forest  and  bushes,  20% 
to  23%  by  grass,  and  the  remainder  bare.  There 
are  no  gullies,  but  rejuvenation  of  the  main  channel 
is  evident  by  the  steep-walled  inner  channel  which 


becomes  considerably  deeper  toward  the  mouth.  In 
each  of  the  two  main  tributaries,  small  levee-like 
features  extend  along  the  channel  for  nearly  a  quar¬ 
ter  of  a  mile,  and  are  probably  the  result  of  a  late 
Pleistocene  earth-flow  or  solifluction  episode.  The 
basin  was  probably  never  glaciated,  since  no  spur 
trimming  or  other  alteration  of  the  valley,  and  no 
cirque  or  nivation  cirque  is  present.  The  broad, 
poorly-drained  upland  in  the  highest  part  of  the 
basin  was  affected  greatly  by  periglacial  phenomena, 
such  as  frost  heaving,  and  even  now  soil  moisture 
may  be  frozen  a  few  feet  below  the  surface  during 
much  of  the  year.  The  P-E  index  is  estimated  to  be 
79  to  85,  depending  on  elevation  and  exposure. 

Moon  Gulch  (Nederland  Quadrangle)  is  in  Secs. 

2,  3,  10,  11,  Twp.  2  S.,  Ra.  73  W.,  about  four  miles 
SSE  of  Hicks  Gulch,  at  about  the  same  elevation. 
Underlying  the  basin  is  the  same  granite-gneiss 
and  schist  as  in  the  two  basins  farther  north  in  this 
quadrangle,  and  the  planes  of  schistosity  have  a 
general  northward  dip.  The  greater  height  of  the 
southern  part  of  the  basin  is  due  to  the  presence  of 
a  fine-grained  granite  stock,  resistant  to  erosion. 
Around  this  stock  are  numerous  inactive  mine  ex¬ 
cavations.  Considerable  portions  of  the  basin  have 
been  burned  over,  some  rather  recently,  and  are 
gullied  extensively.  About  50%  of  the  basin  area  is 
in  forest  cover,  35%  burned,  and  the  remainder  is 
grassland.  The  grasslands  are  primarily  moderately 
grazed  meadows  near  the  main  channel,  and  the 
range  is  in  good  to  excellent  condition;  elsewhere 
the  range  is  poor,  and  not  grazed.  The  main  channel 
has  produced  a  broad  floodplain,  bounded  by  ter¬ 
races  of  20  to  30  feet  in  height.  The  floodplain  nar¬ 
rows  rapidly  upstream,  and  is  not  apparent  more 
than  a  mile  above  the  mouth. 

The  P-E  index  is  estimated  to  be  78  to  85,  de¬ 
pending  on  elevation  and  exposure.  According  to  the 
nearest  records  of  the  Silver  Lake  weather  station, 
the  months  with  greatest  precipitation  are  March 
through  July,  when  54%  of  the  annual  amount  occurs. 

Tributaries  to  Saruche  Canyon  and  Gallinas 
Creek  (Starkville  Quadrangle)  comprise  two  basins, 
one  third-order  and  one  fourth-order,  in  the  Park 
Plateau,  Raton  section,  Great  Plains  province.  The 
tributary  to  Gallinas  Creek  lies  in  Secs.  27,  34, 

Twp.  34  S.,  Ra.  64  W.;  the  tributary  to  Saruche 
Canyon  lies  partially  in  the  southwest  quarter  of 
Sec.  32,  Twp.  34  S.,  Ra.  64  W.  The  rock  is  alternat¬ 
ing  thick- bedded  sandstone  and  thinner  shale,  and 
coal  strata.  The  beds  are  flat-lying,  and  cliffs  are 
common  on  the  sandstone.  Considerable  coal  min¬ 
ing  operations  are  carried  on  in  the  area,  though 
not  in  these  basins.  The  estimated  P-E  index  is  48. 
The  nearest  source  of  weather  records  is  Madrid, 

10  miles  northwest,  where  66%  of  the  annual  pre¬ 
cipitation  occurs  during  the  period  April  through 
August. 
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BASINS  IN  NEW  MEXICO 

Cottonwood  Canyon  (Allie  Canyon  Quadrangle) 
lies  mainly  in  Secs.  1,  2,  3,  Twp.  16  S.,  Ra.  12  W., 
in  the  Mexican  Highland  section,  Basin  and  Range 
province.  This  area  was  not  included  in  the  field 
study.  The  P-E  index  is  38,  estimated  from  the 
elevation.  Weather  data  are  available  from  the 
Mimbres  Ranger  Station,  one  mile  away.  The  great¬ 
est  rainfall  occurs  in  July  through  September, 
when  52%  of  the  annual  precipitation  occurs. 

Tributaries  to  Sapillo  Creek  (Allie  Canyon 
Quadrangle)  are  two,  small,  fourth-order  basins 
in  Secs.  20,  21,  29,  Twp.  15  S.,  Ra.  12  W.,  in  the 
Mexican  Highland  section,  Basin  and  Range  prov¬ 
ince.  These  areas  were  not  visited  for  field  study. 
The  estimated  P-E  index  is  40,  and  the  rainfall 
characteristics  are  assumed  to  be  the  same  as 
above  for  the  Mimbres  Ranger  Station,  six  miles  to 
the  southeast. 

Sawmill  Canyon  and  Tributaries  (Big  Burro 
Mountains  Quadrangle)  are  three  fourth-order  ba¬ 
sins  in  Secs.  11,  12,  13,  Twp.  20  S.,  Ra.  16  W.,  and 
Secs.  6,  7,  18,  Twp.  15  W.,  Ra.  20  S.,  in  the  Mexi¬ 
can  Highland  section,  Basin  and  Range  province. 
They  extend  over  a  large  portion  of  the  Big  Burro 
Mountains,  where  the  area  has  a  relief  of  2,300 
feet.  The  Big  Burro  Mountains  are  an  isolated 
granite  mass,  which  is  structurally  complex,  con¬ 
taining  many  dikes  of  fine  to  coarse  texture  rang¬ 
ing  in  composition  from  granitic  to  ultramafic. 
Wherever  exposed,  the  granite  is  deeply  and  thor¬ 
oughly  weathered.  Forest  of  pines  and  oak  brush 
covers  about  40%  of  the  area,  the  remainder  has 
sparse  bunch  grass  or  is  completely  bare,  and 
the  range  is  poor.  Much  gullying  has  occurred  on 
the  slopes,  and  some  in  the  small  channels.  The 
estimated  P-E  index  is  40  to  55,  varying  with  ex¬ 
posure  and  elevation.  From  the  daily  weather  rec¬ 
ords  kept  for  seven  years  at  the  Smithsonian  Insti¬ 
tution  Astrophysical  Observatory  sunspot  labora¬ 
tory  (Tyrone  Station)  on  top  of  Jacks  Mountain, 
near  the  head  of  these  basins,  the  P-E  index  was 
54.  Greatest  precipitation  is  from  July  through 
September,  and  a  lesser  maximum  in  December 
through  February;  74%  of  the  yearly  precipitation 
occurs  during  those  months. 

Walnut  Creek  (Big  Burro  Mountains  Quadrangle) 

is  in  the  northeast  quarter  of  Sec.  7,  Twp.  20  S., 

Ra.  15  W.,  on  the  northeast  slope  of  the  Big  Burro 
Mountains.  These  two  third-order  basins  differ 
from  Sawmill  Canyon  in  their  greater  relative  re¬ 
lief  and  abundance  of  vegetation;  forest  covers 
about  85%  of  the  area,  the  remainder  is  bare.  The 
estimated  P-E  index  is  50  to  55.  In  other  respects 
these  basins  are  essentially  like  the  others  in  the 
Big  Burro  Mountains. 

Paso  del  Norte  Basin  (Bland  Quadrangle)  occu¬ 
pies  half  of  Sec.  2,  Twp.  18  N.,  Ra.  4  E.,  and  is 


near  the  boundary  between  the  Southern  Rocky 
Mountains  province  and  the  Mexican  Highland  sec¬ 
tion  of  the  Basin  and  Range  province.  It  is  on  the 
edge  of  the  Jemez  Caldera,  and  lies  in  an  area  of 
volcanic  dacite,  basalt,  and  pumice  beds,  dipping 
southeastward.  Immediately  to  the  north  and  west, 
the  Sierra  de  los  Valles  drops  700  feet  into  the 
Valle  Grande,  which  is  the  caldera  proper.  Else¬ 
where  along  the  Sierra,  peaks  reach  elevations  of 
10,200  feet.  Almost  100%  of  the  basin  is  forest 
covered.  The  estimated  P-E  index  is  98.  The  near¬ 
est  weather  station  is  Los  Alamos.  July  and  August 
are  the  months  of  maximum  rainfall,  when  35%  of 
the  yearly  precipitation  occurs. 

Peters  Dome  Basin  (Bland  Quadrangle)  is  a 
fourth-order  tributary  of  Capulin  Canyon  immedi¬ 
ately  east  of  Sec.  20,  Twp.  18  N.,  Ra.  5  E.,  in  the 
San  Miguel  Mountains,  in  the  Mexican  Highland 
section,  Basin  and  Range  province.  The  rocks  are 
nearly  horizontal,  alternating  basalt  and  pumice. 
Though  on  all  sides  of  the  basin  the  basalt  forms 
cliffs  800  to  900  feet  in  height,  the  basin  itself  is 
free  of  larger  cliffs;  only  near  its  mouth  has  it  be¬ 
gun  to  erode  into  the  main  cliff-forming  layer. 
Forest  and  bushes,  partially  burned  50  to  75  years 
ago,  cover  about  94%  of  the  basin;  4%  is  grassland, 
and  the  remainder  is  bare.  The  basin  is  poor  range 
and  is  not  grazed;  however,  it  does  serve  as  a 
main  stock  drive- way  for  herds  from  the  west.  No 
gullies  are  present.  The  estimated  P-E  index  is 
65  to  75,  depending  on  elevation  and  exposure.  The 
nearest  source  of  weather  records  of  the  same 
elevation  is  10  miles  northeast  at  Los  Alamos;  dis¬ 
tribution  of  precipitation  is  probably  the  same  as 
for  Paso  del  Norte  basin. 

Cerro  Pavo  Basin  (Cerro  del  Grant  Quadrangle) 
is  a  large  third-order  basin  in  Secs.  7,  18,  19,  Twp. 
21  N.,  Ra.  4  E.,  in  the  southernmost  end  of  the 
Southern  Rocky  Mountains  province.  The  basin  lies 
atop  a  high  lava- capped  mesa,  which  falls  away 
with  considerable  relief  to  the  west,  north,  and 
east,  and  rises  southward  to  the  smooth,  forested 
peaks  of  the  north  rim  of  the  Jemez  Caldera.  Un¬ 
derlying  the  basin  is  gray  to  pink  dacite  or  rhyolite 
porphyry,  with  occasional  interbedded  scoria,  all 
essentially  flat-lying.  About  80%  of  the  basin  is 
forest  covered;  the  remainder  is  grassland.  The 
range,  grazed  by  sheep  and  cattle,  is  fair  to  good 
in  the  meadows  near  the  main  channels.  No  gullies 
are  present,  and  small  channels  are  poorly  defined. 
The  P-E  index  is  110,  estimated  from  the  elevation. 
The  distribution  of  rainfall  is  assumed  to  be  simi¬ 
lar  to  that  of  Capulin  Ranger  Station,  where  the 
months  of  greatest  rainfall  are  May  through  Septem¬ 
ber,  when  74%  of  the  annual  precipitation  occurs. 

Basins  in  Hurley  East  Quadrangle  are  located 
in  Secs.  15,  22,  23,  Twp.  18  S.,  Ra.  12  W.,  north¬ 
east  of  the  mining  town  of  Hurley.  These  three 
third-order  and  one  fourth-order  basins  were  not 
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visited  in  the  field.  A  published  geologic  map 
(Paige,  1916)  shows  that  Tertiary  basalt  and  ande¬ 
site,  overlying  gravel  and  tuff,  form  the  higher 
parts  of  the  basins.  The  P-E  index  is  31  to  35,  es¬ 
timated  from  elevation  only.  The  nearest  source  of 
rainfall  data,  Ft.  Bayard,  records  that  54%  of  the 
yearly  precipitation  is  in  the  months  July  through 
September. 

Hairpin  Basin  (Jarosa  Quadrangle)  is  a  small 
third-order  basin  in  Secs.  1,  2,  Twp.  21  N.,  Ra.  2 
E.,  and  is  in  the  southernmost  part  of  the  Southern 
Rocky  Mountains  province.  The  main  rock  type  is 
thickly- bedded  red  sandstone,  essentially  hori¬ 
zontal,  cropping  out  in  low  cliffs  on  south-facing 
slopes.  About  90%  of  the  area  is  covered  by  brush 
and  forest,  5%  is  grassland,  and  5%  is  bare  soil 
and  rock.  The  basin  is  moderately  grazed,  and  the 
range  on  north-facing  slopes  is  good  to  fair.  The 
only  gullies  present  are  those  formed  as  a  result 
of  recent  road  building.  The  P-E  index  is  72,  es¬ 
timated  from  the  elevation.  Capulin  Ranger  Station 
is  the  nearest  source  of  weather  data;  there,  74% 
of  the  annual  precipitation  occurs  in  May  through 
September. 

Poleo  Basin  (Jarosa  Quadrangle)  is  a  fourth- 
order  basin  in  Secs.  27,  28,  Twp.  22  N.,  Ra.  2  E., 
in  the  southernmost  part  of  the  Rocky  Mountains 
province,  near  the  northern  limit  of  the  plateau 
comprising  the  San  Pedro  Mountains.  The  topograpy 
drops  rapidly  northward  into  the  Rio  Puerco  Valley. 
The  rock  types  are  reddish,  coarse-to-conglomer- 
atic  arkosic  sandstone,  gray  quartzose  sandstone, 
crinoidal  limestone,  and  shale.  In  the  greater  part 
of  the  area,  the  beds  are  flat- lying,  but  dip  east¬ 
ward  10°  to  20°  in  the  higher  western  part  of  the 
basin.  Trees  and  brush  cover  an  estimated  96%  of 
the  area,  grass  1%,  and  the  remainder  is  bare  soil 
and  rock.  Grazing  is  light  to  moderate;  the  range  is 
poor.  Gullies  are  present  in  the  main  channel  and 
one  of  the  main  tributaries.  They  are  V-shaped, 
one  to  eight  feet  deep,  cut  in  alluvium,  but  not  bed¬ 
rock.  The  estimated  P-E  index  is  77,  and  rainfall 
distribution  is  assumed  to  be  the  same  as  for  Capu¬ 
lin  Ranger  Station. 

Quartz  Hill  Basin  (Paxton  Springs  Quadrangle) 
is  a  fifth-order  basin  in  Secs.  19,  20,  29,  30,  Twp. 

10  N.,  Ra.  11  W.,  on  the  eastern  flank  of  the  Zuni 
Mountains,  in  the  Datil  section  of  the  Colorado 
Plateau  province.  The  rock  underlying  the  entire 
basin  is  a  gneissic  granite,  with  occasional  sept¬ 
like  layers  of  muscovite- chlorite  schist,  dipping 
steeply  southwestward.  Of  several  mines  present, 
only  one  is  still  active.  The  granite  is  deeply 
weathered  in  most  places,  especially  at  the  eastern 
end  of  the  basin  where  remnants  of  arkose  origin¬ 
ally  deposited  directly  on  the  granite  are  found. 

The  greater  decomposition  here  can  be  attributed 
to  pre-arkose  weathering.  Trees  and  bushes  cover 
an  estimated  80%  of  the  basin,  grass  15%,  and  5% 


is  bare  soil  and  rock.  The  area  is  grazed  moder¬ 
ately  to  heavily;  range  is  poor  to  fair.  The  few 
gullies  present  are  in  the  smallest  channels.  The 
P-E  index  is  44,  estimated  from  the  elevation.  The 
distribution  of  rainfall  is  assumed  to  be  similar  to 
that  of  El  Morro  National  Monument,  18  miles  west, 
and  San  Rafael,  15  miles  northeast.  The  months  of 
greatest  precipitation  are  July  through  September, 
when  48%  of  the  annual  rainfall  occurs. 

Camino  Basin  (Santa  Fe  Quadrangle)  is  a  small 
fourth-order  basin  located  in  the  northeast  quarter 
of  Sec.  28,  Twp.  17  N.,  Ra.  9  E.,  in  the  Mexican 
Highland  section  of  the  Basin  and  Range  province, 
in  the  valley  of  the  Santa  Fe  River.  The  lithology  is 
the  Santa  Fe  formation,  indurated  sand  and  gravel  of 
alluvial  origin.  About  15%  of  the  basin  is  covered  by 
trees  and  bushes,  30%  by  bunch  grass,  and  55%  is 
bare.  Grazing  is  very  heavy,  mainly  by  sheep  and 
goats,  resulting  in  very  poor  range  conditions.  The 
main  channels  are  gullied,  with  vertical  sides.  The 
P-E  index  is  between  that  for  Santa  Fe  Airport  and 
Santa  Fe  City,  or  about  23.  The  distribution  of 
rainfall  is  assumed  to  be  about  the  same  as  that  of 
Santa  Fe  City,  where  44%  of  the  annual  precipitation 
is  in  July  through  September. 

Tano  Basin  (Santa  Fe  Quadrangle)  is  a  fourth- 
order  basin  essentially  like  Camino  Basin  in  lithol¬ 
ogy.  It  lies  north  of  the  city  of  Santa  Fe  proper,  but 
within  the  city  limits,  and  is  a  tributary  to  Arroyo 
de  los  Frijoles.  Trees  and  bushes  cover  80%  of  the 
basin  area;  grass  4%,  and  16%  is  bare.  The  range  is 
very  poor  to  fair,  and  has  probably  been  grazed  in 
the  recent  past,  though  not  immediately  prior  to  the 
time  of  field  study.  Most  of  the  channels  are  trenched, 
and  gullies  have  become  established  on  the  slopes. 

The  main  channel  is  deepened  to  15  feet  near  its 
head,  but  is  not  gullied  toward  the  mouth;  there,  a 
broad,  sandy  channel  is  present.  All  gully  sides  are 
nearly  vertical;  the  bottoms  are  usually  flat.  The 
P-E  index  is  estimated  at  33,  somewhat  above  that 
of  Santa  Fe.  The  distribution  of  precipitation  is  the 
same  as  for  Camino  Basin. 

Cedro  Peak  Basin  (Sedillo  Quadrangle)  is  a 
fourth-order  basin  in  Secs.  25,  36,  Twp.  10  N.,  Ra. 

5  E.,  in  the  Mexican  Highland  section,  Basin  and 
Range  province,  in  the  northernmost  extension  of 
the  Manzano  Mountains.  The  rock  is  the  Pennsyl¬ 
vanian  Magdalena  formation,  composed  of  limestone, 
feldspathic  sandstone,  breccia,  and  thin  shale  beds. 

An  anticlinal  axis  passes  north-south  through  the 
western  part  of  the  basin,  from  which  beds  dip  gently 
northeast  and  northwest.  About  97%  of  the  basin  is 
forest  covered,  2%  is  grassland,  and  the  remainder 
is  bare.  Grazing  is  light  and  the  range  is  generally 
very  poor,  with  only  small  areas  of  fair  range.  The 
main  channel  is  gullied  near  its  mouth  to  a  depth 
of  five  feet.  All  other  channels  show  signs  of  intense 
recent  scour,  down  to,  but  not  into  bedrock,  but  are 
not  actually  gullied.  Slopes  give  evidence  of  recent 
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sheet-wash  or  rill  erosion,  but  have  very  few  gul¬ 
lies.  A  single  gully  has  become  established  down  a 
ridge  crest,  erosion  having  occurred  along  a  trail. 
The  P-E  index  is  estimated  to  be  36  to  40,  depend¬ 
ing  on  elevation  and  exposure.  The  distribution  of 
rainfall  is  assumed  to  be  the  same  as  for  the 
Tijeras  Ranger  Station,  three  miles  northwest, 
where  during  the  months  of  May,  July,  August,  and 
September  47%  of  the  annual  precipitation  occurs. 

Tributaries  to  Bear  Canyon  (Twin  Sisters  Quad¬ 
rangle)  are  three  fourth-order  basins  in  Secs.  13, 

14,  23,  24,  Twp.  16  S.,  Ra.  13  W.,  in  the  Mexican 
Highland  section,  Basin  and  Range  province.  The 
area  is  a  rugged,  dissected  section  of  lavas,  mainly 
andesite  and  basalt  (Paige,  1916),  alternating  with 
gravels  and  tuffs;  the  lava  beds  form  prominent 
steps  in  the  hillside.  The  area  was  not  visited  in 
the  field.  The  estimated  P-E  index  is  60.  The  near¬ 
est  source  of  weather  records  is  Pinos  Altos, 
where  50%  of  the  yearly  precipitation  occurs  during 
July,  August,  and  September. 

Little  Cherry  Creek  (Twin  Sisters  Quadrangle) 
is  a  fifth-order  basin  in  Secs.  15,  16,  19,  20,  21,  22, 
Twp.  16  S.,  Ra.  13,  W.,  in  the  Mexican  Highland 
section  of  the  Basin  and  Range  province.  The  esti¬ 
mated  P-E  index  is  51.  In  other  respects,  the  basin 
is  similar  to  the  others  in  Twin  Sisters  quadrangle. 
The  basin  was  not  inspected  in  the  field. 

Mill  Creek  (Twin  Sisters  Quadrangle)  is  a  fourth- 
order  basin  similar  to  the  two  preceding  ones,  ex¬ 
cepting  the  P-E  index,  which  is  estimated  at  45. 

Wilson  Creek  Basins  (Twin  Sisters  Quadrangle) 
are  four  fourth-order  basins  similar  to  the  three 
other  areas  in  Twin  Sisters  quadrangle,  except  that 
the  estimated  P-E  index  is  51.  Forest  covers  100% 
of  the  basin. 


BASINS  IN  UTAH 

Dry  Hollow  (Lehi  Quadrangle)  is  a  fourth-order 
basin  in  Secs.  20,  21,  Twp.  4  S.,  Ra.  1  E.,  on  the 
south  side  of  the  Traverse  Mountains,  five  miles 
from  the  foot  of  the  Wasatch  Front,  in  the  Great 
Basin  section,  Basin  and  Range  province.  The  Tra¬ 
verse  Mountains  have  not  been  affected  by  Recent 
uplift,  as  have  the  Wasatch  Mountains,  resulting  in 
less  rugged  topography  (Marsell,  1953,  p.  20).  The 
prevailing  rock  type  is  hard  buff  quartzite.  The  soil 
contains  considerable  clay  that  is  probably  wea¬ 
thered  from  limestone  and/or  lava,  though  these 
lithologies  are  not  seen  in  the  basin.  Trees  and 
brush  cover  60%  of  the  area,  grass  35%,  and  5%  is 
bare.  The  range  which  is  generally  poor,  is  heavily 
over-grazed  in  places  by  sheep  and  cattle.  While 
a  few  small  tributaries  are  gullied,  the  main  chan¬ 
nel  is  not  and  is  indistinct  in  most  places.  The 
tributaries  usually  end  in  alluvial  cones,  and  do  not 
directly  join  the  main  channel.  The  entire  surface 


of  the  basin  is  covered  by  quartzite  fragments, 
which  cause  instability  of  the  slopes  on  the  steep 
hillsides.  The  estimated  P-E  index  is  48  to  52,  de¬ 
pending  on  elevation.  The  distribution  of  precipita¬ 
tion  is  probably  nearly  the  same  as  that  of  Alpine, 
five  miles  east,  where  except  for  four  dry  months, 
June  through  September,  rain  is  distributed  evenly 
throughout  the  year.  Only  21%  of  the  annual  precipi¬ 
tation  occurs  in  those  months. 

Hog  Hollow  (Lehi  Quadrangle)  is  a  fifth- order 
basin  similar  to  Dry  Hollow,  but  is  nearer  the 
Wasatch  Front.  The  surface  of  the  basin  shows  con¬ 
siderable  andesitic  lava  and  much  less  quartzite. 

A  greater  percentage  of  this  area  is  covered  by 
trees  and  bushes  than  in  Dry  Hollow.  The  estimated 
P-E  index  is  51  to  55,  depending  on  elevation.  The 
distribution  of  rainfall  is  the  same  as  for  Alpine, 
which  is  1  1/2  miles  southeast. 

Maple  Hollow  and  Mercer  Hollow  (Lehi  Quad¬ 
rangle)  are  two  fourth- order  basins  located  in 
Secs.  15  and  16,  21,  respectively,  Twp.  4  S.,  Ra.  1 
E.,  and  are  similar  to  the  two  preceding  areas  in 
this  quadrangle,  with  a  greater  amount  of  lava 
cropping  out  than  in  Dry  Hollow.  The  estimated 
P-E  index  is  50  to  53,  depending  on  elevation.  The 
distribution  of  precipitation  is  assumed  to  be  the 
same  as  above  for  Alpine. 

Pole  Canyon  (Spanish  Fork  Peak  Quadrangle)  is 
a  large  fifth-order  basin  in  Secs.  13,  14,  15,  22, 

23,  24,  27,  Twp.  9  S.,  Ra.  3  E.,  in  the  Wasatch 
Mountains  of  the  Middle  Rocky  Mountains  province. 
Almost  continual  uplift  during  the  late  Pleistocene 
and  Recent  has  caused  rejuvenation,  and  accounts 
for  the  great  ruggedness  of  the  area.  The  main  lith¬ 
ologies  are  limestone  and  quartzite.  Forest  covers 
almost  the  entire  area.  The  estimated  P-E  index  is 
70  to  80,  depending  on  elevation.  The  distribution 
of  precipitation  is  assumed  to  be  similar  to  that  of 
Thistle,  where  the  rainfall  is  fairly  evenly  distrib¬ 
uted  throughout  the  year,,  with  slight  maxima  in  Jan¬ 
uary  and  March,  and  a  minimum  in  June  and  July. 

Deadmans  Hollow  (Springville  Quadrangle)  is  a 
small  third-order  basin  in  Secs.  1  and  36,  Twps.  7 
and  8  S.,  Ra.  3  E.,  in  the  Wasatch  Mountains  of  the 
Middle  Rocky  Mountains  province.  Uplift  has  caused 
rejuvenation  and  rugged  topography.  The  main 
rock  type  is  buff  quartzite  and  sandstone.  Trees 
and  brush  cover  about  20%  of  the  area,  the  remain¬ 
der  has  sparse  bunch  grass  or  is  completely  bare. 
The  estimated  P-E  index  is  45  to  50,  depending  on 
elevation.  The  distribution  of  rainfall  is  probably 
the  same  as  for  Provo  Bench  station,  where  the 
occurrence  of  precipitation  is  fairly  even  through¬ 
out  the  year,  with  a  slight  maximum  in  May,  and  a 
minimum  in  June  and  July. 

Spring  Creek  (Springville  Quadrangle)  is  a  fourth- 
order  basin  in  Secs.  24,  25,  26,  Twp.  7  S.,  Ra.  3  E., 
in  the  Wasatch  Mountains  of  the  Middle  Rocky 
Mountains  province.  Uplift  of  the  mountain  range 
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has  caused  rejuvenation  and  rugged  topography  in 
the  area.  The  main  rock  type  is  gray  limestone 
and  buff  quartzite.  Trees  and  bushes  cover  85%  of 
the  basin,  grass  nearly  15%,  while  only  a  negligi¬ 
ble  part  is  bare.  Cattle  graze  the  basin  heavily 
near  the  roads  and  the  main  channel;  the  range  is 
poor  to  fair.  The  main  channel  is  gullied  to  a  depth 
of  3  feet  near  its  head  and  20  feet  at  the  mouth. 

The  P-E  index  is  estimated  at  70  to  75,  depending 
on  elevation.  The  distribution  of  rainfall  is  prob¬ 
ably  the  same  as  above  for  Provo  Bench  station. 


Tributary  to  Hobble  Creek  (Springville  Quad¬ 
rangle)  is  a  small  third-order  basin  in  Secs.  36 
and  1,  Twps.  7  S.  and  8  S.,  Ra.  3  E.  The  basin  is 
similar  to  Deadmans  Hollow.  The  main  lithology  is 
buff  sandstone,  quartzite,  and  gray  limestone.  The 
basin  is  grazed  heavily,  and  the  range  is  poor  to 
fair.  A  few  new  gullies  are  present  on  the  basin 
slopes,  but  old  channels  are  not  gullied.  The  esti¬ 
mated  P-E  index  is  55  to  60.  Precipitation  distribu¬ 
tion  is  the  same  as  for  Provo  Bench  station,  the 
nearest  source  of  weather  data. 
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Methods  of  Statistical  Analysis 


APPROPRIATENESS  OF  STATISTICAL  METHODS 

The  important  assumptions  inherent  in  statisti¬ 
cal  procedures  are  becoming  well-known.  Depart¬ 
ures  in  practice  from  rigid  adherence  to  the  condi¬ 
tions  imposed  by  those  assumptions  will  produce 
more-or-less  serious  effects  on  validity  of  conclu¬ 
sions  drawn  from  the  results  of  a  test.  In  this  study, 
the  attempt  was  made  to  comply  with  sound  statis¬ 
tical  methods,  but  some  wide  departures  exist. 

Normality  of  parent  populations  is  fundamental 
in  the  parametric  methods  of  analysis  used  here. 
The  number  and  size  of  samples  taken  in  any  single 
category  of  measurement  within  a  given  environ¬ 
ment  are  insufficient  to  allow  detection  of  moderate 
departures  from  normality.  One  exception  is  max¬ 
imum  valley- side  slope,  which  other  investigators 
have  found  to  be  normally  distributed  within  a 
drainage  basin  (Strahler,  1950,  p.  681-685);  this  is 
also  true  of  all  areas  studied  here.  Another  ex¬ 
ception  is  the  drainage  density  of  subordinate  ba¬ 
sins  within  a  single  larger  basin,  or  within  the 
same  rock  type;  these  distributions  do  not  appear 
to  be  significantly  different  from  normal,  although 
total  samples  are  small.  Miller  (1953,  p.  14)  ren¬ 
dered  symmetrical  the  distribution  of  areas  of 
first-order  basins  by  logarithmic  transformation. 

Sample  means  tend  to  be  normally  distributed, 
regardless  of  the  parent  distribution,  when  samples 
are  large,  and  the  parent  populations  have  finite 
variances  (Dixon  and  Massey,  1951,  p.  52).  Most 
data  used  here  are  sample  means,  and  for  that  rea¬ 
son  are  expected  to  be  approximately  normally  dis¬ 
tributed.  The  analysis  of  variance  and  similar  tests, 
are  remarkably  robust  to  departures  from  normal¬ 
ity  (Box,  1953,  p.  318)  for  tests  on  means.  The 
question  of  normality  is  therefore  probably  of  less 
importance,  for  the  purposes  of  this  study,  than 
other  factors. 

Objective  sampling  is  the  sine  qua  non  of  all 
statistical  methods;  this  requirement  of  statistics 
imposes  restrictions  on  the  methods  of  obtaining 
data  suitable  for  analysis.  Choice  of  sampling 
methods  is  the  area  of  greatest  divergence  from 
conventional  methods  used  in  geology.  Practicing 
statisticians  have  long  realized  the  difficulty  in  ob¬ 
taining  unbiased  samples  of  a  population  when  there 
is  freedom  of  choice  in  selecting  the  observations 
to  be  recorded.  The  training  of  a  geologist  is  at 
least  in  part  that  of  learning  to  observe  in  a  man¬ 
ner  that  assumes  a  bias  toward  representation  of 


the  true  situation,  thus  reducing  variability  in 
measurements  and  allowing  treatment  by  non- 
statistical  methods.  In  the  statistical  approach,  it 
must  be  assumed  that  the  true  characteristics  of  the 
universe  are  unknown  and  can  best  be  estimated 
from  random  samples.  Attempts  to  obtain  data  close 
to  some  value  arrived  at  by  “judgment,”  “insight,” 
or  “experience,”  will  nearly  always  result  in 
biased  samples,  and  failure  of  the  method  of  analy¬ 
sis.  In  this  study,  samples  were  taken  both  by  ran¬ 
dom  and  non-random  methods,  and  results  are  ac¬ 
cordingly  subject  to  qualification. 

Study  areas  were  chosen  after  cursory  examina¬ 
tion  showed  them  to  be  satisfactory  in  accessibil¬ 
ity,  absence  of  cliffs,  low  in  roughness,  and  in  the 
mature  stage  of  development.  No  effort  was  made 
to  select  by  random  methods  sample  basins  in  all 
types  of  rock;  the  majority  are  in  granite  and  coarse 
clastic  sediments,  and  few  in  shale  or  limestone. 
Measurements  within  basins  were  largely  on  areas 
selected  to  be  representative  and,  as  noted  above, 
certain  measurements  would  mean  little  on  areas 
differing  greatly  from  the  average  of  the  entire 
basin.  Selecting  measurement  sites  in  a  manner  to 
reduce  the  variance  of  the  observations  may  not 
bias  the  observed  mean.  However,  if  analysis  of 
variance  is  used  to  test  for  significance  of  differ¬ 
ences  among  means  of  a  measurement  in  several 
basins,  too  many  significant  results  will  be  obtained 
if  inferences  are  made  about  the  entire  population 
instead  of  the  portion  that  was  measured.  This  is 
because  the  error  attributable  to  a  single  measure¬ 
ment  will  occur  in  the  denominator  of  the  F  ratio, 
and  will  cause  it  to  exceed  the  critical  value  too 
frequently.  In  regression  analysis,  selection  in  in¬ 
dependent  variables  probably  does  not  seriously 
affect  the  value  of  correlation  and  regression  co¬ 
efficient  unless  either  extreme  or  intermediate 
values  are  eliminated.  In  this  case,  as  extreme 
values  were  frequently  rejected,  the  indicated  cor¬ 
relation  will  be  too  low.  Selection  in  the  dependent 
variables  has  greater  effect.  If  the  selection  tends 
to  reduce  the  variability  of  the  dependent  variable, 
the  indicated  correlation  will  be  too  high  (Ezekiel, 
1930,  p.  266-267). 

Additivity  and  linearity  of  effects  is  an  impor¬ 
tant  assumption  in  multiple  regression  analysis, 
when  several  variables  are  examined  for  their  in¬ 
fluence  on  another.  The  plane  of  best  fit  through  a 
set  of  points  can  always  be  found,  and  in  cases  pre¬ 
sented  here  it  is  used  only  where  the  dependent 
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variable  is  linear  when  plotted  against  each  inde¬ 
pendent  variable.* *  It  is  not  definitely  known  whether 
or  not  effects  of  rock  type,  vegetal  cover,  and  cli¬ 
matic  elements  on  morphologic  properties  are  addi¬ 
tive.  However,  as  will  be  seen  in  considering  effects 
of  these  factors  on  valley- side  slope  and  drainage 
density,  non-linear  two-variable  plots  obtained 
were  rectified  by  taking  logarithms,  indicating 
either  multiplicative  effects  or  non-linear  variables. 

Multiple  regression  and  correlation  can  be  used 
where  Y  is  a  dependent  random  variable  and  X1? 

X2,  ...  Xn  are  independent  random  variables  whose 
values  are  known,  and  not  subject  to  errors  in 
measurement  or  observation.  In  the  cases  to  be 
considered  here,  values  of  all  variables  were  ob¬ 
tained  by  measurement  in  the  field,  from  maps,  or 
by  meterological  observation,  and  hence  are  esti¬ 
mates  subject  to  at  least  some  observational  error. 
The  effect  of  introducing  observational  error  into 
the  independent  variables  is  to  reduce  the  indicated 
correlation  and  regression  coefficient.  The  regres¬ 
sion  is  Y  on  estimated  Xi,  ...  Xn,  rather  than  on 
true  X1?  ...  Xn  (Ezekiel,  1930,  p.  225;  Professor 
Howard  Levene,  Personal  Communication).  The 
combined  influence  of  selection,  tending  to  increase 
indicated  correlation  and  regression  coefficients, 
and  introduction  of  error  into  the  X^,  tending  to 
reduce  the  coefficients,  is  complex,  and  the  effects 
on  the  conclusions  drawn  in  this  study  are  not  known. 

The  theoretical  considerations  discussed  in  the 
previous  paragraphs  are  of  great  importance  in 
confirming  hypotheses  by  statistical  methods,  when 
research  and  experimentation  have  been  properly 
controlled,  when  properties  of  experimental  mate¬ 
rial  are  well-known,  and  when  only  a  few  variables 
are  considered  at  one  time.  The  research  in  this 
study  is  primarily  investigative,  and  it  is  the  wri¬ 
ter’s  opinion  that  it  is  probably  sound  scientific 
practice  to  lean  toward  declaring  too  many  rela¬ 
tions  or  differences  to  be  significant,  rather  than 
risk  failing  to  detect  a  truly  significant  relation.  It 
should  be  understood  that  the  hypotheses  thus  pro¬ 
posed  are  not  considered  to  be  definitely  confirmed 
until  rigorous  investigation  has  produced  sustaining 
evidence. 


classes  and  transferred  to  punch  cards.*  The  data 
thus  recorded  were: 


Morphometric 

properties 

Drainage  density 
Density  of  first-order 
channels 

Channel  frequency 
Frequency  of  first- 
order  channels 
Ratio  of  channel  length 
to  perimeter 
Valley-wall  slope 
Ruggedness  number 
Length  of  perimeter 
Inflection  angle 
Relative  relief 
Basin  order 


Properties  of  soil 
and  cover 

Infiltration  capacity 
Percent  bare  area 
Soil  strength  (wet) 
Roughness  number 
Cover 
Rock  type 
Life  zones 


Climatic  properties 

5-year  1-hour  rainfall 
P-E  index 


The  information  available  for  each  basin  was  re¬ 
corded  on  two  series  of  punch  cards.  Fourth  and 
fifth-order  basins  were  divided  into  third-order 
basins,  with  one  card  for  each.  Cards  were  then 
easily  sorted  into  classes  to  obtain  frequencies  in 
each  cell  in  a  two-way  table.  Usual  methods  of  an¬ 
alysis  of  grouped  data  yielded  correlation  and 
regression  coefficients  with  comparative  speed. 

Of  the  total  number  of  properties  investigated, 
those  proving  to  be  significantly  correlated  were 
considered  for  further  analysis. 

Correlation  analysis  of  grouped  data  is  limited 
by  several  considerations:  (1)  because  of  grouping, 
the  estimated  standard  errors  tend  to  be  too  large, 
and  the  correlation  coefficient  too  small.  It  is  nec¬ 
essary  to  have  a  considerable  number  of  cells  in 
the  correlation  table,  and  this  requirement  is  lim¬ 
ited  by  the  number  of  holes  in  the  cards.  (2)  The 
correlation  table  is  difficult  to  interpret  as  to 
linearity,  degree  of  correlation,  importance  of  ex¬ 
treme  values,  etc.  (Croxton,  1953,  p.  136).  Despite 
these  objections,  the  method  of  correlation  using 
grouped  data  on  punched  cards  proved  valuable  in 
facilitating  rapid  analysis.  Appendix  V  gives  the 
correlation  tables  and  the  correlation  coefficients 
obtained,  for  a  number  of  important  relations. 


REGRESSION  AND  CORRELATION  ANALYSIS  OF 
PAIRED  VARIABLES 


FINAL  MULTIVARIATE  ANALYSIS 


Preliminary  analysis  of  grouped  data:  The  data 

of  greatest  interest  and  promise  were  grouped  into 

*In  view  of  the  differences  in  opinion  prevalent  on  the  ef¬ 
fects,  say,  of  increased  rainfall  on  texture  of  topography, 
when  modifications  of  intensity-frequency  relations,  soil 
properties,  vegetal  cover,  and  stream  regimen  are  con¬ 
sidered,  the  question  of  non -additivity  of  effects  should 
be  the  subject  of  further  research. 


Certain  of  the  significant  relations  discovered 
from  analysis  of  grouped  data  were  further  ana¬ 
lyzed  by  ungrouped  methods,  in  order  to  find  the 
statement  of  most  linear  function  giving  the  best 
fit  and  lowest  scatter.  The  final  results  of  these 
are  given  in  part  B  of  Appendix  V.  These  relations 

*Keysort  cards,  manufactured  by  The  McBee  Company, 
Athens,  Ohio. 
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are  valid  as  they  are,  but  their  usefulness  can  be 
increased  if  combined  with  other  measurements. 
Several  combinations  of  logarithmic  and  arithmetic 
scales  were  tried  for  each  pair  of  variables. 

Multivariate  correlation  analysis  is  a  method  of 
analyzing  the  effect  of  several  independent  variables 
on  the  variation  of  a  dependent  variable.  This  meth¬ 
od  is  necessary  in  determining  whether  each  in¬ 
dependent  variable  acts  independently  of  the  others, 
or  if  an  apparent  significant  correlation  obtained 
on  a  two-variable  analysis  is  actually  due  to  their 
both  being  correlated  with  a  third  variable.  Where 
all  aspects  of  an  experiment  are  not  easily  con¬ 
trolled  and  all  possible  combinations  of  factors  are 
not  found,  multiple  correlation  and  regression  anal¬ 
ysis  will  be  of  greater  use  in  quantitative  geomor¬ 
phology,  than  the  standard  experimental  designs 
employing  the  analysis  of  variance.  (See  analysis 
of  covariance.) 

The  methods  used  in  analyzing  two  to  four  vari¬ 
ables  are  usually  presented  fully  in  elementary 
statistics  textbooks.  The  amount  of  calculation 
needed  to  prepare  a  multiple  correlation  study  in¬ 
creases  as  some  power  of  the  number  of  variables 
considered.  In  this  study,  the  final  analysis  con¬ 
sists  of  determination  of  the  effects  on  a  dependent 
variable  (drainage  density),  of  five  independent 
variables  (infiltration  capacity,  P-E  index,  etc.). 

No  explicit  solution  of  normal  equations  for  this 
many  variables  is  available.  The  normal  equations 
were  first  expressed  in  matrix  form,  the  inverse 
matrix  found  by  the  Doolittle  method,  and  the  in¬ 
verse  then  used  to  obtain  the  desired  solutions. 

(See  Anderson  and  Bancroft,  1952,  p.  168,  and 
Ezekiel,  1950,  pp.  168,  363,  367).  The  full  solution 
is  given  in  Appendix  VII. 

Correlation  analysis  of  six  variables,  to  yield 
stable  solutions,  usually  requires  more  values  on 
each  variable  than  are  used  here.  In  this  study, 
however,  each  of  the  23  values  used  represents 
more  actual  measurements  averaged  to  obtain  it. 
Drainage  density  represents  the  average  stream 
length  per  unit  of  area,  and  was  obtained  from  ten 
to  300  individual  measurements  of  stream  length; 
each  P-E  index  is  obtained  from  12  paired  values; 
infiltration  capacity  is  an  average  of  four  to  11 
readings;  percent  bare  area  is  the  average  of  two 


to  20  measurements;  roughness  number  is  meas¬ 
ured  from  hundreds  of  rock  fragments.  Much  better 
results  would  probably  be  obtained  by  using  the  or¬ 
iginal,  unaveraged  figures,  as  the  gain  in  degrees  of 
freedom  would  be  considerable  in  most  cases.  Un¬ 
fortunately,  to  do  this  would  mean  duplicating  values 
of  the  dependent  variable.  Therefore,  the  value  of 
the  drainage  density  would  have  to  be  repeated  as 
many  times  as  there  are  infiltration  readings  in  the 
basin,  and  there  would  be  no  real  gain  in  degrees 
of  freedom.  A  solution  would  be  the  use  of  analysis 
of  variance  with  unequal  replication  in  each  cell, 
but  this  is  unfeasible  because  all  combinations  of 
the  five  independent  variables  are  not  found. 

Analysis  of  variance  is  a  technique  of  determin¬ 
ing  the  significance  of  difference  among  several 
means  drawn  from  distinct  environments.  This  study 
treats  the  effect  of  rock  type  on  several  morpho¬ 
metric  and  soil-cover  properties  of  the  basins  in 
each  type.  A  simple  one-way  classification  is  used, 
and  six  rock  categories  are  distinguished.  Where  sig¬ 
nificant  differences  are  found  to  exist  among  the 
means,  Scheffe’s  method  of  judging  contrasts  is  used 
to  test  specific  pairs  or  groups  of  means  for  signifi¬ 
cant  differences  (Scheffe,  1953,  p.  87).  The  analysis 
of  variance  tables  are  presented  in  Appendix  VI, 
and  contrasts  are  analyzed,  where  appropriate. 

The  analysis  of  covariance  is  essentially  similar 
to  the  analysis  of  variance,  with  the  addition  that 
an  associated  measurement  is  made  for  each  item, 
in  order  to  determine  how  much  of  the  differences 
in  the  means  of  the  primary  measurements  is  the 
result  of  variation  in  the  related  variable.  The 
means  are  adjusted  to  what  they  would  be  were  the 
associated  measure  the  same  for  each.  Then,  if  sig- 
nigicant  differences  are  present,  these  are  due  to 
the  differences  in  the  main  variable  of  classification 
(Dixon  and  Massey,  1951,  p.  173).  One  aim  of  the 
study  was  to  determine  if  a  consistent  difference  in 
drainage  density  could  be  shown  to  exist  among  ba¬ 
sins  with  gullied  slopes,  trenched  channels,  and  no 
gullies.  Infiltration  capacity  was  the  associated 
measurement,  and  variations  in  drainage  density 
due  to  differences  in  infiltration  capacity  were  re¬ 
moved.  The  analysis  of  covariance  table  is  given  in 
Appendix  VI  B,  and  the  results  discussed  in  the  con¬ 
clusions  section. 
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Inferences  Drawn  From  Results 
of  Statistical  Tests 


RESTRICTIONS  ON  INFERENCES 

Certain  restrictions  must  be  placed  on  the  con¬ 
clusions  because  of  the  nature  of  the  statistical 
methods  used.  Perhaps  most  important  is  the  in¬ 
fluence  of  subjective  sampling  on  the  values  of  the 
correlation  coefficients  obtained  in  the  analysis. 
This  influence  may  be  great,  but  little  can  be  said 
of  it  before  much  more  detailed  work  is  done  on 
limited  aspects  of  this  general  study.  Rather  high 
values  for  the  critical  level  (e.g.  .05,  or  .10)  were 
used,  so  that  by  chance  alone,  one  would  expect 
some  of  the  differences  or  the  correlations  to  be 
sufficiently  great  to  be  judged  significant  erron¬ 
eously. 


RESULTS  OF  ANALYSIS  OF  VARIANCE  AND 
COVARIANCE  TESTS 

Rock  type  and  drainage  density:  The  type  of 
rock  has  an  apparent  influence  on  drainage  density, 
according  to  the  high  value  obtained  for  F  in  the 
following  analysis  of  variance  table: 


Source  of 
Variation 

S.  S. 

d.f. 

M.  S.  F 

Among  rock 
types 

Within  rock 

57,371 

5 

11,474  56.7 

types 

12,140 

60 

202.3 

Total 

69,511 

65 

1,069 

F.001=  4.76 


Six  categories  of  lithology  are  distinguished:  clastic 
(sandstone  and  alluvium),  shale,  limestone,  schist, 
granite  (includes  gneiss),  and  acid  volcanic  (in¬ 
cludes  pumice).  Inspection  of  the  means  of  D  for 
each  rock  type  shows  that  some  are  decidedly  dif¬ 
ferent  from  others  (Appendix  VI  A).  The  estimated 
effect  of  each  lithology  on  D  is: 

Shale:  150.5  Schist:  37.9 

Clastic:  -  7.8  Granitic:  -  3.86 

Limestone:  -12.5  Acid  volcanic:  -12.28 

These  values  are  probably  meaningless,  and  it 
should  be  noted  that  they  do  not  add  to  zero,  as  the 


sample  size  was  different  in  each  case.  The  impor¬ 
tant  point  is  that  shale  and  schist  have  average 
drainage  density  well  above  the  mean,  whereas 
limestone  and  acid  volcanic  rocks  are  well  below 
the  mean. 

Scheffe’s  method  of  judging  contrasts  was  used 
to  find  which  of  the  means  produce  the  significant 
results  (Scheffe  ,  1953).  Significant  differences  in 
drainage  density  exist  between  the  average  for  shale 
and  schist,  the  two  lithologies  with  the  highest 
values,  and  the  average  for  granitic,  clastic,  vol¬ 
canic,  and  limestone,  and  also  between  the  average 
for  shale  and  the  average  for  schist.  No  significant 
difference  in  drainage  density  was  found  among 
granitic,  volcanic,  clastic,  or  limestone  (see  Appen¬ 
dix  VI  A). 

It  is  apparent  that  rock  type  has  perceptible  ef¬ 
fect  on  drainage  density  under  the  general  mountain- 
type  climate  in  this  study.  However,  the  sampling 
of  shale  and  limestone  basins  is  poor,  and  all  ba¬ 
sins  from  both  these  lithologies  are  taken  from 
areas  with  extremes  of  other  conditions.  For  in¬ 
stance,  all  the  shale  basins  lie  in  a  very  dry  cli¬ 
mate.  The  limestone  basins  are  in  the  Wasatch 
Range,  and  were  subjected  to  late  Pleistocene  up¬ 
lift  and  rejuvenation.  Tests  of  differing  rock  types 
within  the  same  general  climatic  type,  and  in  areas 
with  identical  histories  by  Miller  (1953,  p.  24,  and 
Table  1)  show  the  importance  of  lithology  in  deter¬ 
mining  drainage  density  in  a  humid  climate.  Further 
such  studies  may  show  that  elimination  of  variation 
due  to  climate  and  differing  geologic  histories  will 
reduce  the  error  variance  enough  to  allow  detection 
of  a  significant  difference  among  granite,  clastic 
rocks,  and  acid  volcanic  rocks,  that  in  this  study 
show  no  differences. 

Rock  type  and  infiltration  capacity:  Individual 
infiltration  capacity  readings  (not  basin  averages) 
were  tested  against  rock  type.  The  rock  types  dis¬ 
tinguished  are  the  same  as  for  drainage  density. 

The  analysis  of  variance  table  is: 


Source  of 
Variation 

S.  S. 

d.f. 

M.  S. 

F 

Among  rock  types 
Within  rock  types 

67.89 

249.4 

5 

135 

13.58 

1.85 

7.43 

Total 

317.3 

140 

2.27 

F  001(5,135)  =  4.10 
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The  value  of  F  is  significant,  and  indicates  that 
differences  in  infiltration- capacity  means  exist 
among  the  various  rock  types.  Examination  of  the 
rock-type  means  (Appendix  VI  A)  shows  that  some 
differ  more  widely  from  the  grand  mean  than 
others.  The  estimated  effect  of  each  lithology  on 
infiltration  capacity  is: 

Shale:  -1.43  Schist:  -1.00 

Clastic:  -0.53  Granitic:  0.81 

Limestone:  0.52  Acid  volcanic:  -0.31 

Much  reliance  cannot  be  placed  on  these  values.  It 
is  significant,  perhaps,  that  shale  and  schist  have 
infiltration  capacities  below  the  mean,  whereas 
they  have  drainage  densities  well  above  the  mean. 
Limestone  and  granitic  rocks  have  infiltration 
capacities  above  the  mean,  but  drainage  density 
below  the  mean.  Acid  volcanic  rocks  probably  have 
infiltration  capacity  not  significantly  different  from 
the  mean,  though  drainage  density  below  the  mean, 
and  thus  act  somewhat  anomalously. 

The  high  value  of  F  is  the  result  of  the  high  in¬ 
filtration  rates  on  granitic  and  limestone  litholo¬ 
gies.  The  contrast  (granitic  and  limestone)  vs. 
(clastic  and  volcanic)  is  significant  at  .05.  Also, 
the  contrast  (clastic  and  volcanic)  vs.  (shale  and 
schist),  formed  from  the  groups  with  the  lowest  in¬ 
filtration  rates,  was  judged  not  significant  at  .05 
(see  Appendix  VI  A). 

The  question  arises  as  to  whether  lithology 
affects  values  of  drainage  density  other  than 
through  the  value  of  the  infiltration  capacity.  A 
possible  method  of  approach  may  be  through  an  ex¬ 
periment  using  analysis  of  covariance,  correcting 
values  of  drainage  density  for  differences  in  infil¬ 
tration  capacity,  and  selecting  basins  with  the 
same  general  climate  and  history,  but  from  dif¬ 
ferent  lithologies.  The  data  obtained  in  this  study 
are  inadequate  to  attempt  this  analysis. 

Rock  type  and  valley-side  slopes:  The  individ¬ 
ual  slope  measurements  within  each  area  studied 
were  arranged  into  groups  corresponding  to  rock 
type,  without  regard  to  basin.  The  estimated  effect 
of  each  lithology  on  valley-slide  slope  is: 


Shale: 

7.9 

Schist: 

-4.4 

Clastic: 

1.5 

Granitic: 

-1.3 

Limestone: 

9.1 

Acid  Volcanic: 

-3.9 

Limestone  is  well  above  the  general  mean,  prob¬ 
ably  because  the  only  samples  are  from  basins  in 
the  Wasatch  Mountains  which  have  undergone  late 
Pleistocene  rejuvenation.  Shale  is  also  high,  per¬ 
haps  because  of  the  aridity  of  the  region  sampled. 
Schist  and  acid  volcanic  rocks  are  low.  To  test 
whether  these  differences  are  significant,  an  anal¬ 
ysis- of- variance  test  was  made: 


Source  of 

Variation  S.  S.  d.  f.  M.  S.  F 


Among  rock 

types  11,983.48  5  2396.69  18.99 

Within  rock 

types  110,298.35  874  126.19  . 


Total  122,281.83  879 


F001(5,874)  =  4.10 

The  value  of  F  obtained  is  significant  at  .001,  and 
the  evidence  points  to  the  conclusion  that  rock  type 
affects  valley-side  slopes. 

Because  the  average  slopes  on  shale  and  lime¬ 
stone  are  10°  greater  than  on  the  other  lithologies, 
it  is  reasonably  suspected  that  the  high  values  are 
the  result  of  some  other  factors,  specifically  arid¬ 
ity  and  rejuvenation,  that  become  manifest  through 
the  poor  sampling  of  shale  and  limestone  basins, 
mentioned  above.  Therefore,  the  contrast  (clastic 
and  granitic)  vs.  (schist  and  volcanic)  was  analyzed, 
and  was  found  to  be  significant  at  .05.  It  can  be 
stated  that  lithology  probably  does  influence  valley- 
side  slopes  (see  Appendix  VI  A). 

Infiltration  capacity,  drainage  density,  and 
gullies:  Basins  for  which  infiltration- capacity  read¬ 
ings  are  available  were  grouped  according  to  whe¬ 
ther  they  had  gullied  slopes,  trenched  channels,  or 
unmodified  drainage  lines.  Infiltration  capacity  was 
used  as  a  supplementary  measurement  to  adjust  the 
drainage  densities  of  the  basins.  The  purpose  of  the 
test  was  to  determine:  (1)  whether  basins  with  gul¬ 
lies  show  consistently  higher  drainage  densities 
than  basins  with  trenched  channels  or  unmodified 
drainage  lines,  when  variation  due  to  infiltration 
capacity  is  removed;  (2)  whether  basins  with  gullied 
slopes  show  consistently  lower  infiltration  capacity. 
The  analysis  of  covariance  table  is: 


Source  of 


Variation 

d.f. 

£x2 

£xy 

sy2 

d.f. 

sy’2 

M.S. 

Among 

means 

2 

1.4 

-.74 

.42 

2 

.50 

.25 

Within 

groups 

20 

10.5 

-12.03 

24.44 

19 

10.68 

.56 

Total 

22 

11.9 

-12.77 

24.86 

21 

11.18 

For  X,  F 

=  1.33;  for  Y,  F  = 

0.17; 

for  YT,  F  = 

0.45, 

where  X  is  infiltration  capacity,  Y  unadjusted 
drainage  density,  and  Y’  the  drainage  density  ad¬ 
justed  for  differences  in  infiltration  capacity.  The 
critical  value  of  F  05(2,2O)  is  3.49,  and  for  F 05 (2, 19) 
is  3.52.  None  of  the  values  of  F  obtained  is  signifi¬ 
cant.  Though  there  is  an  increase  in  F  from  Y  to 
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Yf,  it  is  not  enough  to  exceed  the  critical  value.  As 
can  be  seen  from  the  effects  (see  also  Appendix 
VI  B): 


Effects  of  gullying 
on  Loge  D 


Channels 

trenched:  -0.237 

Slopes  gullied:  0.167 

No  gullies  -0.035 


Effects  of  gullying 
on  infiltration 
capacity 

Channels 

trenched:  0.15 

Slopes  gullied:  -0.28 
No  gullies:  0.27 


The  mean  drainage  density  for  gullied  basins  is 
higher  than  for  trenched  unmodified  basins,  and  the 
mean  infiltration  capacity  is  lower  than  for  the 
other  classes;  basins  with  trenched  channels  have 
lower  drainage  density  than  either  gullied  basins 
or  untrenched  and  ungullied  basins.  However  these 
differences  are  not  great  enough  to  be  judged  sig¬ 
nificant  at  .05  considering  the  amount  of  variability 
in  the  samples.  These  results  fit  theoretical  con¬ 
siderations,  and  that  may  be  reason  to  believe  that 
with  larger  samples,  and  restricted  sampling,  the 
results  would  prove  significant. 

There  is  present  a  significant  decrease  in  within- 
groups  mean  square  (variance)  when  the  effect  due 
to  differences  in  infiltration  is  removed.  The  unad¬ 
justed  mean  square  is  24.44/20  =  1.22,  and  the 
adjusted  mean  square  is  0.56;  F  =  1.22/0.56  = 

2.18,  which  is  greater  than  the  critical  value  for 
.05.  This  indicates  that  a  significant  effect  of  in¬ 
filtration  capacity  on  drainage  density  is  present. 


RESULTS  OF  SIMPLE  CORRELATION  AND 
REGRESSION  ANALYSES 


Relations  among  morphometric  properties  of 
relief:  Preliminary  analysis  with  punched  cards 
indicated  that  significant  correlations  exist  among 
many  of  the  morphometric  properties  of  the  basins 
studied.  The  correlation  tables  thus  obtained  are 
given  in  Appendix  V  A.  Among  properties  involving 
relief,  it  was  found  that  maximum  valley-side  slopes 
are  correlated  positively  with  relative  relief  and 
ruggedness  number;  slopes  are  negatively  corre¬ 
lated  with  channel  frequency  and  drainage  density. 

A  scatter  plot  of  valley-side  slope  (ordinate) 
against  ruggedness  number  (abcissa)  showed  that  a 
large  number  of  points  representing  third-order 
basins  lay  above  the  average  line,  and  points  repre¬ 
senting  fifth-order  basins  lay  consistently  below 
the  average,  indicating  that  the  ruggedness  number 
was  affected  by  the  greater  absolute  relief  of  fifth- 
order  basins.  For  purposes  of  estimating  valley- 
side  slopes,  a  better  fit  is  obtained  by  decreasing 


the  value  of  H  by  an  amount  depending  on  the  basin 
order.  This  procedure  can  be  explained  by  consider¬ 
ing  a  model  area  having  uniform  drainage  density 
throughout,  and  with  uniform  valley-side  slopes, 
for  all  orders  of  basins.  Then  the  absolute  relief 
within  a  basin  will  normally  be  greater  for  basins 
of  higher  orders,  simply  because  they  are  larger. 
Whereas  the  drainage  density  remains  constant 
for  all  sizes  and  orders,  the  ruggedness  number 
will  increase  with  order.  Therefore,  to  obtain  the 
best  relation  with  valley-side  slope,  which  is  con¬ 
stant,  the  value  of  H  obtained  for  each  basin  should 
be  divided  by  some  number  that  increases  with 
order,  so  that  a  constant  for  the  area  is  obtained. 

In  this  case  the  order  itself  was  used,  though  some 
function  of  the  order  number  might  give  better  re¬ 
sults,  perhaps  u2  or  u^.  The  plots  of  0  against 
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log  —  and  log  R  are  given  in  Appendix  V  A.  The 

relation  of  valley- side  slopes  with  these  variables 
is  considered  below  in  connection  with  climate  and 
soil  properties. 

The  correlations  of  ©  with  log  F  and  log  D  indi¬ 
cate  that  an  inverse  relationship  exists  between 
slope  steepness  and  topographic  texture  (Appendix 
V  A)  The  absolute  value  of  the  correlation  co¬ 
efficient  of  ©  with  log  D  (-.43)  is  greater  than  for 
©  with  log  F  (-.34),  indicating  that  perhaps  part 
of  the  decrease  in  texture  with  increasing  valley- 
side  slopes  is  the  result  of  foreshortening  of 
steeper  channels  by  projection  onto  the  map.  Fore¬ 
shortening  should  not  affect  the  number  of  channel 
segments,  and  as  channel  frequency  is  negatively 
correlated  with  ©,  the  decrease  in  fineness  is 
probably  real.  A  possible  mathematical  model 
fitting  this  would  have  valley- side  slopes  inversely 
proportional  to  the  logarithm  of  both  drainage  den¬ 
sity  and  channel  frequency.  Conversely,  topographic 
texture  would  vary  as  a  negative  exponent  of  aver¬ 
age  valley-side  slope. 

In  addition,  D  and  F  are  negatively  correlated 
with  R,  as  would  be  expected,  since  R  is  highly 
correlated  with  ©  (Appendix  V  A).  Because  corre¬ 
lation  is  low,  it  is  not  possible  to  determine  which 
type  of  function  gives  the  best  linear  scatter  plot. 

Which  of  the  morphometric  properties,  among 
D,  F,  R,  and  ©  can  be  considered  as  controlling 
factors,  if  any  are,  is  not  definitely  known.  A 
plausible  hypothesis  is  that  R  and  either  F  or  D, 
or  both,  are  relatively  independent.  D  and  F  de¬ 
pend  on  rock  type  and  climate,  and  R  probably  de¬ 
pends  greatly  on  past  diastrophic  history,  Then 
H  H 

©  and  H  or  —  are  dependent.  Log  —  is  correlated 

positively  with  log  R,  (r  =  .666,  b  =  0.237),  as 
expected,  suggesting  the  mathematical  model  that 
ruggedness  number  corrected  for  basin  order  var¬ 
ies  about  as  the  one-fourth  power  of  relative  re¬ 
lief. 
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RELATIONS  AMONG  PLANIMETRIC  PROPER¬ 
TIES  OF  BASINS 


Simple  correlation  of  two  variables  by  grouped 
data  showed  that  high  correlations  exist  among 
logarithms  of  D,  D]_,  F,  and  Fi.  Three  power  func¬ 
tions  were  obtained  as  the  regression  lines  relating 
these  properties: 


(1)  LogeDvs.  Log10  F:  r  =  .94;  b  =  0.83, 
^  =  1.92;*  F  =  0.60D1-92. 

(2)  LoggDj  vs.  LogeD:  r  =  .88;  b  =  0.921; 
Dx  =  0.758  D0-921. 

(3)  Log  F!  vs.  Log  F:  r  =  .97;  b  =  .975; 

Fx  =  .99F-975. 


A  test  of  significance  of  b  in  equation  (1)  above 
shows  that  it  is  not  significantly  different  from 
.857  (.40  >P  >  .30),  and  2.303-b  is  not  significantly 
different  from  2.00,  when  logarithms  are  removed. 
The  relation  F  =  kD2  can  therefore  be  used  as  a 
mathematical  model.  With  increasingly  fine  tex¬ 
tured  topography,  F  increases  proportionally  to  D2, 
showing  that  the  increase  in  total  length  of  drain¬ 
age  lines  is  due  to  an  increase  in  number  of  chan¬ 
nel  segments,  and  not  to  longer  segments  of  each 
order.  An  increase  in  drainage  density  is  then,  ac¬ 
cording  to  this  hypothesis,  a  change  of  scale  of 
topography,  resulting  in  smaller  basins  and  shorter 
channel  segments  of  each  order.  If  a  model  is  con¬ 
sidered  in  which  the  area  of  a  certain  basin  is 
constant,  that  is,  erosion  is  not  allowed  to  enlarge 
the  initial  drainage  area  by  extending  the  major 
divide,  then  the  equation  F  =  kD  2  can  be  trans- 

,  with  the  restriction  that 


formed  into  lu  = 


k£lu 


increases  of  £lu  are  due  solely  to  increases  in  D. 
Then  the  average  channel  lengths  of  all  orders  are 
inversely  proportional  to  the  total  length  of  chan¬ 
nels  at  a  given  instant,  during  the  time  that  the  in¬ 
crease  in  D  is  taking  place. f 

The  exponents  of  D  and  F,  in  equation  (2)  and 
(3)  respectively,  are  not  significantly  different 
from  1.00.  This  suggests  that  the  density  of  first- 
order  channels  increases  linearly  with  the  total 
drainage  density,  and  the  frequency  of  first-order 
channels  increases  linearly  with  the  total  channel 
frequency;  D}  =  k]_D,  and  Fj.  =  k2F.  These  equa¬ 
tions  can  be  transformed  into  li  =  k].£  lu,  and 
ni  =  k2£nu,  where  increases  in  £  lu  and  £1^  are 
restricted  to  increases  in  fineness  of  topographic 
texture.  Combining  these  two  model  equations  gives 


=  2.303  -  loge10. 

fThe  same  may  be  said  of  the  relations  among  several 
basins  of  the  same  planimetric  size,  but  differing  in 
drainage  density. 


ll  =  k3lu.*  According  to  this  hypothesis,  the  aver¬ 
age  length  of  first-order  channels  is  a  constant 
proportion  of  the  average  length  of  channels  of  all 
orders;  the  model  implies  then,  combining  with  the 
first  equation  obtained,  that  average  length  of  first- 
order  channels  is  inversely  proportional  to  the 
total  length  of  channel,  when  the  latter  increases 

—  k  a 

during  an  increase  of  drainage  density,  lj  =  ■  , 

if  area  is  constant.  A  corollary  of  this  hypothesis 
is  that  the  average  maximum  length  which  first- 
order  channels  attain  before  bifurcating  is  less  in 
areas  of  high  drainage  density  than  in  areas  of  low 
density  (Schumm,  1954,  p.  16). 

The  values  which  the  constants  of  proportion¬ 
ality  assume  in  the  model  equations  stated  above 
are  abstract,  fundamental  properties  of  basins, 
and  probably  are  influenced  by  many  of  the  govern¬ 
ing  agents  in  basin  regimen.  These  quantities 
should  be  of  great  interest  for  further  investiga¬ 
tion. 


RELATION  OF  VALLEY-SIDE  SLOPES  TO 
OTHER  BASIN  CHARACTERISTICS 

The  effects  of  seven  variables  on  slope  values 
are  dealt  with  as  simple  two- variable  correlation 
and  regression  relations,  and  also  combined  in 
multiple  correlation  coefficients.  The  variables 
are  (logarithms  taken  where  necessary  for  linear 
plot,  Y  vs.  Xi): 

Xi  wet  soil  strength/1000 
X  2  infiltration  capacity 
X3  Log  Ruggedness  number/order 
X  4  Log  (P-E) 

X  5  Log  relative  relief 
X  g  Roughness  number 
X7  5-year,  1-hour  rain  +  5 
Y  Average  maximum  valley-side  slope. 

The  simple,  two-variable  correlation  coeffi¬ 
cients  obtained  from  analysis  of  both  grouped 
and  ungrouped  data  are  given  in  Appendix  V  and 
below: 

ryi  =  -.423;  ©  =  -.  756  y^.+  28.26;  sy>x  =  7.4 

r y2  =  .746;  ©  =  1.312  f  +  20.7;  sY  x  =  4.9 

r  Y3  =  .777;  ©  =23.71  log  ^  +  30.84;  sy>x  =  4.6 

r  y4  =  .289;  ©  =  9.90  log  (P-E)  +  10.8;  sY .  x  =  7.0 
r  y5  =  .625;  0=  18.69  log  R  +  13.1;  sY>x  =5.7 


♦True  for  third  and  fourth -order  basins;  experimentally 
verified  since  writing. 
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r  Y6  =  -.356;  6  =  -.091M  +  25.46;  s  y  x  =7.6 
(not  significant) 

r Y7  =  _.472;  ©  -4.36q  +  40.96;  sy.x  =  7.2 

r  35  =  -666  r27  =  -.918  r  24  =  .750 

r  23  =  .386  r  12  =  -.0946. 

Multiple  correlation  coefficients  of  three  vari¬ 
ables  were  obtained:* 

R  Y-23  =  .916;  R  y*24  =  .651;  Ry.27  =  *918. 

Three-variable  partial  correlation  coefficients 
were  found  by  combining  the  simple  correlation 
coefficients: 


Y3.5 

=  .620 

r  Y3.2 

=  .797 

r  yi  2  —  —  •  5 3 2 

Y5.3 

=  .230 

r  Y7.2 

=  .804 

r  Y4.2  =  “-614 

Y2.3 

=  .769 

r  Y2.7 

=  .893 

rY2  4  -  .836 

These  coefficients  are  not  part  of  a  complete 
multiple  correlation  study.  They  are  used  only  as 
estimates  of  the  population  correlation  coefficients. 
Thus,  some  of  the  coefficients  calculated  from  lar¬ 
ger  samples  are  judged  significant  though  smaller 
than  ry6  which  is  not  significant. 

The  only  variable  not  judged  significant  is  the 
roughness  number.  However,  the  negative  value  of 
the  regression  coefficient  obtained  is  consistent 
with  a  theory  to  be  developed  below  to  explain  the 
inverse  relation  between  slope  angle  and  wet  soil 
strength;  both  roughness  and  strength  depend  on 
the  soil  texture  and  adequate  sampling  should  de¬ 
tect  a  positive  relation  between  them. 

The  highest  correlations  were  found  between  © 

|| 

and  log  — ,  also  ©  and  f .  These  have  lowest  stand¬ 
ard  error  of  estimate,  that  being  on  the  order  of 
the  usual  standard  deviation  of  measurements  with- 

|| 

in  individual  basins.  Log  —  alone  explains!  60.4% 

of  the  variation  in  © ;  f  alone  explains  55.6%.  The 
third  highest  correlation  was  found  to  be  with  log 
R,  agreeing  with  previous  findings  (Schumm,  1954, 


♦The  computation  of  multiple  correlation  coefficients 
was  restricted  here  because  not  all  intercorrelations 
among  the  X  ^  were  computed.  The  multiple  coefficients 
chosen  for  presentation  show  the  importance  of  the  in¬ 
fluence  of  infiltration  capacity,  log  ruggedness  number/ 
order,  log(P-E)  index,  and  runoff  intensity  on  valley- 
side  slopes.  These  four  variables  are  furthermore  of 
greatest  interest  in  this  study,  summarizing  hydrologic, 
diastrophic,  and  climatic  factors. 

fBy  “explains,”  it  is  meant  that  knowledge  of  X  enables 
us  to  predict  Y  with  error  equal  to  (1  -  r2)  of  previous 
error.  Above,  60.4  =  100r2  ,  etc. 


p.  23,  and  fig.  15).  Log  R  alone  explains  39.1%  of 
the  variation  in  ©.  When  expressed  in  degrees, 
valley-side  slope  varies  proportionally  to  the  loga¬ 
rithm  of  the  relative  relief,  the  logarithm  of  rug¬ 
gedness  number  corrected  for  order,  and  infiltra¬ 
tion  capacity.  The  lowest  correlation  was  found 
with  log  (P-E),  but  it  was  judged  significant  be¬ 
cause  it  was  based  on  a  large  sample. 

Three  multiple  correlation  coefficients  were  cal¬ 
culated  from  the  simple  coefficients  (Croxton,  1953, 
p.  172).  The  values  obtained  were  rather  high,  show¬ 
ing  that  a  large  portion  of  the  variation  in  the  valley- 
side  slopes  is  explained  by  consideration  of  only 
two  independent  variables.  For  instance,  infiltra¬ 
tion  capacity,  and  the  five-year,  one-hour  rainfall 
amount  minus  the  infiltration  capacity  plus  five  (q) 
(which  take  into  account  infiltration  and  the  relation 
of  infiltration  to  average  intensity  of  rain  storms), 
together  explain  84.3%  of  the  variation  in  the  valley- 
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side  slope.*  Infiltration  capacity  and  —  together 

explain  83.8%  of  the  variation;  infiltration  capacity 
and  log  (P-E)  together  explain  72.4%  of  the  varia¬ 
tion  in  ©. 

Nine  of  the  possible  partial  correlation  coeffi¬ 
cients  were  computed,  and  excepting  rys.3  ,  all 
are  significantly  different  from  zero.  This  means 
that  log  R  does  not  improve  the  estimation  of  © 
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when  log  —  is  already  used.  However,  consideration 

of  log  —  when  log  R  is  already  used  does  signifi¬ 
cantly  improve  the  estimation  of  ©.  Consideration 
of  q,  when  f  is  already  taken  into  account,  explains 
about  64.4%  of  the  variation  in  ©  unexplained  by  f 
alone.  In  this  connection,  a  correlation  exists  be¬ 
tween  f  and  log  (P-E),  (r  24  =  .750),  a  relation  well 
known,  and  probably  best  explained  by  the  increased 
amounts  of  humus  and  bacteria  in  soils  withmoister 
environments. 

The  results  of  greatest  interest  here  are  the  re¬ 
lations  between  ©  and  f,  and  ©  and  q.  It  is  apparent 
that  greater  average  amounts  and  intensities  of  run¬ 
off  are  associated  with  lower  slopes.  Conversely, 
areas  with  high  infiltration  capacity  have  low  run¬ 
off,  and  ceteris  paribus,  low  rates  of  slope  erosion; 
valley-side  slopes  in  these  areas  tend  to  become 
and  remain  steep.  In  discussing  a  continuous  down- 
slope  profile,  Gilbert  (1877,  p.  117)  stated  that  the 
declivity  of  slope  is  inversely  proportional  to  the 
amount  of  water  flowing  over  it,  an  observation  that 

♦The  multiple  correlation  coefficient  and  predicted  slope 
values  obtained  from  the  use  of  infiltration  capacity  (f) 
and  5 -year,  1-hour  rainfall  amount  (P)  as  independent 
variables  would  be  identical  with  those  obtained  from  any 
two  independent  linear  combinations  of  f  and  P,  such  as 
f  and  q  =  (P  -f+5),  used  here.  The  reason  for  using  q 
is  that  it  has  physical  meaning,  viz.,  the  intensity  of  run¬ 
off  expected  every  five  years  from  a  one-hour  rain. 
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appears  to  hold  in  general.  Furthermore,  these  data 
show  that  the  relation  between  slope  and  f  or  q  is 
arithmetically  linear: 

0  =  aj  -  k]_  q;  0  =  a2  +  k2  f. 

The  following  working  hypothesis  is  proposed: 
where  infiltration  capacity  is  low,  greater  erosion 
accompanies  greater  runoff,  and  channels  are 
eroded  slowly  in  relation  to  the  adjacent  slopes, 
because  much  debris  from  the  slopes  is  brought 
down  into  them.  As  a  result,  the  slopes  are  lowered 
more  rapidly  than  the  channels,  and  the  angle  of 
slope  declines.  With  higher  infiltration  capacity, 
little  runoff  occurs,  and  the  valley-side  slopes  are 
not  eroded  so  rapidly  as  in  the  previous  case.  Con¬ 
centration  of  the  runoff  in  channels  greatly  in¬ 
creases  its  eroding  efficiency.  An  increment  of 
water  to  channel  flow  comes  by  interflow  (subsur¬ 
face),  and  this  addition  serves  to  increase  erosion 
in  the  channels,  but  the  water  thus  added  did  not 
have  the  opportunity  to  erode  the  valley  sides.  Thus 
when  the  channel  below  a  segment  of  slope  is 
deepened,  a  proportionate  amount  of  erosion  has 
not  occurred  on  the  slope  above.  The  tendency  is 
thus  to  remove  the  support  of  the  surface  of  the 
valley  side  afforded  by  the  toe,  thereby  allowing 
an  acceleration  of  creep,  and  consequent  steepen¬ 
ing.  The  quantity  of  interflow  is  affected  by  the 
structure  of  the  soil,  and  is  favored  by  shallow 
soils,  with  presence  of  parent  rock  or  other  im¬ 
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Figure  8.  Hypothetical  channels  (A)  with  low  infiltration 
capacity,  high  runoff  and  erosion  rates,  low  interflow, 
debris -laden  channel  and  low  valley-side  slopes;  (B)  with 
high  infiltration  capacity,  low  runoff  and  erosion  rates, 
high  interflow,  eroded  channel,  and  high  valley-side 
slopes. 


permeable  layer  a  short  distance  below  the  sur¬ 
face  of  the  ground.  This  would  be  the  case  in  much 
of  the  semiarid  and  arid  west  (Linsley,  Kohler, 
and  Paulhus,  1949,  p.  388). 

Climatic  conditions  enter  further,  because  as 
the  average  rainfall  intensity  increases,  the  aver¬ 
age  amount  and  intensity  of  runoff  increases.  Fur¬ 
thermore,  greater  P-E  index  tends  to  increase  in¬ 
filtration  capacity,  thus  producing  steeper  slopes. 

The  effect  of  soil  strength  on  the  average  valley - 
side  slope  maintained  in  an  area  of  mature  topog¬ 
raphy  could  be  deduced  from  several  reasonable 
premises,  but  the  results  would  conflict  unless  the 
quantitative  importance  of  all  pertinent  variables 
is  known.  A  theory  will  be  constructed  that  is  con¬ 
sistent  with  the  results  obtained  by  statistical 
analysis  of  the  measurements  made  in  this  study, 
but  confirmation  of  the  theory  depends  on  further 
hydrologic  research.  We  know  (1)  that  ©  varies 
inversely  with  Sw  (wet  soil  strength),  because 
r y i  =  -.423  is  significantly  less  than  zero;  (2)  that 
the  correlation  between  ©  and  Sw  is  improved 
when  the  effect  of  infiltration  capacity  is  removed 
(r y i .  2  =  -  - 532);  (3)  that  wet  soil  strength  is  greater 
for  coarser  soils  (confirmed  by  a  one-way  analysis- 
of-variance  test,  not  reproduced  here).  Excluding 
cliffs  from  consideration,  the  effect  of  soil  strength 
must  be  through  erosive  processes  rather  than 
stability,  because  weaker  soils  tend  to  be  stable  at 
lower  slope  angles,  in  opposition  to  (1)  above. 

Analysis  of  covariance  of  the  data  in  Table  4 
gives  the  following  table: 


d.f. 

Ex2 

£xy 

CM 

Ml 

d.f. 

sy’i2 

M.  S. 

Among 

means 

Within 

2 

3.25 

40.82 

714.0 

2 

521.2 

260.6 

groups 

20 

12.59 

18.02 

686.4 

19 

660.6 

34.7 

Total 

22 

15.84 

58.84 

1400.4 

21 

1181.8 

where  X  =  f ,  Y  =  0 .  The  above  was  obtained  by 
grouping  0  and  f  into  three  classes  depending  on 
soil  texture.  Fine-grained  soils  give  high  values  of 
0,  and  coarse-grained  soils  give  low  values  of  0. 
The  F- ratio  for  slope  means  unadjusted  for  infil¬ 
tration  capacity  is  10.4  (significant  at  .001),  and 
the  F-ratio  for  slope  means  after  adjustment  for 
infiltration  capacity  is  7.52  (significant  at  .01). 
Therefore,  soil  texture  influences  valley- side 
slope  through  some  agency  other  than  infiltration 
capacity  alone.  The  correlation  between  0  and  Sw 
is  probably  due  to  this  agency. 

In  the  following  proposal  of  a  theory  to  explain 
the  inverse  relation  between  0  and  Sw,  it  is  as¬ 
sumed  (1)  that  the  hypothetical  areas  under  consid¬ 
eration  have  identical  climates,  (2)  that  the  bed 
rock  near  the  surface  is  essentially  homogeneous 
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and  does  not  differ  greatly  in  mechanical  compo¬ 
sition  from  the  top  soil  layer,  and  (3)  that  the  area 
does  not  lie  down-stream  from  a  source  of  debris 
of  contrasting  caliber  to  that  produced  within  the 
area. 

In  an  area  of  low  Sw,  silt  and  clay  predominate 
in  the  soil,  and  gravel  is  absent.  Debris  brought  in¬ 
to  the  channels  by  runoff  from  the  slopes  remains 
in  suspension  and  is  immediately  transported  out 
of  the  area.  Consequently,  the  channels  remain  free 
from  bed  load  and  are  subject  to  immediate  erosion 
by  the  concentrated  stream  of  water.  While  runoff 
lasts,  the  channels  are  continually  deepened  and/or 
widened,  and  are  not  refilled  as  discharge  wanes— 
silt  and  especially  clay  remaining  in  suspension 
even  in  low  velocities.  Because  of  the  greater  ef¬ 
ficiency  for  erosion  of  the  concentrated  stream 
over  slope  wash,  the  channels  are  lowered  faster 
than  the  slopes,  the  toes  of  the  slopes  are  removed, 
and  the  adjacent  slopes  are  steepened  until  a  steady 
state  is  achieved  through  a  reduction  of  the  rate  of 
erosion  in  the  channels  relative  to  the  rate  of 
erosion  of  the  slopes.  The  self-limiting  mechan¬ 
ism  is  the  result  of  increased  erosion  and  creep 
on  the  steepened  slopes,  causing  channel  gradient 
to  be  lowered  to  that  value  necessary  to  transport 
the  small  quantity  of  fine-caliber  bed  load.  In  such 
areas,  at  steady  state  the  slope  orthogonals  will 
intersect  the  channel  at  high  angles,  and  the  inflec¬ 
tion  angles  will  be  acute. 

In  an  area  of  high  Sw  gravel  predominates  in  the 
soil.  Slope  runoff  sweeps  debris  into  the  subjacent 
channel  where  coarse  material  accumulates  as  bed 
load.  Assuming  a  case  in  which  valley-side  slopes  are 
initially  steeper  than  equilibrium  requires  (e.g. 
where  channels  are  trenched  or  rejuvenated),  the 
greater  declivity  of  the  slopes  permits  movement 
of  gravel  and  cobbles  by  rolling  and  sliding  under 
relatively  low  down-slope  discharges,  but  in  the 
channel,  higher  velocities  and  discharges  are  re¬ 
quired  for  transportation  because  the  down-stream 
gravitational  component  is  low.  Down- cutting  of  the 
bed-rock  channel  cannot  occur  until  the  bed  load 
is  removed,  and  except  where  the  channel  is  sup¬ 
plied  by  ground  water,  the  bed  load  is  replenished 
by  material  brought  down  from  the  slopes  when¬ 
ever  sufficient  discharge  occurs  to  transport  all 
of  the  bed  load  already  present.  Thus,  channel  down¬ 
cutting  is  a  relatively  rare  event,  but  the  slopes  can 
be  eroded  frequently.  As  channel  discharge  wanes, 
after  a  runoff-producing  rain,  any  depressions  that 
had  formed  in  the  channel  are  refilled  by  removal 
of  debris  from  suspension.  Erosion  of  the  slope 
toes  by  lateral  corrasion  can  occur  at  the  outside 
of  sharp  channel  bends,  with  consequent  local  steep¬ 
ening  of  the  cut  bank  slope,  but  a  compensating 
lowering  of  opposite  slopes  results  in  no  quantita¬ 
tively  important  change  over  the  basin  as  a  whole 


(Strahler,  1950,  pp.  812,  813).  The  slopes  are  there¬ 
fore  lowered  at  a  greater  rate  than  the  channels, 
and  reach  a  steady  state  at  some  relatively  low 
angle,  when  the  down- slope  component  of  gravity  is 
so  adjusted  as  to  allow  only  that  amount  of  coarse  de¬ 
bris  to  enter  the  channel  which  can  be  transported 
by  stream  flow.  The  channel  gradient  necessary 
to  transport  the  large  amount  of  coarse  bed  load  is 
relatively  steep;  thus  the  angle  between  slope  or¬ 
thogonals  and  channel  will  be  very  acute,  and  the 
inflection  angles  obtuse. 

In  extreme  cases,  if  the  amount  of  debris  sup¬ 
plied  to  the  channels  is  quite  great,  it  seems  prob¬ 
able  that  the  channel  will  aggrade  rapidly  and  as¬ 
sume  a  braiding  habit.  The  slope  toes  along  the 
entire  course  of  the  channel  may  then  be  cut  by 
lateral  corrasion  and  the  slopes  be  maintained  at  a 
steep  angle  until  the  maximum  breadth  possible  for 
the  drainage  area  has  been  attained.  Basins  with 
this  environment  were  not  included  in  the  areas 
studied. 

The  mechanism  proposed  above  for  determining 
the  relative  rates  of  slope  and  channel  erosion,  i.e. 
the  abundance  of  bed  load,  is  the  same  as  that  pro¬ 
posed  to  explain  the  relation  of  slope  angles  to  in¬ 
filtration  rates.  However,  areas  most  favorable  for 
the  accumulation  of  bed  load  because  of  low  infiltra¬ 
tion  capacity-high  runoff  intensity  are  usually  least 
favorable  for  the  production  of  coarse  debris  be¬ 
cause  of  the  prevalence  of  silt  and  clay.  Therefore, 
the  effect  of  infiltration  capacity  is  generally  op¬ 
posite  the  effect  of  debris  availability,  though  both 
depend  on  soil  texture  and  act  through  the  same 
mechanism.  This  explains  the  indicated  increased 
correlation  between  0  and  Sw  when  the  effect  of  in¬ 
filtration  capacity  is  removed.  Infiltration  capacity 
depends  also  on  climate,  so  the  net  effect  of  soil 
texture  cannot  be  judged  qualitatively  except  per¬ 
haps  in  extreme  cases.  The  equilibrium  valley- 
side  slopes  for  given  conditions  of  climate,  infiltra¬ 
tion  capacity,  soil  texture,  etc.  are  best  estimated 
quantitatively  by  use  of  regression  formulas. 

Much  further  investigation  is  necessary  before 
definite  hypotheses  can  be  formulated  to  include 
more  than  a  few  of  the  variables  at  one  time  for  the 
entire  environmental  conditions  of  soil,  climate, 
history,  structure,  and  slope. 


TEXTURE  OF  TOPOGRAPHY 
AS  AFFECTED  BY  CLIMATIC 
AND  MANTLE  CHARACTERISTICS 

Information  from  23  basins  subjected  to  full  field 
investigation  was  arranged  for  multiple  regression 
and  correlation  analysis.  Drainage  density  (log  D) 
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was  the  dependent  variable, 
dependent  variables:* 

with 

the  following  in¬ 

Xx: 

log  (P-E) 

X4: 

log  roughness 

X2: 

log  infiltration 

number 

X3: 

percent  bare 
area/100 

X5: 

log  q. 

The  estimating  equation  found  by  simultaneous 
solution  of  the  five  normal  equations  is  given  in 
Appendix  VII  and  below: 

Y  =  4.82  -  1.40  Xx  -  0.859  X2  +  0.0904 
X3  -  0.074  X4  -  1.59  X5. 

The  standard  error  of  estimate  is  0.1224.  The 
simple  correlation  coefficients  are: 


rYl  = 

-.943 

r12  - 

.808 

r24  = 

-.150 

II 

CM 

f-l 

-.793 

II 

CO 

r—H 

u 

-.903 

r  25  = 

-.906 

rY3  = 

.900 

r  14  = 

-.318 

r34  = 

.241 

II 

>1 

u 

.161 

r  15  = 

-.748 

r  35  = 

.689 

r  Y5  = 

.658 

r  23  = 

-.848 

r  45  = 

.266 

The  multiple  correlation  coefficient  is  Ry.  12345  = 
0.965,  and  thus  93.2%  of  the  variation  in  D  is  ex¬ 
plained  by  including  only  climatic  and  surface 
properties.  If  the  effects  of  certain  morphologic 
features,  such  as  R  or  ©  were  considered  in  addi¬ 
tion,  possibly  a  greater  amount  of  variation  would 
be  explained.  An  analysis  of  variance  test  for  the 
significance  of  Ry  12345  (Anderson  and  Bancroft, 
1952,  p.  163,  172)  'is: 


Source  of 
Variation 

d.f. 

S.  S. 

M.  S. 

F 

Regression 

5 

4.489 

0.898 

46.29 

Error 

17 

.329 

0.0194 

Total 

22 

4.818 

0.219 

ry.  12345  is  significant  at  .001. 


♦These  particular  independent  variables  are  used  be¬ 
cause  they  summarize  the  climatic,  hydrologic,  and  sur¬ 
face  properties  that  determine  the  topographic  texture 
of  a  region.  From  two-variable  plots,  it  was  found  that 
log  D  is  a  linear  function  of  the  logarithms  of  f,  M,  and 
q,  so  the  log  transformation  is  used  for  these.  Log  D  is 
linear  with  b/100,  and  therefore  the  logarithmic  trans¬ 
formation  was  not  used  for  it. 


Multiple  correlation  coefficients  involving  only 
log  (P-E),  log  f,  and  log  q  are: 

ry.15  =  -946 
R Y. 125  =  *960 

Selected  partial  correlation  coefficients  are:* 
a)  Three  variables: 


rY1.2  " 

-.842 

r 13 . 2  “ 

-.698 

rY1.3  = 

-.969 

r15.2  = 

-,064(N.S.) 

rY1.5  = 

-.902 

r35.2  = 

-.353(N.S.) 

rY2.1  = 

-  .158(N.S.)t 

r  24. 1  = 

.191(N.S.) 

rY3.1  = 

.341(N.S.) 

r Y3.2  = 

.704 

rY4.1  = 

-.440 

r  Y5. 1  = 

-  .214(N.S.) 

r  Y5.2  = 

-  .234(N.S.) 

Four  variables: 

r  Y1.23  = 

-.689 

r15 . 23  = 

-.463 

r  Y1.25  = 

-.883 

r  Y2. 15  = 

-.518 

r  Y3.12  = 

,306(N.S.) 

r  Y4. 12  = 

-.424(.10) 

rY5.12  = 

-.533 

r  Y5.23  = 

.0223(N.S.) 

c)  Five  variables: 
r  Yl.  325  =  “-766 
r  Y3. 125  =  .0139(N.S.) 

The  following  inferences  may  be  made  from  the 
quantities  tabulated  and  presented  above: 

1.  Log  D  =  Y  shows  significant  correlations 
with  (1)  log  (P-E)  index,  (2)  log  percent  bare  area/ 
100,  (3)  log  infiltration  capacity,  and  (4)  log  runoff 
intensity,  in  decreasing  order  of  importance  as  de¬ 
termined  from  the  simple  correlation  coefficients. 
Log  roughness  number  alone  appears  to  have  no 
correlation  with  Y.  These  variables,  together,  ex¬ 
plain  all  but  6.8%  of  the  variation  in  Y,  and  we  may 
state  that  the  climatic,  hydrologic,  and  surface 
properties  of  an  area  are  of  overwhelming  impor- 

*The  partial  correlation  coefficients  presented  were 
chosen  mainly  to  clarify  the  relative  importance  of  the 
P-E  index,  infiltration  capacity,  and  runoff  intensity  in 
determining  drainage  density. 

t(N.S.)  denotes  not  significant  at  .10;  (.10)  denotes  sig¬ 
nificant  at  .10  but  not  .05;  all  others  significant  at  .05  or 
less. 
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tance  in  determining  topographic  texture.  Further¬ 
more,  their  influences  can  be  detected  by  relatively 
simple  methods  of  field  measurements. 

2.  The  value  of  Ry.i2  5  indicates  that  of  the  93.2% 
of  the  variation  in  Y  explained  by  consideration  of 
all  five  independent  variables,  92%  is  explained  by 
Xj_,  X2,  and  X5;  percent  bare  area  and  surface 
roughness  are  very  minor  influences,  when  the 
others  are  already  considered  (note  also,  rY3  195  = 
.0139). 

3.  The  regression  coefficients,  which  are  par¬ 
tial  derivatives  of  Y  with  respect  to  each  independ¬ 
ent  variable,  with  the  others  held  constant,  gives 
the  following  decreasing  order  of  importance  to 
the  independent  variables:  (1)  log  runoff  intensity, 

(2)  log  (P-E)  index,  (3)  log  infiltration  capacity, 

(4)  log  percent  bare  area/100,  and  (5)  log  rough¬ 
ness  number.  The  discrepancy  in  order  of  impor¬ 
tance  results  in  part  from  intercorrelations  among 
the  independent  variables,  and  in  part  from  the 
scale  of  the  measures. 

4.  The  partial  correlation  coefficients  indicate 
that  when  X4  is  already  considered,  the  only  vari¬ 
able  that  is  significant  is  X4  (roughness).  If  this 
is  not  a  spurious  result,  then  when  the  P-E  index 
is  held  constant,  drainage  density  is  inversely  re¬ 
lated  to  roughness.  (Or,  areas  of  greater  rough¬ 
ness  are  associated  with  lower  drainage  density.) 

It  should  be  noted  that  X4  shows  the  lowest  simple 
correlation  with  X4  of  all  the  independent  variables. 
Conversely,  when  the  effects  of  X2,  X3,  X5,  X2  and 
X3,  X2  and  X5,  and  X2,  X3,  and  X5  are  removed, 

Xi  still  has  a  significant  correlation  with  drainage 
density.  The  conclusion  drawn  is  that  the  P-E  index 
influences  drainage  density  through  the  agencies 

of  infiltration  capacity,  amount  of  cover,  and  run¬ 
off  intensity,  but  probably  not  through  surface 
roughness.  In  addition,  the  P-E  index  influences 
drainage  density  by  other  mechanisms,  not  ex¬ 
plained  here. 

5.  The  infiltration  capacity  of  an  area  is  deter¬ 
mined  by  the  P-E  index  of  the  region,  the  surface 
cover,  and  by  the  type  of  rock  from  which  the  soil 
was  developed.  The  amount  of  surface  cover  is  de¬ 
termined  by  the  P-E  index  also.  The  runoff  inten¬ 
sity  depends  on  average  rainfall  intensity  (a  func¬ 
tion  of  the  same  factors  that  determine  P-E  index), 
infiltration  capacity,  and  hence  on  surface  cover 
and  the  P-E  index.  In  equilibrium  conditions,  the 
drainage  density  will  attain  and  stay  at  a  value  de¬ 
pendent  on  all  these  variables,  and  others  not  con¬ 
sidered  here.  However,  if  say,  the  infiltration 
capacity  and  percent  bare  area  were  changed  by 
fire,  grazing,  cultivation,  etc.,  to  values  not  in 
agreement  with  the  P-E  index,  a  change  in  drainage 
density  would  result,  and  the  new  value  would  no 
longer  agree  with  that  expected  from  consideration 
of  the  P-E  index  alone. 


GENERALIZATIONS  ON  INVESTIGATIVE 
PROCEDURES 


The  following  are  generalizations  which  seem 
appropriate  concerning  the  philosophy  of  the  inves¬ 
tigation  as  a  whole: 

1.  It  has  been  shown  that  when  various  proper¬ 
ties  of  drainage  basins  are  analyzed,  occurrences 
are  found  that  would  not  reasonably  be  expected  to 
result  from  operation  of  chance,  or  sampling  error 
alone. 

2.  It  was  proposed  that  these  non-chance  occur¬ 
rences  represent  in  some  manner  true  laws  of  na¬ 
ture,  hence  mathematical  models  were  suggested 
consistent  with  what  is  known  about  these  laws 
from  sampling  and  statistical  analysis. 

3.  It  is  the  nature  of  statistical  laws,  such  as 
obtained,  that  a  single  observation  contrary  to  the 
supposed  law  will  not  refute  it,  as  the  element  of 
random  error  is  always  present  in  any  observation 
or  measurement,  and  is  here  presumed  to  exist  in 
nature  when  a  large  system  is  considered. 

4.  Because  it  is  impossible  to  consider  all  as¬ 
pects  of  the  environment  of  a  drainage  system,  or 
any  physical  system  whatsoever,  we  will  never 
know  that  all  factors  influencing  observations  are 
included  in  the  analysis.  It  seems  most  unlikely 
that  enough  factors  could  ever  be  considered  at 
once  to  predict,  completely  without  error,  the 
lengths,  gradients,  and  positions  of  all  channels  in 
a  basin;  adherence  to  the  notion  that  perfect  pre¬ 
diction  could  be  obtained  by  finite  extension  of 
purely  analytical  differential  equations  is  then 
rather  useless.  A  better  point  of  view  is  to  acknowl¬ 
edge  that  some  natural  agencies  operate  in  a  man¬ 
ner  that  cannot  be  predicted  except  by  the  use  of 
probability  theory,  and  then  build  a  theory  based  on 
statistical  mechanics  rather  than  on  Cartesian- 
Newtonian  mechanics. 

5.  Some  of  the  non- random  occurrences,  the 
correlations  obtained  above,  may  reflect  the  agency 
of  cause  and  effect.  However,  because  of  the  com¬ 
plexity  of  the  interconnections  of  causative  factors 
and  the  presence  of  random  error,  simple  cause 
and  effect— meaning  relations  of  forces,  motions, 
and  bodies— seems  inappropriate,  and  instead  a 
complex  feedback  among  cause-like  agencies  and 
effect-like  agencies  seems  more  pertinent  (in  this 
case,  properties  of  climate,  lithology,  history, 
vegetation,  mantle,  and  form),  when  the  totality  of 
environment  is  included. 

6.  This  does  not  mean  that  physical  laws  of 
cause  and  effect  do  not  operate,  or  that  methods  of 
analysis  by  differential  equations  cannot  be  applied. 
Simple  mechanical  models  of  a  natural  system  are 
satisfactory  only  in  those  very  small  parts  limited 
both  in  time  and  physical  boundaries. 


Summary  of  Conclusions 


The  topographic  texture  and  average  valley- 
side  slopes  of  fluvially- controlled  landscapes  are 
related  to  the  causative  and  determinative  factors 
of  climate  and  lithology.  However,  it  is  difficult, 
or  even  meaningless,  to  consider  isolated  elements 
of  the  environment  of  a  drainage  basin  in  connec¬ 
tion  with  these  geometric  properties  of  the  basin 
because  the  geologic  and  climatic  elements  act  in 
complex  ways,  through  the  agencies  of  vegetal 
growth,  soil  formation,  runoff  and  erosion,  infil¬ 
tration,  and  soil  creep.  Thus  in  this  study,  as  many 
elements  as  practicable  were  analyzed  concur¬ 
rently. 

Drainage  density  is  a  negative  power  function 
of  the  P-E  index,  a  measure  of  the  availability  of 
moisture  to  vegetation,  but  also  depends  on  the 
percentage  of  bare  area,  infiltration  capacity,  soil 
strength,  and  intensity  of  runoff.  These  in  turn  de¬ 
pend  on  the  P-E  index  and  are  the  agents  by  which 
climate  influences  drainage  density.  Lithology  in¬ 
fluences  drainage  density  by  determining  the 
amount  of  clay  in  the  soil  (though  this  depends 
also  on  the  climate),  and  so  affects  the  infiltration 
capacity  and  soil  strength.  Perhaps,  also,  lithol¬ 
ogy  indirectly  affects  climate,  because  in  the  areas 


studied  moist  climates  are  confined  to  mountain¬ 
ous  regions,  and  thus  are  found  primarily  on  re¬ 
sistant  rocks.  Lastly,  recency  of  uplift  influences 
drainage  density  and  valley-side  slopes  because 
areas  of  great  immediate  relief  have  large  relief 
ratios  and  generally  coarse  texture.  To  develop 
fine  texture  in  such  an  area  would  require  impos¬ 
sibly  steep  slopes  for  equilibrium  conditions. 

The  valley-side  slopes  of  drainage  basins  vary 
with  lithology,  being  steepest  on  shale  and  perhaps 
limestone,  and  gentlest  on  schist  and  acid  volcanic 
lithologies.  Steep  slopes  are  associated  with  large 
relief  ratio,  high  infiltration  rate,  low  wet  soil 
strength,  low  runoff  intensity,  and  high  P-E  index. 
This  suggests  that  steep  slopes  are  produced  and 
maintained  in  areas  favoring  channel  deepening 
over  slope  erosion.  In  such  areas,  water  passing 
into  the  soil  by  infiltration  has  no  opportunity  to 
erode  the  slope  and  contribute  detritus  to  the  chan¬ 
nel  load,  but  emerges  directly  into  the  channel  and 
is  immediately  available  for  channel  erosion.  Thus, 
basins  with  steep  slopes  should  be  expected,  on  the 
average,  to  have  low  drainage  density.  This  rela¬ 
tion  was  found  in  the  correlation  analysis  of  the 
two  variables. 
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Appendix  I 


DERIVATION  OF  RELATIONS  AMONG  GEOME¬ 
TRICAL  ELEMENTS  OF  VALLEYS 

Valley- side  slopes,  inflection  angle,  and  channel 
gradient  are  related  by  laws  governing  lines  and 
planes,  assuming  that  the  valley  sides  are  planar, 
and  therefore  meet  in  a  straight  line  that  is  the 
channel.  Consider  the  contour  map  (a)  showing  a 
portion  of  a  basin  with  valley  sides  7r1  and  772  ,  with 
unequal  slopes  O^and  02  ,  intersecting  in  a  line 
(channel)  A3. 

The  projection 
of  the  channel, 

A3,  on  the  map 
plane  is  taken 
as  the  Y  axis, 
the  Z  axis  is 
vertical,  and 
the  X  axis  is 
perpendicular 
to  the  YZ  plane, 
as  shown  in  (a). 

If  r^  and  r2  are 
drawn  in  the 
XY  plane  from 
a  point  0  on  the  Y  axis  perpendicular  to  a  contour 
line,  say  200,  on  each  side  of  the  valley,  two  right 
triangles  are  formed,  ABO  and  ACO,  and  the  fol¬ 
lowing  relations  are  true: 

ei  =  90  - 

(1)  e2  =  90  -  £2 


©!  ,  the  slope  of  77^,  is  the  angle  between  pi  and 
the  Z  axis,  ^  is  the  projection  of  pi  on  the  XY 
plane,  and  (q  is  the  angle  between  r^  and  the  Y 
axis.  If  a ±  and  b1  are  the  angles  between  p l  and 
the  X  and  Y  axes  respectively,  then: 

cos  a  1  =  Aj  =  sin  sin 
(2)  cos  bx  =  Pi  =  sin  ©!  cos  |i, 
cos  ©1  =  Ui  =  cos  ©!• 

Xy,  pi,  and  iq  ,  are 
the  direction  co¬ 
sines  of  77  1 . 

If  similarly  th 
perpendicular,  p2, 
to  77  2  is  drawn,  as 
shown  in  figure  (c), 
the  projection  of 
p2  on  the  XY  plane 
is  r2,  angle  a2  is 
between  p2and  the 
negative  end  of  the 
X  axis,  angle  b2  is 
between  p2  and 
the  Y  axis,  and  angle  ©2  is  between  p2  and  the  Z 
axis.  By  geometry,  the  direction  cosines  of  772 


are: 

11 

CM 

-  sin 

02  sin 

(3) 

V-2  = 

sin 

©2  cos 

v2  = 

cos 

©2 . 

^  =  ei  +  e2  =  180  -  (£1  +  &), 

where  e-^  ,  e2  are  the  angles  the  contour  line  makes 
with  the  projected  channel  (Y  axis),  and  ^  ,  |2  are 
the  acute  angles  which  r1  and  r2  make  with  the  Y 
axis  (axil  angles),  and  is  the  inflection  angle  of 
the  contour  line. 

If  a  perpendicular  to  77 1 ,  p9  be  drawn,  as  in  fig¬ 
ure  (b), 


y 


By  (1)  it  is  possible  to  eliminate  £  \  and  £2  from 
the  direction  cosines,  and  obtain  expressions  in 
e  1  and  e2  ,  O^nd  02  .  For  771: 


Al  = 

sin 

01 

cos 

(4) 

M  !  = 

sin 

©1 

sin 

v\  = 

cos 

©li 

1 
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and  for  77  2  : 

(5) 


A 


2 
M2 
v  2 


-sin  ©2  COS  62  , 
sin  ©2  sin  e  2  , 
cos  ©2  . 


A 3,  the  intersection  of  tt  1  and  772  must  lie  in  the 
YZ  plane,  because  by  definition  the  Y  axis  is  the 
projection  of  A 3.  The  direction  cosines  of  the  YZ 
plane  are: 


(6) 


Ayz  ~  1 
Pyz  =  0 
v  YZ  =  0. 


Since  A3  is  perpendicular  to  both  pi  and  the  X  axis, 
then  the  direction  numbers  of  A3  can  be  found  by 
evaluating  the  determinants: 


(7) 


(8)  m3,= 


(9)  n3  = 


*3  ~ 


Pi  VI 

Ml 

PYZ  UYZ 

II 

0 

0 

=  0, 


v1  Ax 

cos  ©3  sin  ©i  cos  ei 

V\Z  A YZ 

0  1 

*1 

Mi 

^  YZ 

Myz 

=  COS  ©x , 


sin  ©jcos  e1  sin  ©jsinei 


0 


=  -sin  ©  jsin  ei 
(Olmsted,  1947,  p.  23). 

The  direction  cosines  of  A3are: 

(10)  *3=  0, 

COS  ©1 

(11)  m3  = 


(12)  ^  3  = 


+  (cos2  ©1  +  sin2©i  sin2  e  1  )  ^ 
-  sin  ©j^sinej 

+  (cos  2  ©  1  +  sin  2  ©1  sin2  ei  ) 


X 


Since  p3  is  the  cosine  of  the  angle  between  A3  and 
the  Y  axis,  cos*1  p3  is  the  gradient  of  the  channel, 

y. 


(13) 


cos  © 


cos  y  -  n-  = 


1 


±  (cos2©1+sin  2©j^  sin  2ei  )  ^ 
Squaring  and  simplifying, 


cos  © 


cos  2y  = 


cos  2©1  +sin2©1  sin2 e  1 


sin 2e  =  cot2  ©2 


sin2y 
cos  2y 


(14)  sin  e 1  =  ± 


+  tan  y 


tan© 


1 


The  direction  numbers  of  A3  can  also  be  ex¬ 
pressed  in  terms  of  the  direction  cosines  of  772 
and  the  YZ  plane,  and  by  analogy: 


(15) 


sin  e  ©  =  - 


tany 
tan  ©2 


Assuming  that  ei  and  e2  will  always  be  acute, 
which  is  invariably  true  of  natural  channels,  only 
the  positive  signs  can  be  used  in  (14)  and  (15).  From 
(14): 

<16>  =  tanlp&rip 


and  from  (15) 


(17) 


tan  e0  = 


tany 


2  tan  ©2  cos e  2  ’ 


Adding  (16)  and  (17),  and  solving  for  tan  y, 

tan  y  tan  y 


tan  e  1  +  tan  e  0  =  *r - ^ -  -r  t - ^ - 

1  2  tan  ©icosei  tan  ©2  cos  62 


tany  = 


(18) 


(tane  1 +tane 2 )  tan  ©itan  ©2cose  1  cose2 
tan  ©2  cos e 2  +  tan  ©1  cose  1 

tanei  +  tane  2 


tany  = 


cot  ©x  cot  ©2 


cose 


1 


cose , 


Equation  (18)  is  the  general  expression  of  the 
relations  among  channel  gradient,  valley-side 
slopes,  and  the  components  of  the  inflection  angle 
ej  and  e  2  ,  where  e  1  +  e  2  =  ^  • 

In  the  special  case  where  ©1  =  ©2,  that  is,  the 
two  valley  sides  have  the  same  slope,  ei  =  e2  =lA 
&  ,  and  (18)  can  be  simplified  to: 


(19) 


tan  y  =  sin  tan  0. 
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If  £  is  the  dihedral  angle  between  the  planes 
7r ^  and  77 2  ,  the  relation  among  £  ,  €, ,  e2  ,  ©i  ,  and 
©2  is  easily  obtained.  By  a  theorem  in  solid  analy¬ 
tic  geometry, 

(20)  cos  ^  =  X^  X2  "i"  M  2.  1^2  ^1  ^2  * 

where  X,  p,  and  y  are  the  direction  cosines  of  the 
planes.  From  (4)  and  (5), 

cos  £  =  -  sin  ©2  sin  ©2cos  £2  cos €2  +sin 

©2  sin  ©2sine2Sine2  +  cos  ©2C0S  ©2. 

By  collecting  terms  and  converting  to  functions  of 


£2  +  e 2  >  and  ©2  +  ©2>  the  following  expression  is 
obtained: 

o 

COS  £  =  COS (©2+  ©2) 'COS  V2(£2+e2)+COS(©2-©2)’ 

sir2  i/2(e1+e2); 
or 

(21)  cos£  =  cos  (©2 +©2  )-cos2  V2i^+cos(©2  -  ©2)* 

sin2 

In  the  case  that  ©1  =  ©2,  (21)  can  be  simplified  to 

(22)  cos£  =  cos2©  -  cosi^sin2©. 


L 
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Appendix  II 


INDEX  OF  STUDY  AREAS 

The  quadrangles  listed  below  correspond  to  numbered  squares  within  each  state,  as 
given  on  the  index  maps. 

Arizona 

1.  Chinle  Badlands  (Special  Map) 

2.  Prescott 

3.  Saguaro  National  Monument  (Special  Map) 

4.  Sonoita 

5.  Mt.  Hughes 

6.  Harshaw 

Colorado 

1.  Nederland 

2.  Black  Hawk 

3.  Morrison 

4.  Starkville 

New  Mexico 

1.  Jarosa 

2.  Cerro  del  Grant 

3.  Bland 

4.  Santa  Fe 

5.  Paxton  Springs 

6.  Sedillo 

7.  Twin  Sisters 

8.  Allie  Canyon 

9.  Hurley  East 

10.  Big  Burro  Mountains 

Utah 

1.  Lehi 

2.  Springville 

3.  Spanish  Fork  Peak 
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INDEX  MAP  OF  STUDY  AREAS 


Appendix  III 


DERIVATION  OF  RELATION  BETWEEN  DIAME¬ 
TER  OF  IMPRINT  AND  VOLUME  OF  IMPRINT, 
AS  USED  IN  SOIL  STRENGTH  TEST 

The  imprint  of  the  shot  on  a  planar  surface  is  a 
spherical  segment  with  depth  h,  diameter  d,  and 
volume  V.  The  radius  of  the  shot,  r,  is  0.188  feet. 
An  expression  for  V  in  r  and  d  must  be  found. 


(4)  V  =  -  7r  (  - — dy 

yj  (r2-y2)* 

is  the  volume  of  a  segment  below  a  line  at  X1;  with 
a  radius  of  base  +  yx  .  The  integration  of  (4)  is 
easily  done:  let  u  =  r2  -  y2;  then  y2  =  r2  -  u, 
and  du  =  -2  y  dy. 


X7 


(a) 


Figure  (a)  is  a  semi-circle  with  a  center  at  0,  and 
radius  r.  Since  only  the  portion  in  the  positive  side 
of  the  Y  axis  is  considered,  its  equation  is 

(1)  X  =  +  (r2-y2)^. 

For  a  circle, 


I  =  -j-j  y  2(r2-y  2)  (-2  y  dy)  =  y|  (r  2  -u)u'^du. 

By  usual  methods: 

(5)  I  =  77(r2u'^  -  1/3  u'4 )  +  C, 

and  converting  back  to  y’s 


(6)  I  =  77 


r  2(r  2 -y  2 )%  -  1/3  (r  2 -y  2)34 


Then  from  (4), 


V  =  7 7 


r2(r2 -y2  )1/2  -  1/3  (r  2 -y'z) 


2  „  2\3/2 


=  7 t  J  2/3r3 


+  C. 


0 


yi 


r2(r2-y2)K2  _l/3(r  2-y  2)34 


or 


(7)  V  =  2/3  77  r3  -  77  r2  (r2  -y\)Vl  +  j  (r2-y  2)3/' . 

Since  the  relation  between  diameter  and  volume  of 
the  imprint  is  the  main  concern,  let 


(2) 


ex  =  .  ydy  =  .  ..ydy 

X  (r2-y2)^ 


and  the  element  of  volume,  dV,  is 


d  =  2  y, 


y 


d 

2’ 


(3) 


dV  =  77y2  dX  =  - 


77y3 

(r2  -y2)Vl 


since  a  sphere  is  being  considered,  which  is  ob¬ 
tained  by  rotating  the  circle  about  the  X  axis.  By 
integration, 


then 

(8)  V  =  2/3?7r3-  77r2(r2-  |  j1 *  +  |(r2“  |  )^* 

For  the  particular  shot  used,  r  =  0.188  feet,  and 
(8)  becomes 
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- 
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(9) 


V 


=  0.0138  -  0.1107  (.0352 

(  d2y4 

+  1.047  y.0352  -  ^  ) 


d2yA 

4  )  + 


In  the  soil  strength  test,  the  strength,  S  is 


(10) 


Equation  (9)  can  be  used  to  construct  a  graph  of  S 
in  terms  of  d.  (See  Figure  6). 


Appendix  IV 


MAPS  AND  DRAINAGE  NETS  OF  STUDY  AREAS 


Mt.  Hughes  Quadrangle,  Arizona 


Tributaries  to  Harshaw  Creek 

S  1/2  17,  T  22  S,  R  16  E  (Fourth  order) 
NE  1/4  20,  T  22  S,  R  16  E  (Third  order) 
W  1/2  16,  T  22  S,  R  16  E  (Fourth  order) 


Harshaw  Quadrangle,  Arizona 

Finley  and  Adams  Canyon  I,  II  (Fourth  order) 
Sycamore  Canyon  I,  II  (Fourth  order) 
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Sonoita  Quadrangle,  Arizona 
Hog  Canyon  (Fourth  order) 


Prescott  Quadrangle,  Arizona 
Whipple  Basin  (Fourth  order) 


L 
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Black  Hawk  Quadrangle,  Colorado 

Dory  Hill  Basin  and  Mesa  Gulch  (Fourth  order) 


Morrison  Quadrangle,  Colorado 
Cabrini  Gulch  (Third  order) 


L 
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Netherland  Quadrangle,  Colorado 
Moon  Gulch  (Forth  order) 


Starkville  Quadrangle,  Colorado 

Above:  Tributary  to  Galliilas  Creek  (Third  order) 
Right:  Tributary  to  Saruche  Canyon  (Fourth  order) 
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Allie  Canyon  Quadrangle,  N.M. 
Above: 

Cottonwood  Canyon(Fourth  order) 
Right: 

Tributaries  to  Sapillo  Creek 
E  1/2  20,  T  15  S,  R.  12  W 
(Fourth  order) 

W  1/2  21,  T  15  S,  R  12  W 
(Fourth  order) 


Big  Burro  Mountains  Quadrangle,  New  Mexico 


Tributary  to  Sawmill  Canyon  (Fourth  order) 
Tributary  to  Arrastre  Gulch  (Fourth  order) 


Sawmill  Canyon  (Fourth  order) 

Tributaries  to  Walnut  Creek  I,  II  (Third  order) 
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Cerro  del  Grant  Quadrangle,  New  Mexico 
Cerro  Pavo  Basin  (Third  order) 


Bland  Quadrangle,  New  Mexico 

Above:  Paso  del  Norte  Basin  (Fourth  order) 
Below:  Peters  Dome  Basin  (Fourth  order) 


Hurley  East  Quadrangle,  New  Mexico 
W  1/2  22,  T  18  S,  R  12  W,  I,  II  (Third  order) 
Martin  Canyon 

S  1/2  15,  T  18  S,  R  12  W  (Fourth  order) 

E  1/2  22,  T  18  N,  R  12  W  (Third  order) 
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Jarosa  Quadrangle,  New  Mexico 

Above:  Hairpin  Basin  (Third  order) 
Below:  Poleo  Basin  (Forth  order) 


Paxton  Springs  Quadrangle,  New  Mexico 
Quartz  Hill  Basin  (Fifth  order) 


Above:  Camino  Basin  (Fourth  order) 
Below:  Tano  Basin  (Fourth  order) 


Sedillo  Quadrangle,  New  Mexico 
Cedro  Peak  Basin  (Fourth  order) 
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Twin  Sisters  Quadrangle,  New  Mexico 
Bear  Canyon  I,  II,  III,  (Fourth  order) 


Twin  Sisters  Quadrangle,  New  Mexico 
Wilson  Creek  I,  II,  m,  IV  (Fourth  order) 


Twin  Sisters  Quadrangle,  New  Mexico 

Little  Cherry  Creek  (Fifth  order) 

Mill  Creek  (Fourth  order) 


Lehi  Quadrangle,  Utah 

Hog  Hollow  (Fifth  order) 
Mercer  Hollow  (Fourth  order) 


Springville  Quadrangle,  Utah 
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Above  left:  Spring  Creek  (Fourth  order) 

Above  right:  Deadman’s  Hollow  (Third  obder) 

Tributary  to  Hobble  Creek  (Third  order) 


Spanish  Fork  Peak  Quadrangle,  Utah 
Pole  Canyon  (Fifth  order) 


Appendix  V 


TWO- VARIABLE  CORRELATION  AND  REGRESSION  ANALYSES 
A.  Correlation  tables  of  relations  analyzed  by  grouped-data  methods. 

1.  Loge  Drainage  Density  vs.  P-E  Index. 

2.  Loge  Drainage  Density  vs.  Infiltration  Capacity. 

3.  Loge  Drainage  Density  vs.  Percent  Bare  Area. 

4.  Loge  Drainage  Density  vs.  Valley- side  Slope. 

5.  Log  Channel  Frequency  vs.  Infiltration  Capacity. 

6.  Log  Channel  Frequency  vs.  P-E  Index. 

7.  Log  Channel  Frequency  vs.  Percent  Bare  Area. 

8.  Log  Channel  Frequency  vs.  Log  Wet  Soil  Strength. 

9.  Log  Relative  Relief  vs.  Log  Channel  Frequency. 

10.  Log  Frequency  of  First-order  Channels  vs.  Log  Channel  Frequency. 

11.  Log  Ruggedness  Number  Corrected  for  Order  vs.  Log  Relative  Relief. 

12.  Loge  Density  of  First-order  Channels  vs.  Loge  Drainage  Density. 

13.  Loge  Drainage  Density  vs.  Channel  Frequency. 

14.  Valley-side  Slope  vs.  Infiltration  Capacity,  Fourth-order  Basins. 

15.  Valley-side  Slope  vs.  Log  P-E  Index. 

16.  Valley-side  Slope  vs.  Wet  Soil  Strength. 

17.  Valley-side  Slope  vs.  Log  Channel  Frequency. 

18.  Valley-side  Slope  vs.  Log  Ruggedness  Number  Corrected  for  Order. 

19.  Valley-side  Slope  vs.  Log  Relative  Relief. 
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r=  0.7^;  P<.00|  |r=  —0.43  >  ,oi  ^  P  y .005 

.  Loge Drainage  Density  vs.  Percent  Bare  Area.  4.  Loge  Drainage  Density  vs.  Valley-side  Slope. 
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10.  Log  Frequency  of  First-order  Channels  11.  Log  Ruggedness  Number  Corrected  for  Order 

vs.  Log  Channel  Frequency.  vs.  Log  Relative  Relief. 
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Appendix  V  (Continued) 


B.  Scatter  Plots  and  Regressions  of  Selected  Relations;  Ungrouped  Data. 

1.  Drainage  Density  vs.  P-E  Index. 

2.  Drainage  Density  vs.  Runoff  Intensity. 

3.  Channel  Frequency  vs.  Percent  Bare  Area.  Each  point  represents  one  major 

basin. 

4.  Regression  of  Channel  Frequency  on  Drainage  Density.  Each  point  repre¬ 

sents  one  third-order  basin. 

5.  Regression  of  Valley-side  Slope  on  Infiltration  Capacity. 

6.  Regression  of  Valley-side  Slope  on  Runoff  Intensity. 

7.  Regression  of  Valley-side  Slope  on  Relative  Relief. 

8.  Regression  of  Valley-side  Slope  on  Ruggedness  Number  Corrected  for  Order. 

9.  Regression  of  Valley-side  Slope  on  Wet  Soil  Strength. 

10.  Regression  of  Valley-side  Slope  on  Ratio  of  Mean  January  Precipitation  to 

1/12  Mean  Annual  Precipitation. 

11.  Regression  of  Drainage  Density  on  Percent  Bare  Area. 
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Left  :  1.  Drainage  Density 
vs.  P-E  Index. 

Right:  3.  Channel  Frequency 

vs.  Percent  Bare  Area. 
Each  point  represents  one 
major  basin. 
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Left  :  2.  Drainage  Density 

vs.  Run-off  Intensity  3 

Right:  4.  Regression  of  Channel  |(J 

Frequency  on  Drainage 
Density. 

Each  point  represents  one 
third- order  basin. 
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10.  Regression  of  Valley-side  Slope  on  Ratio  of 
Mean  January  Precipitation  to  1/ 1 2  Mean 
Annual  Precipitation. 
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11.  Regression  of  Drainage  Density  on  Percent 
Bare  Area. 


Appendix  VI  A 

ANALYSES  OF  VARIANCE 

Effect  of  lithology  on  drainage  density:  All  fourth-order  basins,  and  third-order  basins  if  they  are  not 
included  in  a  fourth-order  basin  are  sources  of  data  for  this  study. 


Observed  Drainage  Densities  Grouped  by  Lithology 


Shale 

179.2 

171.2 

350.4 
N  =  2 
X  =  175.2 


Clastic 

Limestone 

Schist 

Granitic 

Acid  Volcanic 

24.1 

14.4 

80.8 

17.4 

15.8 

20.6 

10.0 

86.6 

19.5 

14.1 

11.2 

24.4 

112.7 

15.2 

13.8 

13.1 

22.9 

17.5 

14.9 

15.1 

N  =  2 

10.4 

16.4 

14.4 

16.2 

313.4 

11.0 

9.9 

12.8 

X  =  12.2 

5.7 

4.4 

11.4 

N  =  5 

7.5 

13.1 

33.5 

7.5 

15.0 

17.1 

X  =  62.7 

5.0 

11.9 

17.1 

14.8 

9.2 

14.7 

13.7 

13.8 

12.8 

14.9 

12.9 

219.7 

17.1 

13.2 

15.9 

12.1 

36.9 

11.2 

N  =  13 

58.1 

13.2 

_ 

39.3 

12.6 

X  =  1  b .  9 

32.2 

13.2 

39.9 

13.1 

32.9 

12.0 

438.4 

12.3 

10.4 

N  =  21 

286.5 

X  =  20.9 

N  =  23 

X  =  12.5 
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N 

Sx. 

X 

XI 

1 

XI 

EX^ 

(sx£)  2/n 

S  2 

Shale 

2 

350.4 

175.2 

150.5 

61,422.08 

61,390.08 

32.00 

Schist 

5 

313.4 

62.7 

37.9 

27,362.06 

19,643.91 

1,929.54 

Clastic 

13 

219.7 

16.9 

-7.8 

4,173.47 

3,712.93 

38.38 

Granitic 

21 

438.4 

20.9 

-3.9 

12,947.18 

9,152.12 

189.75 

Limestone 

2 

24.4 

12.2 

-12.5 

307.36 

297.68 

9.68 

Acid  Volcanic 

23 

286.5 

12.5 

-12.3 

3,693.73 

3,568.79 

5.68 

Totals 

66 

1,632.8 

-- 

151.9 

109,905.88 

97,765.51 

— 

(ssx ij  )  2/n  - 

40,394.48 

Analysis- of- Variance  Table 

Source  of  Variation 

S.  S. 

d.f. 

M.  S. 

F 

Among  column  means 

57,371.03 

5 

11,474.21 

56.7 

Within  columns 

12,140.37 

60 

202.34 

Total 

69,511.40 

65 

1,069.41 

Ft  01(5,  60)  =  3.339 

Therefore,  a  significant  difference  among  means  of  drainage  density  for  each  lithology  exists  at  the 
.01  level  of  significance,  and  also  at  the  .001  level. 

Contrasts*  judged  significant  at  the  .05  level: 

(Granitic,  clastic,  volcanic,  and  limestone)  vs. 

(shale  and  schist)  (-97.97  <0<  -58.85) 

Shale  vs.  Schist  (71.58  <  0  <  153.46) 

Contrasts  judged  not  significant  at  the  .05  level: 

Clastic  vs.  granitic.  (-21.24  <  ©<  29.20)  - 

Limestone  vs.  granitic.  (-44.92  <  0<  27.52) 

Limestone  vs.  volcanic.  (-35.33  <  0<  36.81) 

Clastic  and  granitic  vs.  limestone  and  volcanic 

(-5.98  <  0<  19.82) 

Effect  of  lithology  on  infiltration  capacity:  Individual  measurements  of  infiltration  capacity  grouped  ac¬ 
cording  to  lithology,  regardless  of  basin  or  climate,  are  sources  of  data.  (Individual  infiltration  meas¬ 
urements,  grouped  by  basin,  are  given  in  Table  7.) 


*A  contrast  is  a  difference  in  means  or  averages  of  means  of  several  groups.  In  Scheffe’s  notation,  the  true  (popula¬ 
tion)  contrast  is  denoted  by  0 ,  and  the  estimated  contrast  calculated  from  sample  data  is  denoted  by  (3.  His  notation 
is  used  in  this  section,  and  must  not  be  confused  with  the  identical  notation,  0,  used  for  valley-side  slopes  in  the 
previous  sections. 
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N 

SXi 

X 

X  -  X  .. 

XX .  2 

(SXi)2/N 

s  2 

Shale 

4 

4.4 

1.1 

-1.4 

4.86 

4.84 

0.007 

Schist 

7 

10.7 

1.5 

-1.0 

19.53 

16.36 

0.53 

Clastic 

53 

106.1 

2.0 

-0.5 

257.07 

212.40 

0.86 

Granitic 

52 

173.7 

3.3 

0.8 

744.01 

580.22 

3.21 

Limestone 

8 

24.4 

3.0 

0.5 

96.00 

74.42 

3.08 

Acid  Volcanic 

17 

37.8 

2.2 

-0.3 

100.22 

84.05 

1.01 

Totals 

141 

357.1 

-1.9 

1,221.69 

972.29 

(SLXij^/N  = 

904.40 

Analysis-of-Variance  Table 

Source  of  Variation 

S.  S. 

d.f. 

M.  S. 

F 

Among  column  means 

67.89 

5 

13.58 

7.43 

Within  columns 

249.4 

13  5 

1.85 

Total 

317.29 

140 

F>001(5,  135)  =  4.42. 

A  significant  difference  exists  among  the  infiltration  means  for  the  various  types  of  lithologies  at  the  .001 
level  of  significance. 

The  contrast  (granitic  and  limestone)  vs.  (clastic  and  volcanic),  formed  from  the  groups  with  highest 
and  intermediate  infiltration  capacities,  was  judged  significant  at  the  .05  level  (0.450  <  0  <  2.042,  in 
Scheff^’s  notation).  The  contrast  (clastic  and  volcanic)  vs.  (shale  and  schist),  formed  from  the  groups  with 
intermediate  and  lowest  infiltration  capacities,  was  judged  not  significant  at  the  .05  level  (-0.780  <  ©  < 
2.160). 

Effect  of  lithology  on  valley- side  slope:  Individual  field  measurements  of  slopes,  separated  according 
to  lithologic  type  irrespective  of  basin  boundaries,  are  sources  of  data. 


N 

XX, 

X  X 

-  X.. 

XX . 2 

1 

(XX^2  /N 

s2 

Shale  82 

2,390.2 

29.2 

7.9 

74,687.75 

69,671.41 

61.93 

Schist  113 

1,898.0 

16.8 

-4.4 

38,135.50 

31,879.68 

55.85 

Clastic  316 

7,162.5 

22.7 

1.5 

193,572.75 

162,346.22 

99.13 

Granitic  226 

4,495.0 

19.9 

-1.3 

149,979.00 

89,402.76 

269.22 

Limestone  20 

607.0 

30.4 

9.1 

18,594.00 

18,422.45 

9.02 

Acid  Volcanic  123 

2,128.0 

17.3 

-3.9 

43,868.00 

36,816.13 

57.80 

Totals  880 

18,680.7 

8.9 

518,837.00 

408,538.65 

(XXX^)  2/N  =  396,555.17 

Analysis-of-Variance  Table 

Source  of  Variation 

S.  S. 

d.f. 

M.  S. 

F 

Among  groups 

11,983.48 

5 

2,396.69 

18.99 

Within  groups 

110,298.35 

874 

126.19 

Total 

122,281.83 

879 

F .  o  o  i  °°)  =  4.10. 
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A  significant  difference  among  the  valley-side  slope  means  for  each  lithology  exists  at  the  .001  level. 

The  contrast  obtained  from  the  four  lithologies  with  most  nearly  similar  slopes  (granitic  and  clastic) 
vs.  (schist  and  volcanic),  is  significantly  different  from  zero  at  the  .05  level  (1.91  <  0  <  7.75). 


Appendix  VI  B 


DETERMINATION  OF  THE  EFFECT  OF  GULLYING  AND  TRENCHING  ON  DRAINAGE  DENSITY,  BY 
THE  ANALYSIS  OF  COVARIANCE 


Source  of  data:  drainage  density  and  mean  infiltration  capacity  for  each  basin  for  which  data  on  chan¬ 
nel  condition  are  available. 


1  2  3 

Channels  trenched  Slopes  gullied  Unmodified  basins 


Xi 

Yi 

x2 

y2 

x3 

Y3 

f 

LogeD 

f 

Loge  D 

f 

Loge  D 

3.05 

2.303 

2.60 

3.183 

2.74 

3.428 

1.99 

2.306 

2.75 

1.609 

1.10 

5.142 

2.23 

3.131 

1.65 

2.798 

2.44 

2.292 

1.35 

3.026 

1.95 

2.550 

3.70 

2.015 

1.88 

3.447 

1.70 

2.839 

1.10 

5.187 

1.90 

3.512 

1.00 

4.392 

2.15 

2.550 

1.00 

4.462 

2.10 

1.482 

1.00 

4.727 

1.00 

4.180 

1.78 

2.839 

Total 

12.4 

17.725 

16.43 

33.579 

15.33 

22.096 

Mean 

2.07 

2.954 

1.64 

3.358 

2.19 

3.156 

N 

6 

6 

10 

10 

7 

7 

EEXir 

=  44.16,  X.. 

=  1.92. 

zzxi3-- 

=  73.400,  Y.. 

=  3.191. 

N  = 

23. 

sxL  = 

27.20, 

(sxu) 

=  25.63, 

EXilYil  = 

35.83,  sY2x 

=  53.80, 

(sYipVn,  = 

52.36. 

^X-2  = 

30.91, 

(SXi2)2/n2 

=  26.99, 

exi2yi2  = 

49.93,  SY22 

=  122,19, 

(sYi2)2/n2  = 

112.75. 

SX,23  = 

38.60, 

<SXi3)7n3 

=  33.57, 

SXi3Yi3  = 

42.397,  SY23 

=  83.11, 

(SYi3)2/n3  = 

69.75. 

ssx2. 

1J 

=  96.71 

,  SSY ?.  =  259.10,  S((SX..  )2/n.)  =  86.19, 

-j  -J  J  J 

S(  (S  Y  ij )  2/n .  =  234.86, 

(ssx^)2/ N  = 

84.79,  (SZY..)2/N  =  234J.24.1 

Total  sum  of  products  =  -12.77. 

Within-groups  sum  of  products  =  -12.03. 

Among-means  sum  of  products  =  -0.74. 

Total  sum  of  squares:  for  X’s  11.92,  for  Y’s  24.86. 
Among-means  sum  of  squares:  for  X’s  1.4,  for  Y’s  0.42. 
Within-groups  sum  of  squares:  for  X’s  10.52,  for  Y’s  24.44. 
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Sum  of  squares  about  the  regression  line: 

Total  S.  S.  =  24.86  -  (~12-77^  =  11.18. 

11.92 

Within-groups  S.  S.  =  24.44  -  ^ =  10-68. 

The  dispersion  of  the  Y  values  about  the  regression  line  with  the  coefficient  bt  (total)  is  therefore  11.18. 
The  dispersion  of  the  Y  values  about  the  regression  line  with  coefficient  bw  (within)  is  therefore  10.68. 
The  reduction  in  dispersion  is  0.50,  attributable  to  dispersion  of  the  means.  The  low  value  of  F  signifies 
that  the  reduction  is  not  great  enough  to  be  judged  significant. 


Analysis-of-Covariance  Table 


Source  of  Variation 

d.f. 

Ex2 

Sxy 

sy 2 

d.f. 

sy'2 

M.S. 

F 

Among  means 

2 

1.4 

-  0.74 

0.42 

2 

0.50 

0.25 

0.45 

Within  groups 

20 

10.5 

-12.03 

24.44 

19 

10.68 

0.56 

Total 

22 

11.9 

-12.77 

24.86 

21 

11.18 

' 


Appendix  VII 


Line 

Desigroj 

tion 

Column  Designation 

at, 

Oii 

Ck3 

a* 

ai5 

9t 

Check 

Normc 

jl  E 

.quatlo 

ns 

i 

n 

in 

u 

JZ 


2.031219 

.130534 

•1.767067 

-1.476415 

-.467847 


.130554 
.653(60 
'.9402791  (. 
-.406014 
i3?!374l 


1.767067 

-.940279 

.883174 

1-043488 

.41502.1 


-1.42841! 

-.406014 

1.0434*8 

9.950TO 

•368794 


!5hA67847 

-.321524 

.4(5012. 

.368794 

.192702 


-2.950!% 
-1.406244 
2.712221 
1. 113502 
HiS9i5. 


-3.651832 

-(410217 

3.346559 

10650200 

0.82062 


5»‘  *  _ 

4.818804 


Front  Solution 


I 

r 

3T 

n-i 

r 

m 

m-i 

m-2 

nr 

uc 

32-1 

ET-Z 

rz-3 

& 

EC' 

"2 

32- 1 
Y-2 
Y-3 
1-4 

^3 

1 Z' 

2.0312.19 

-1.0OO0OO 

.930554 

-.438126 

.633160 

-.426311 

-1.767067 

♦WW4 

-310273 

-.809539 

-1428495 
f. 7052  70 

-406014 

.654431 

-.467847 

♦.230318 

-.321324 

.2(433  3 

-2-950196 

I.45M& 

-i.iofciai 

L35I56I 

-3.651832 

1. 717*32 

1.490297 

L672396 

.226846 

-1.000000 

-.(30740 

♦.576330 

1.885174 

-1537267 

-.075349 

.2484(7 

4.095076 

1.043488 

-1.242725 

.143170 

-I0«J9I 

*.471640 

.415022 

-.407005 

r-061662 

-.054833 

*.241716 

2.712221 

-2.566535 

-.031602 

.182702 

-.805390 

3.34*359 

-3176926 

.105297 

.270558 

HOOOOOO 

-.056067 

.207227 

9.958933 

-1.004618 

-.272035 

-.011619 

-.053645 

.198275 

.368794 

-.329023 

JII7I63 

-.011117 

.1 14084- 
-.42IG62 

1.113502 

-2.074784 

.060046 

.023641 

.27493C 

1.0IGI53 

I0.65020( 

£568224 

-.200073 

.056973 

8670661 

-LOOOOOC 

.145817 

1016817 

.192702 
-.(0/ /58 
-050461 
i0l0636 
-.002452 

-.877595 

.101214 

.633915 

-.67950 
-.025861 
.022620 
.01475  8 

7.938864 

-915605 

£11262 

-841119 

.086170 

.054512 

-.133508 

.021395 

1.000000 

-.03408 

159294 

-012683 

.592802. 

Back  Solution 

Eq-Tf 

Check 

9t 

Compu- 

tationof 

SSR-Z.b.1' 

b, 

b. 

bs 

b* 

b5 

1.452426 

-.366899 

-.052542 

.078643 

.393717 

1.399307 

1241716 

1751294 

.081502 

7)52100 

-.859408 

►.421662 

1315840 

10(5423 

.090399 

-.101212 

.026786 

1074426 

-1.59294 

159294 

\9Z70l 

.368794 

.415022 

-.321324 

-.467847 

-.306963 

-.027448 

.037518 

.276148 

.654662 

.633915 

1113502 

2.71222 

-1.40639) 

-2.950196 

-1.009788 

-.082873 

.2451*2. 

1.208644 

4.128230 

EqT -.633915 

.633917 

4.4894 

f*  =0.9317  55E  =  S*‘ -  55R=  0.3292 


R=  .9652 
a  =  110.8099 


3*-  SSE/N-tfO  -0.01934 
s-0.1391 


Solution  of  Normal  Equations  Relating  Drainage 
Density  to  Soil  Properties  and  Climate  Elements 
by  the  Doolittle  Method. 
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TABLE  1 


DEFINITIONS  OF  SYMBOLS 


Measurements  from  maps. 

Units 

Dimensions 

lu 

Total  length  of  u’th  order  channels; 

Miles 

L 

21U 

Total  length  of  channels  of  all  orders; 

Miles 

L 

P 

Perimeter  of  basin; 

Miles 

L 

r 

Total  relief  of  basin; 

Feet 

L 

z 

Elevation  of  mouth; 

Feet 

L 

a 

Area  of  basin  (planimetric); 

Square  miles 

L2 

“u 

Number  of  u’th  order  channels; 

Enumerative 

0 

Snu 

Total  number  of  channels  of  all  orders; 

Enumerative 

0 

nu 

Minimum  possible  number  of  channels  of  order 
u  in  a  given  basin; 

Enumerative 

0 

u 

Order  of  channel  or  basin; 

Enumerative 

0 

»/> 

Inflection  angle  of  contour  lines  in  valleys; 

Degrees 

0 

£ 

Axil  angle  of  tributaries  (angle  of  entry); 

Degrees 

0 

B.  Quantities  derived  from  map  measurements. 

Y 1  • 

D  Drainage  density, - -’ 


D1  Density  of  first-order  channels, 


U  ; 


F  Stream  frequency, 


Sn,,; 


Fu  Frequency  of  u’th  order 

3. 

D#  r 

H  Ruggedness  number,  ’ 

R  Relative  relief,  ’ 

’  5280p 

P  Ratio  of  channel  lengths  to  perimeter, 


siu; 


C  Circularity  of  basin 


B  u :  u+  i  Bifurcation  ratio,  - si—’ 

nntl 

Lu;u_1  Length  ratio,  ’ 


area  circle  with  same  p  ’ 


n 


1 


u-l 


Miles 
Sq.  miles 

Miles 
Sq.  miles 

Number 
Sq.  miles 

Number 
Sq.  miles 


L-1 


L  '2 


0 

0 

0 

0 

0 

0 


r 
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e 

Maximum  valley-side  slopes; 

Degrees 

0 

Dihedral  angle  between  valley  sides; 

Degrees 

0 

y 

Gradient  of  stream  channel; 

Degrees 

0 

Field  measurements. 

f 

Infiltration  capacity  of  soil; 

Inches/hr. 

L  T"l 

s 

Strength  of  soil  under  instantaneous  impact; 

Lbs./sq.  ft. 

M  L-lrT-2 

b 

Percent  of  area  that  is  bare  of  all  cover; 

Percent 

0 

c 

Number  of  rock  fragments  in  a  circle  of  one  foot 
radius; 

Enumerative 

0 

m 

Average  mid-diameter  of  rock  fragments; 

Inches 

L 

M 

Roughness  number,  average  of  total  length  of 
pebble  diameters  for  each  one-foot  circle  in  a 
basin; 

Inches 

L 

e 

Maximum  valley-side  slopes; 

Degrees 

0 

Climatic  Data. 

P 

Mean  precipitation,  either  monthly  or  annual; 

Inches 

L 

T 

Monthly  average  mean  daily  temperatures; 

Degrees 

6 

I 

Precipitation-  Effectiveness  index; 

- 

0 

q 

Runoff  intensity:  5-year,  1-hour  rainfall  amount- 
infiltration  capacity  +  5; 

Inches /hr. 

LT'1 

LIST  OF  MORPHOMETRIC  PROPERTIES  OF  BASINS 


Table  2 


ARIZONA 


Cameron 

r 

z 

U 

0 

P 

1, 

1? 

U 

U 

U 

<U<1 

n, 

n? 

n, 

n+ 

ns 

l , 

la 

u 

U 

Bi  p 

Bpt 

B3.-4 

04-5 

L  2ll 

L  3.2 

L-4'3 

1-5:4 

D 

D, 

F 

F| 

H 

R 

P 

C 

Y 

e 

CHINLE  BADLANDS,  N. 

31 

— 

4 

300903 

1436 

0783 

0356 

0018 

0462 

1618 

23 

7 

2 

i 

33 

00340 

00508 

00090 

0462 

3.2857 

3.50 

2.0 

- 

1.49 

.177 

5/  3 

- 

179  2 

86  71 

36.540 

25,470 

1  07 

4.15 

1  13 

55 

21  2 

CHINLE  BADLANDS.  N.  I 

10 

— 

3  . 

000069 

0356 

0114 

0118 

000 928 

- 

- 

0241 

6 

2 

1 

- 

- 

9 

0019 

00590 

000928 

- 

_ 

3. 

20 

- 

- 

3  10 

.157 

- 

- 

34  9 

165.2 

130,400 

86,960 

66/ 

5  32 

.677 

69 

212 

CHINLE  BADLANDS,  N.  U 

10 

— 

3 

700075 

0364 

01093 

01083 

00097 

_ 

_ 

02273 

5 

2 

1 

- 

- 

8 

00219 

00542 

000 97 

_ 

_ 

2.5 

2  0 

_ 

_ 

247 

.179 

- 

- 

303 

145.7 

106,700 

66,670 

•574 

5.20 

624 

68 

~ 

21  2 

CHINLE  BADLANDS,  5. 

37.5 

— 

3 

OOIII 

1544 

1202 

0374 

0323 

- 

- 

1900 

19 

4 

1 

- 

- 

24 

00533 

00935 

0323 

- 

- 

47  5 

4.0 

- 

- 

1  48 

3  45 

- 

- 

171 

108-3 

21,620 

17,120 

1.22 

460 

1  23 

58 

31.1 

Harshaw  Quadrangle 

FINLEY  AND  ADAMS  CANYON  I 

760 

5490 

4 

143 

1  704 

1.515 

.379 

.360 

227 

- 

2. 481 

20 

7 

2 

l 

- 

30 

0758 

0541 

180 

22  7 

- 

2  8571 

350 

2  0 

- 

714 

3  33 

126 

- 

17.  35 

1059 

209.8 

139  9 

2.49b 

64  5 

/  457 

.  62 

133 

31  8 

FINLEY  AN  ADAMS  CANYON  n 

510 

54  90 

4 

.665 

4.09 

8  790 

2.312 

1.554 

303 

- 

12.96 

98 

14 

3 

/ 

- 

116 

■0897 

.165 

■518 

.303 

_ 

7.0 

4-666 

3.0 

_ 

1  84 

3  14 

.585 

- 

19  48 

13  22 

174. 4 

1474 

5.5  7 

6  99 

3.165 

50 

113 

30  0 

SYCAMORE  CANYON  I 

655 

5245 

4 

2905 

2  58 

2.80 

.72 

45 

45 

- 

4- 42 

21 

6 

2 

I 

- 

30 

.133 

.1 20 

225 

45 

- 

3.5 

3.00 

2.0 

- 

.9 02 

1.88 

2.00 

- 

15.22 

9  638 

103  3 

72.28 

4.77 

1 2  15 

1.71 

55 

III 

31.9 

SYCAMORE  CANYON  U 

655 

5245 

4 

3125 

2  56 

3.45 

14 

■83 

.04 

5  46 

34 

8 

2 

/ 

- 

45 

.101 

142 

415 

■  04 

- 

4.25 

4-00 

2  0 

1.40 

2  92 

.096+ 

~ 

17.47 

II  04 

144  0 

100  8 

5  48 

12  24 

2  13 

60 

108 

33.7 

Mt.  Hughes  Quadrangle 

TRIBUTARY  TD  HARSHAW  CREEK 

5  17  T22S,  RIG  E 

1691 

4 450 

3 

.312 

2.765 

3  1 44 

492 

.758 

4- 394 

26 

3 

1 

30 

.121 

164 

.758 

86666 

3.0 

1.36 

462 

14  003 

10  08 

96  15 

83.33 

5  044 

12  95 

1  59 

51 

107 

37.0 

TRIBUTARY  TO  HARSHAW  CREEK 

N£ .  Va  20  7  22  3,  R  lb  E 

TRIBUTARY  TO  HARSHAW  CREEK 
W/2  16  T 22  S.  RI6E 

1475 

2144 

4450 

4197 

3 

4 

.1496 

.628 

1.837 

4-/48 

1269 

6.761 

606 

.458 

189 

.193 

530 

- 

2  064 

9  942 

15 

57 

2 

If 

/ 

2 

/ 

- 

18 

71 

0846 

.1/9 

.303 

.132 

18  9 

556 

.530 

_ 

75 

5.1818 

2  0 

5.50 

2.0 

_ 

358 

i.n 

.6  24 

452 

.889 

: 

13  800 

15631 

8483 

10  76 

120-3 

113  0 

1003 

90  76 

3.855 

6-428 

1521 

9  79 

/  15 

2  40 

57 

46 

107 

IOO 

37.0 

350 

Prescott  Quadrangle 

WHIPPLE  BASIN 

325 

5225 

4 

16  2 

150 

158 

68 

.10 

-22 

- 

2  66 

25 

8 

2 

/ 

36 

0632 

0850 

.09 

.22 

- 

3.125 

4.0 

20 

- 

1.34 

1.06 

244 

- 

16.42 

9.753 

222,2 

154  3 

(.01 

4  10 

1-77 

91 

21.1 

WHIPPLE  I 

275 

5275 

J 

.05/6 

■  928 

.6  20 

.252 

.130 

- 

- 

1.002 

12 

5 

1 

- 

- 

18 

051 7 

■0504 

.130 

_ 

- 

2.40 

50 

- 

- 

.975 

2  58 

- 

- 

19.34 

11.97 

347  5 

23/7 

(.007 

5  61 

100 

.76 

~ 

2/  1 

WHIPPLE  n 

235 

5275 

3 

.0595 

.961 

.521 

269 

■01 6 

- 

- 

•  006 

8 

2 

1 

- 

- 

II 

0651 

134 

.016 

- 

4.0 

2-0 

~ 

2.06 

■  1 19 

- 

13  55 

8.756 

104  9 

134  4 

,603 

463 

.039 

81 

" 
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Sonoita  Quadrangle 

HOG  CANYON 

515 

4760 

4 

22/ 

275 

2  90 

1.73 

20 

49 

- 

5- 32 

98 

15 

3 

/ 

- 

II 7 

0296 

115 

-0666 

49 

- 

6  5333 

50 

3.0 

- 

3.88 

.579 

7.36 

- 

24  07 

1 3/2 

529.4- 

443  4 

235 

3.55 
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37 
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HOG  CANYON  X 
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4865 
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.0056 

32 

1 25 
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- 

- 

24 

6 

2 

/ 

- 

9 
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.022 
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- 

_ 

3.00 

20 

_ 

- 
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- 

- 
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1607 
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68 
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- 

36 
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_ 

- 
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_ 

- 
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- 

- 

10  28 
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40 
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3 

0440 
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60 
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- 
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22 

2 

1 

- 

25 

.0273 
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- 

- 
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20 

- 
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- 
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Tucson  Quadrangle 
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- 
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- 
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- 
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3 
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II 

00414 

-00265 

.0148 

- 

- 
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-640 
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28,950 

21,050 
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TABLE  3 

VEGETATION  CHARACTERISTICS 


Quadrangle  and  Basin  Measured 

Range  of  Elevation  Vegetation  Species  Bare  Area(%) 


ARIZONA 

Chinle  Badlands 
about  4000 

Harshaw  Quadrangle 
Finley  and  Adams 
Canyon 
7000-5490 
and 

Sycamore  Canyon 
6900-5245 

Mt.  Hughes  Quadrangle 
Tributaries  to 
Harshaw  Creek 
6341-4200 


Prescott  Quadrangle 
Whipple  Basin 
5550-5225 


Sonoita  Quadrangle 
Tributary  to 
Hog  Canyon 
5275-4760 


Tucson  Quadrangle 
Saguaro  National 
Monument 
about  2800 


COLORADO 

Black  Hawk  Quadrangle 
Mesa  Gulch  and 
Dory  Hill  Basin 
9651-8300 

Morrison  Quadrangle 
Cabrini  Gulch 
7165-6650 


Green  Mountain 
Basin 

6603-6090 


Nederland  Quadrangle 
Delonde  Gulch 
10215-8580 


No  vegetation. 


100 


South-facing:  shrub  live  oak  chaparral, 
manzanita. 

North-facing:  emory  oak,  silverleaf  oak,  netleaf 
oak,  Arizona  white  oak,  alligator  juniper,  manzanita, 
madrone,  Mexican  pinon,  Mexican  blue  oak, 
abundant  grass. 


Below  4800:  catclaw  acacia. 

4800-5800:  shrub  live  oak  chaparral,  manzanita, 
ocotilla,  yucca. 

Above  5800:  Mexican  blue  oak,  Arizona  white 
oak,  alligator  juniper,  Mexican  pinon. 


Arizona  white  oak,  shrub  live  oak,  one-seed  56 

juniper,  alligator  juniper,  hairy  mountain  ma¬ 
hogany,  singleleaf  pinon,  ponderosa  pine  (small), 
cliffrose,  cholla. 

South -facing:  mainly  grass.  38 

North-facing:  emory  oak,  Arizona  white  oak,  3 

Mexican  blue  oak,  century  plant,  jumping  cholla, 
manzanita,  alligator  juniper. 

Saguaro,  mesquite,  paloverde,  ocotilla,  jumping  83 

cholla,  staghorn  cholla,  creosote  bush,  prickly 
pear,  fishhook  barrel  cactus,  brittle  bush,  Christ¬ 
mas  cactus,  sparse  bunch  grass. 


Lodgepole  pine,  blue  spruce,  Douglas-fir. 


Grass,  sage,  yucca,  poison  ivy,  wild  rose,  juni-  40 

per,  prickly  pear,  chokecherry,  hairy  mountain 
mahogany,  gooseberry,  bigtooth  maple,  barrel 
cactus,  Russian  thistle. 

Birchleaf  mountain  mahogany,  prickly  pear,  50.2 

grass,  sage,  yucca,  common  chokecherry  (on 
north  slopes  and  gullies),  barrel  cactus,  sun¬ 
flowers,  gooseberry. 


Aspen,  lodgepole  pine,  ponderosa  pine,  limber 
pine,  white  fir,  blue  spruce,  mountain  common 
juniper,  grass,  Douglas-fir,  willows. 


Life  Zone(s) 


Upper  Sonoran 
Upper  Sonoran 

Upper  Sonoran 

Lower  Sonoran 
Upper  Sonoran 

Upper  Sonoran 

Upper  Sonoran 

Lower  Sonoran 
Upper  Sonoran 

Lower  Sonoran 


Canadian 

Upper  Sonoran 

Upper  Sonoran 

Canadian  to 
Hudsoni  an 
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TABLE  3  (2) 

Quadrangle  and  Basin 

Range  of  Elevation 

Vegetation  Species 

Measured 

Bare  Area(%) 

Life  Zone(s) 

COLORADO 


Hicks  Gulch 

10210-8640 

Aspen,  lodgepole  pine,  ponderosa  pine,  limber 
pine,  white  fir,  blue  spruce,  mountain  common 
juniper,  grass,  Douglas-fir,  willows. 

9.3 

Canadian  to 
Hudsonian 

Moon  Gulch 

10625-8415 

Lodgepole  pine,  limber  pine,  white  fir,  blue 
spruce,  aspen,  mountain  common  juniper,  sage, 
bunch  grass,  kinnikinnik. 

27.8 

Canadian  to 
Hudsonian 

Starkville  Quadrangle 
Tributaries  to 

Gallinas  Creek 
7590-6920 
and 

Tributaries  to 

Saruche  Canyon 
7625-6875 

Pinon-juniper  woodland. 

- 

Upper  Sonoran 

NEW  MEXICO 

Big  Burro  Mountains 
Quadrangle 
Tributaries  to 

Sawmill  Canyon 
7755-6745 

Shrub  live  oak,  Gambel  oak,  alligator  juniper, 
ponderosa  pine,  yucca,  barrel  cactus,  bunch 
grass,  prickly  pear,  cholla,  Douglas-fir. 
South-facing:  shrub  live  oak  chaparral. 

North-facing:  ponderosa-pine  forest. 

- 

Transition 

Tributary  to 

Arrastre  Gulch 
8035-6745  and 
Sawmill  Canyon 
7820-6280 

South-facing:  shrub  live  oak  chaparral,  Gambel 
oak,  alligator  juniper,  yucca,  barrel  cactus, 
bunch  grass,  prickly  pear,  cholla. 

North-facing:  ponderosa  pine,  Douglas-fir,  shrub 
live  oak,  Gambel  oak,  alligator  juniper. 

Upper  Sonoran 
to  Transition 

Tributaries  to 

Walnut  Creek 

8035-7188 

Ponderosa  pine,  shrub  live  oak,  Gambel  oak, 
Douglas-fir,  walnut,  alligator  juniper,  prickly 
pear,  bunch  grass,  cholla. 

Transition 

Bland  Quadrangle 

Paso  del  Norte 

Basin 

9470-8520  and 

Peters  Dome  Basin 
8460-7690 

Ponderosa  pine,  shrub  Gambel  oak,  white  fir, 
Douglas-fir,  aspen,  Rocky  Mountain  juniper, 
pinon,  limber  pine,  wavy  leaf  oak,  grass. 
North-facing:  ponderosa-pine  and  lower  Douglas- 
fir  forest. 

South-facing:  ponderosa-pine  forest. 

3.3 

Canadian 

Transition  to 
Canadian 

Cerro  del  Grant 
Quadrangle 

Cerro  Pavo  Basin 
10310-9460 

Blue  spruce,  Douglas-fir,  white  fir,  Engelmann 
spruce,  ponderosa  pine,  mountain  common  juni¬ 
per,  one-seed  juniper,  aspen,  grass,  alpine  fir. 

13.2 

Canadian  to 
Hudsonian 

Jarosa  Quadrangle 
Hairpin  Basin 
8665-8295 

Ponderosa  pine,  shrub  Gambel  oak,  aspen,  white 
fir,  birchleaf  mountain  mahogany,  pinon,  locust, 
grass. 

12.8 

Transition 

Poleo  Basin 

9183-8038 

Ponderosa  pine,  shrub  Gambel  oak,  aspen,  white 
fir,  Douglas-fir,  Rocky  Mountain  juniper,  birch¬ 
leaf  mountain  mahogany,  mountain  common  juni¬ 
per,  pinon. 

25.8 

Transition  to 
Canadian 

TABLE  3  (3) 
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Quadrangle  and  Basin 
Range  of  Elevation 

Vegetation  Species 

Measured 

Bare  Area(%) 

Life  Zone(s) 

NEW  MEXICO 

Paxton  Springs 
Quadrangle 

Quartz  Hill  Basin 
8033-7480 

Ponderosa  pine,  pinon,  one-seed  juniper,  Rocky 
Mountain  juniper,  wavy  leaf  oak,  Gambel  oak, 
yucca,  grass,  mountain  mahogany. 

54.1 

Transition 

Santa  Fe  Quadrangle 
Camino  Basin 

6890-6740 

Pinon,  one-seed  juniper,  rabbit  brush,  prickly 
pear,  yucca,  cholla,  grass. 

60 

Upper  Sonoran 

Tano  Basin 

7331-7050 

Pinon,  one-seed  juniper,  prickly  pear,  grass, 
rabbit  brush,  yucca,  hairy  mountain  mahogany. 

60.8 

Upper  Sonoran 

Sedillo  Quadrangle 

Cedro  Peak  Basin 
7767-6525 

Pinon,  one-seed  juniper,  wavy  leaf  oak,  Gambel 
oak  (shrubs  and  small  trees),  Rocky  Mountain 
juniper,  cholla,  prickly  pear,  birchleaf  mountain 
mahogany,  rabbit  brush,  yucca,  grass,  (small) 
barrel  cactus,  Russian  thistle. 

39.6 

Upper  Sonoran 

Twin  Sisters  Quadrangle 
Wilson  Creek 

8020-6862 

North-facing:  ponderosa  pine  and  alligator  juni¬ 
per. 

South-facing:  juniper-Gambel  oak  woodland, 
bunch  grass. 

- 

Transition 

UTAH 

Lehi  Quadrangle 

Dry  Hollow 

6682-5340 

Shrub  Gambel  oak,  sage,  bigtooth  maple,  grass, 
prickly  pear,  cliffrose,  birchleaf  mountain 
mahogany . 

42.8 

Upper  Sonoran 

Hog  Hollow 

6350-5135 

Gambel  oak,  bigtooth  maple. 

- 

Transition 

Maple  Hollow 

6682-5095 

Gambel  oak,  bigtooth  maple. 

Transition 

Mercer  Hollow 
6350-5120 

Gambel  oak,  bigtooth  maple. 

- 

Transition 

Springville  Quadrangle 
Deadmans  Hollow 
6250-4960 

Gambel  oak  (small  trees  and  shrubs),  sparse 
bunch  grass,  sage.  Northslopes;  bigtooth  maple, 
netleaf  hackberry. 

- 

Upper  Sonoran 
to  Transition 

Spring  Creek 

9325-6050 

Gambel  oak  (trees  and  shrubs),  bigtooth  maple 
(trees  and  bushes),  aspen,  common  chokecherry, 
curleaf  and  birchleaf  mountain  mahogany,  bunch 
grass,  Utah  serviceberry. 

28.9 

Transition 

Tributary  to 

Hobble  Creek 

7320-4920 

Gambel  oak  (small  trees  and  shrubs)  sparse 
bunch  grass,  sage,  Northeast  slopes:  bigtooth 
maple,  netleaf  hackberry. 

40* 

Upper  Sonoran 
to  Transition 

♦Indicates  estimate 
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TABLE  4 


SOIL  AND  SURFACE  PROPERTIES 


Infiltration 

Soil  Strength, 

Roughness 

State,  Quadrangle, 

Capacity, 

lbs. /sq.ft. 

3. 

Number, 

Basin 

Inches/hr. 

Wet  Dry 

Soil  Description 

Inches 

ARIZONA 
Cameron 
Chinle  Badlands 


North 

Chinle  Badlands 

1.1 

500 

15,600 

Blue-gray  silt  and  clay, 
trace  coarse  to  fine  sand; 

3.99 

South 

Harshaw  Quadrangle 

1.1 

500 

16,960 

very  plastic. 

0 

Finley  and  Adams 

Canyon 

Prescott  Quadrangle 

3.80 

Sieve  analysis:  91.1%  sand, 
median  at  .373  mm.  Coarse 
to  fine  (+)  sand,  little  med¬ 
ium  to  fine  (+)  gravel,  trace 
(+)  silt. 

Whipple  Basin 

Sonoita  Quadrangle 

Tributary  to 

1.65 

9270 

24,900b 

Coarse  to  fine  (+)  sand, 
trace  coarse  to  fine  gravel, 
trace  silt.  Residual  granite 
soil. 

34.4 

Hog  Canyon 

Tucson  Quadrangle 

Saguaro  National 

2.6 

12,000 

46,300 

Coarse  to  fine  sand  88%, 
little(-)  fine  gravel,  trace 
silt. 

45.9 

Monument,  North 

Saguaro  National 

1.0 

14,000 

44,300 

Angular  coarse  to  fine  sand, 
little  (+)  medium  to  fine 
gravel,  trace(+)  silt  and  clay. 

83.5 

Monument,  South 

COLORADO 

Morrison  Quadrangle 

1.0 

10,750 

54,300 

Angular  coarse  to  fine  sand, 
little  medium  to  fine  gravel, 
trace(+)  silt  and  clay. 

94.2 

Cabrini  Gulch 

Green  Mountain 

2.23 

12,667 

80,400 

Coarse  to  fine  sand  and  silt, 
little  medium  to  fine  gravel. 

16.9 

Basin 

1.35 

11,800 

66,300 

Coarse  to  fine  sand,  little 
silt,  trace(+)  medium  to  fine 
gravel. 

19.1 

a 

Average  soil  descriptions  for  all  infiltrometer  test  plats 
within  each  basin. 

bMoist  soil. 
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State,  Quadrangle, 
Basin 


Infiltration 
Capacity, 
Inches /hr. 


TABLE  4  (2) 

Soil  Strength, 
lbs. /sq.ft. 
Wet  Dry 


Soil  Description 


COLORADO 


Nederland  Quadrangle 


Hicks  Gulch 

3.7 

8,700 

15,900 

Coarse  to  fine  sand,  some 
silt,  little  medium  to  fine 
gravel.  (Residual  soil.) 

Moon  Gulch 

NEW  MEXICO 

Bland  Quadrangle 

2.75 

11,000 

27,450 

Coarse  to  fine  sand,  little 
medium  to  fine  gravel, 
little  silt.  (Residual  soil.) 

Peters  Dome  Basin 

Cerro  del  Grant 

Quadrangle 

2.44 

6,330 

21,300 

Coarse  to  fine  sand  and  silt, 
little  fine  gravel. 

Cerro  Pavo  Basin 

Jarosa  Quadrangle 

2.1 

17,000 

69,800 

Sand  and  silt,  trace(+)  fine 
gravel. 

Hairpin  Basin 

2.15 

5,420 

30,400 

Coarse  to  fine(+)  sand  and 
silt. 

Poleo  Basin 

Paxton  Springs 

Quadrangle 

1.88 

5,790 

43,300 

Coarse  to  fine  sand,  and  silt, 
little  medium  to  fine  gravel. 

Quartz  Hill  Basin 

Santa  Fe  Quadrangle 

2.74 

5,330 

18,600 

Coarse  to  fine  sand,  some 
medium  to  fine  gravel,  trace 
(+)  silt  and  clay.  (Residual 
granite  soil.) 

Camino  Basin 

1.9 

5,112 

32,500 

Coarse  to  fine  sand,  little 
silt,  little  medium  to  fine 
rounded  gravel. 

Tano  Basin 

Sedillo  Quadrangle 

1.78 

6,300 

17,000 

Coarse  to  fine  sand,  some 
(-)  silt,  trace(+)  medium  to 
fine  gravel. 

Cedro  Peak  Basin 

1.99 

4,850 

17,500 

Coarse  to  fine  sand,  little 
silt  and  clay,  trace(+)  med¬ 
ium  to  fine  gravel. 

Roughness 

Number, 

Inches 


18.0 

31.4 

16.46 

1.29 

.44 

2.90 

23.0 

83.7 

28.1 

28.4 
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State,  Quadrangle, 
Basin 

UTAH 

Lehi  Quadrangle 

Dry  Hollow 

Springville  Quadrangle 

Spring  Creek 

Tributary  to 
Hobble  Creek 


Infiltration 

Capacity, 

Inches/hr. 


TABLE  4  (3) 

Soil  Strength, 
lbs  ./sq.ft. 
Wet  Dry 


Soil  Description 


Roughness 

Number, 

Inches 


1.70 

9,400 

30,200 

Fine  sand,  some  silt  and 
clay,  little  fine  gravel. 

66.3 

3.05 

2,650 

26,200 

Coarse  to  fine  sand,  and  silt 
and  clay,  trace(+)  gravel. 

14.21 

1.95 

1,665 

32,100 

Fine  sand,  and  silt  and  clay, 
little(-)  medium  to  fine 

52.48 

gravel . 
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TABLE  5 

SUMMARY  OF  CLIMATIC  DATA 


A:  Average  number  of  days  per  year  with  rain. 

B:  P-E  index. 

C:  5-year  1-hour  rain,  in  inches. 

D:  January  precipitation  in  terms  of  1/12  mean  annual  precipitation. 
E:  Average  length  of  growing  season. 

F:  Elevation  of  station. 

*  Indicates  estimate. 


State,  Region, 
Station 


Mean  Annual  Mean  Annual 

Precipitation  Temperature  A  B 


C  D  E 


F 


ARIZONA 


North  Central 


Prescott 

18.61 

52.9 

54 

39.3 

1.57 

1.14 

143 

5700 

Northeast 

Tuba  City 

Wupatki 

6.88 

55.1 

37 

12.8 

.77 

.95 

179 

4936 

National  Monument 

7.78 

57.5 

50* 

11.0 

.77 

.71 

180 

4908 

Southeast 

Benson 

10.63 

63.8 

35 

- 

1.40 

.64 

225 

3500 

Canelo  RS 

18.96 

58.5 

71 

30.4 

1.54 

.78 

180 

5000 

Crown  C  Ranch 

16.08 

- 

55* 

23.6* 

1.82 

.84 

210 

4775 

Ft.  Huachuca 

16.96 

61.4 

49 

25.6 

1.54 

.94 

232 

5100 

Pantano-Vail 

11.65 

66.6 

26 

15.4 

1.68 

.89 

286 

3300 

Rosemont 

17.73 

- 

45 

- 

1.82 

.91 

227 

4800 

Santa  Rita 

Exp.  Range 

19.73 

63.8 

58 

“ 

1.96 

.99 

210 

4300 

Tucson  Magnetic 
Observatory 

10.98 

66.8 

50* 

14.6 

1.75 

.92 

240 

2526 

COLORADO 


North  Central 


Allenspark 

22.31 

40.3 

80* 

71.5 

.84 

.67 

90 

8500 

Boulder 

18.88 

50.9 

74 

- 

1.54 

.48 

161 

5404 

Denver 

13.43 

50.1 

85 

30.0 

1.19 

.45 

171 

5221 

Edgewater 

16.17 

50.0 

65 

- 

1.26 

.27 

146 

5450 

Frances 

25.02 

39.9 

132 

82.2 

1.05 

.53 

108 

9300 

Golden 

18.61 

- 

80* 

- 

1.26 

.34 

- 

5799 

Gold  Hill 

25.56 

- 

90* 

- 

1.26 

.51 

139 

8630 

Hawthorne 

21.21 

- 

78 

- 

1.54 

.49 

120 

5923 

Idaho  Springs 

15.27 

43 .2 

78 

- 

.98 

.29 

125 

7543 

Longs  Peak 

21.51 

37.6 

89 

78.3 

.81 

.48 

59 

9000 

Morrison 

18.14 

47.8* 

52 

40.3* 

1.33 

.30 

- 

6120 

Mt.  Morrison 

15.69 

- 

- 

- 

1.09 

.99 

- 

6000 

Silver  Lake 

28.17 

29.3* 

94 

128.0* 

.84 

.80 

30 

10200 

South  Central 

Madrid 

14.75 

49.5* 

81 

27.6* 

1.54 

.17 

130 

6364 

Mule  Shoe  Lodge 

16.94 

- 

- 

- 

1.33 

1.24 

130 

8800 

North  Lake 

21.34 

- 

88 

- 

- 

.56 

- 

8800 

Stonewall 

17.97 

- 

81 

- 

1.40 

.13 

110 

8000 

Tercio 

17.21 

- 

- 

- 

1.26 

.66 

100 

8040 

Trinidad  AP 

14.02 

50.8 

55 

- 

1.65 

.42 

142 

5746 

Trinidad  CY 

16.20 

51.4 

69 

29.8 

1.65 

.47 

167 

6030 
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State,  Region, 
Station 


TABLE  5  (2) 

Mean  Annual  Mean  Annual 

Precipitation  Temperature  A  B  C  D  E  F 


COLORADO 


Lake  Maloya,a 

21.77 

49.16 

- 

52.4 

- 

.51 

140 

7400 

New  Mexico 

NEW  MEXICO 

Northern  Rio  Grande 

Aspen  Grove  Ranch 

23.61 

38.17* 

92 

97.2* 

- 

.85 

104 

9500 

Bandelier 

15.60 

50.22 

45 

31.2 

1.40 

.54 

- 

6061 

Bland 

19.08 

- 

79 

-  ; 

1.40 

.75 

- 

7050 

Capulin  RS 

15.94 

45.96 

- 

31.05 

1.05 

.31 

143 

7200 

Edgewood 

11.08 

- 

- 

- 

- 

.55 

- 

6800 

Lazy  Ray  Ranch 

18.70 

- 

- 

- 

- 

1.12 

- 

8172 

Lee  Ranch 

22.45 

38.9 

50 

69.9 

- 

.36 

98 

9000 

Los  Alamos 

18.32 

47.7 

73 

78.9 

1.40 

.59 

168 

7324 

Otto  AP 

11.58 

49.2 

- 

21.5 

1.36 

.47 

120 

6226 

Red  River 

22.04 

39.9 

94 

77.5 

- 

.71 

88 

8956 

Sandia  Park 

19.20 

48.9 

- 

43.5 

- 

.89 

- 

7011 

Santa  Fe  AP 

10.80 

50.4 

- 

19.1 

1.13 

.36 

- 

6312 

Santa  Fe  CY 

14.29 

49.1 

86 

26.3 

1.12 

.53 

178 

7045 

Tijeras  RS 

14.57 

49.1 

63 

32.4 

1.47 

.83 

136 

6300 

Wolf  Canyon 

22.35 

39.79 

- 

80.5 

- 

.83 

- 

8150 

Albuquerque  CY 

9.05 

55.8 

47 

11.3 

.95 

.51 

198 

5130 

Albuquerque  AP 

8.48 

56 .6 

- 

- 

.98 

.39 

- 

5310 

Zuni  Mountains 

Bluewater 

10.35 

47.6 

44 

19.3 

.91 

.51 

122 

6650 

El  Morro  National 

13.07 

47.6 

- 

30.1 

1.26 

.85 

- 

7218 

Monument 

Grants 

8.46 

51.5 

- 

13.8 

.98 

.60 

- 

6500 

San  Rafael 

14.31 

57.2 

56 

26.6 

.98 

.30 

- 

6509 

Gila  and  Southwest 

Bear  Creek  Ranch 

13.08 

- 

- 

- 

- 

1.14 

- 

5300 

Buckhorn 

11.02 

- 

- 

- 

- 

1.14 

161 

4900 

Cliff 

14.35 

57.9 

52 

- 

- 

.78 

171 

4800 

Cureton  Ranch 

13.11 

- 

- 

- 

- 

1.01 

- 

5200 

Ft.  Bayard 

16.38 

55.0 

68 

27.1 

1.4 

.63 

194 

6152 

Mimbres  RS 

17.36 

- 

50 

- 

1.33 

.82 

162 

5754 

Pinos  Altos 

23.08 

51.42* 

95 

45.2* 

1.54 

.70 

127 

7000 

Silver  City 

16.83 

53.4 

67 

30.7 

1.54 

.74 

180 

6051 

Tyrone 

14.98 

- 

49 

- 

- 

.70 

190 

5923 

Tyrone  Station 

24.28 

49.99 

- 

54.0 

1.61 

.86 

180 

7986 

White  Signal 

13.59 

- 

- 

- 

1.58 

1.13 

210 

6070 

Whitewater 

7.32 

- 

- 

- 

- 

.33 

- 

5150 

UTAH 

North  Central 

Alpine 

15.94 

46.9* 

57 

48.0* 

.77 

1.14 

100 

5000 

Alta 

57.33 

36.7° 

- 

128.0 

.77 

1.79 

90 

8760 

Birdseye 

12.37 

- 

- 

- 

- 

1.16 

- 

5740 

Deer  Creek  Dam 

22.10 

43.2 

- 

101.9 

- 

1.42 

- 

5285 

East  Portal 

20.82 

46.7 

- 

- 

- 

1.45* 

- 

7606 

^Borders  on  south  central  Colorado  area. 
No  May -June  records. 

No  December-April  records. 
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State,  Region, 
Station 


TABLE  5  (3) 

Mean  Annual  Mean  Annual 

Precipitation  Temperature  A  B  C  D  E  F 


UTAH 


North  Central 


Lower  American  Fork 

16.73 

51.9 

Maplewood 

20.18 

- 

Provo  Bench 

17.33 

48.7 

Provo  Radio 

14.76 

48.3 

Soldier  Summit 

14.04 

39.03 

Spanish  Fork 

17.55 

52.0 

Thistle 

17.16 

46.3 

Timpanogos  Cave 

24.92 

49.2 

Timpanogos  Summit 

32.92 

- 

Utah  Lake  -  Lehi 

12 .25 

48.5 

West  Portal 

28.79 

- 

- 

42.1 

- 

1.22 

168 

5063 

- 

- 

- 

1.06 

- 

4890 

- 

49.6 

.70 

1.06 

140 

5000 

66 

42.1 

.70 

1.25 

122 

4650 

- 

65.3 

.84 

1.22 

60 

7490 

83 

- 

.70 

.98 

162 

4711 

69 

57.0 

.70 

1.50 

90 

5050 

- 

71.8 

- 

1.10 

- 

5523 

- 

- 

- 

2.03 

- 

8300 

- 

33.6 

.77 

1.10 

135 

4497 

- 

- 

- 

1.10 

- 

7550 

TABLE  6 


LIST  OF  WITHIN-GROUPS  VARIANCES 

1.  Maximum  Valley -Side  Slope 

A.  Field  measurements* 

1)  Within  basins 

N  =  880,  k  =  22,  s2  =  19.57,  s  =  4.4° 

2)  Within  lithologies 

N  =  880,  k  =  6,  s2  =  126.19,  s  =  11.2° 

B.  Map  measurements 

Within  basins 

N  -  413,  k  -  24,  s2  =  31.44,  s  =  5.6° 

2.  Percent  Bare  Area 

Within  basins  (individual  field  measurements) 

N  =  222,  k  =  21,  s2  =  206.65,  s  =  14.4% 

3.  Infiltration  Capacity 

A.  Within  basins  (individual  field  measurements) 

N  =  122,  k  =  22,  s2  =  1.16,  s  =  1.08  in/hr. 

B.  Within  lithologies  (individual  field  measurements) 

N  =  122,  k  =  6,  s2  =  1.31,  s  =  1.14  in/hr. 

4.  Roughness  Number 

Within  lithologies 

N  =  22,  k  =  6,  s2  =  628.8,  s  =  25.08  inches 

5.  Drainage  Density 

A.  Within  homogeneous,  contiguous  areas 

N  =  100,  k  =  24,  s2  =  68.09,  s  =  8.25  per  mile 

B.  Within  lithologies 

N  =  106,  k  =  6,  s2  =  407.33,  s  =  20.18  per  mile 

6.  Relative  Relief 

Similar  basins  within  homogeneous,  contiguous  areas 
N  =  106,  k  =  25,  s2  =  3.24,  s  =  1.80 

7.  Ruggedness  Number 

Similar  basins  within  homogeneous,  contiguous  areas 
N  =  108,  k  =  26,  s2  -  0.596,  s  =  0.772 

*N  =  Total  number  of  observations 
k  =  Number  of  basins,  lithologic  classes,  etc. 
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TABLE  7 

INDIVIDUAL  INFILTRATION-CAPACITY  MEASUREMENTS 


BASIN 

INFILT. 

CAP. 

SLOPE  OF 
PLAT  (°) 

TYPE  OF  SURFACE 

Arizona 

CHINLE  BADLANDS  SPECIAL  MAP 

N.  Basin  (shale) 

1.1 

13 

Bare,  U* 

S.  Basin 

1.2 

26 

Bare,  U 

1.1 

27 

Bare,  U 

1.0 

7 

Bare,  U 

PRESCOTT  QUADRANGLE 

Whipple  Basin  (Granite,  alluvium) 

1.6 

4 

Bare,  U 

1.5 

2 

Sparse  grass,  U 

1.7 

6 

Bare,  U 

1.8 

9 

Bare,  U 

SAGUARO  MON.  SPECIAL  MAP  (schist) 

N.  Basin 

0.9 

10 

Bare,  U 

1.1 

12 

Bare,  U 

S.  Basin 

0.9 

7 

Bare,  U 

1.1 

2 

Bare,  U 

SONOITA  QUADRANGLE 

Hog  Canyon  Tributary  (alluvium) 

1.6 

1 

Loose  leaves,  R| 

2.5 

1 

V  »  » 

2.3 

1 

n  »  » 

1.5 

1 

Bare,  U 

1.9 

1 

Loose  leaves,  R 

6.1 

1 

A0  removed 

1.5 

13 

Loose  leaves,  R 

5.2 

1 

Sparse  grass,  U 

3.1 

1 

Bare,  U 

1.0 

0 

Bare,  U 

2.0 

0 

Loose  leaves,  R? 

Colorado 

MORRISON  QUADRANGLE 

Cabrini  Gulch  (schist) 

2.1 

8 

Bare,  U 

2.8 

12 

Grass  litter,  U 

1.8 

14 

Grass,  U 

Green  Mountain  (alluvium) 

1.3 

8 

Grass, U 

1.7 

18 

Grass, U 

1.1 

4 

1/2  grass,  1/2  bare, 

1.3 

16 

Grass,  U 

NEDERLAND  QUADRANGLE 

Hicks  Gulch  (granitic  gneiss) 

4.4 

35 

Grass, U 

3.3 

4 

Pine  litter,  U 

3.0 

12 

”  ”  ,  R 

4.5 

22 

Bare,  U 

4.7 

32 

Grass,  pine  litter,  U 

5.3 

18 

Pine  duff  (A  ),  R 

4.2 

22 

»  »  r>  » 

2.7 

17 

»  »  »  v 

2.8 

19 

»  »  r?  » 

*U  =  undisturbed  surface 

1.8 

4 

Y)  r>  r>  » 

fR  =  removed 
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TABLE  7  (2) 


BASIN 

INFILT. 

CAP. 

SLOPE  OF 
PLAT  (°) 

TYPE  OF  SURFACE 

Moon  Gulch  (gneiss,  granite) 

2.2 

8 

Grass, U 

1.8 

12 

Pine  duff  (Aco  ),  R 

3.1 

9 

n  v  99  » 

4.1 

14 

v  79  y>  99 

2.4 

1 

Bare,  U 

1.2 

7 

Pine  duff  (AQO),  R 

3.4 

3 

79  99  99  99 

3.8 

5 

y>  r>  n  r> 

New  Mexico 

BLAND  QUADRANGLE 

Peters  Dome  Basin  (pumice,  rhyolite) 

2.9 

9 

Pine  duff  (AQO),  R 

3.0 

13 

Fir  duff,  R 

2.9 

17 

Pine  duff  (Aq,-,  ),  R 

2.3 

14 

99  79  99  99 

1.9 

12 

99  99  99  99 

1.5 

9 

79  79  99  99 

2.6 

24 

99  99  99  99 

CERRO  DEL  GRANT  QUADRANGLE 

Cerro  Pavo  Basin  (scoria) 

5.3 

7 

Fir  needles,  R 

2.3 

9 

Grass,  U 

1.9 

3 

Grass,  U 

0.5? 

3 

Fir  needles,  R 

1.9 

8 

Loose  leaves,  R 

1.9 

2 

Grass,  U 

1.5 

2  1/2 

Pine  duff,  R 

2.0 

9 

Grass,  U 

1.7 

10 

Leaves, R 

1.7 

8 

Grass, U 

JAROSA  QUADRANGLE 

Hairpin  Basin  (sandstone) 

2.0 

17 

Pine  duff,  R 

2.9 

10 

Grass;  leaves  R 

1.9 

6 

Grass;  U 

1.8 

6 

Pine  duff,  R 

Poleo  Basin  (sandstone,  limestone) 

1.7 

3 

Pine  duff,  R 

2.0 

10 

79  99  99 

3.6 

12 

79  79  99 

1.7 

11 

79  99  79 

1.5 

9 

79  79  99 

1.2 

18 

Grass, U 

1.4 

9 

79  79 

1.7 

8 

97  79 

2.1 

11 

Pine  duff,  R 

1.9 

9 

79  79 

PAXTON  SPRINGS  QUADRANGLE 

Quartz  Hill  Basin  (granite) 

2.0 

11 

Bare,  U 

1.8 

6 

79  99 

1.8 

5 

79  79 

3.3 

8 

99  77 

3.5 

12 

Pine  duff,  R 

1.9 

7 

Bare,  U 

6.5 

13 

79  79 
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TABLE  7  (3) 


BASIN 

INFILT. 

CAP. 

SLOPE  OF 
PLAT  (°) 

TYPE  OF  SURI 

SANTA  FE  QUADRANGLE 

Camino  Basin  (alluvium) 

2.0 

6 

Bare,  U 

1.4 

10 

Grass,  U 

2.3 

3 

Bare,  U 

Tano  Basin  (alluvium) 

1.6 

8 

Bare,  U 

2.8 

2 

99  99 

1.4 

5 

99  99 

1.6 

4 

Sparse  grass,  U 

1.5 

1 

t>  99  99 

SEDILLO  QUADRANGLE 

Cedro  Peak  Basin  (sandstone,  shale) 

2.1 

10 

Bare,  U 

1.2 

11 

Bare,  U 

2.0 

15 

Pine  duff,  R 

1.4 

6 

79  99  99 

1.5 

4 

79  99  99 

2.9 

14 

79  79  99 

2.5 

6 

99  99  99 

2.3 

28 

99  99  99 

Utah 

LEHI QUADRANGLE 

Dry  Hollow  (quartzite,  volcanic) 


SPRINGVILLE  QUADRANGLE 
Hobble  Cr.Trib.  (sandstone) 

Spring  Creek  (limestone) 


0.9 

26 

Grass, U 

2.2 

20 

Bare,  U 

2.5 

5 

99  99 

0.9 

7 

99  99 

1.6 

12 

99  99 

2.1 

21 

Grass,  U;  oak  lit.  R 

2.3 

22 

Grass,  U 

1.6 

41 

Bare,  U 

2.9 

14 

Grass, U 

2.6 

30? 

99  99 

4.3 

28 

99  99 

6.8 

20 

Leaves,  grass,  U 

1.8 

26 

Grass, U 

2.7 

24 

Leaves,  grass,  U 

1.9 

19 

Leaves,  weeds,  U 

1.4 

? 

Grass; leaves,  R 

Additional  measurements  on  residual  granite  soils: 
Arizona 

HARSH  AW  (PATAGONIA  MTNS.)  QUAD,  (granite) 


IRON  SPGS.  (SIERRA  PRIETA)  QUAD,  (granite) 


2.3 

1.9 
5.7 

5.3 

2.5 

2.9 
1.8? 
2.0 

2.6 
6.2 
4.5 

1.9 


? 

Bare 

99 

Grass,  R 
Pine  duff,  R 

7)  yy  99 

Bare,  U 
Pine  duff,  R 
? 

Pine  duff,  R 
Bare,  U 
Pine  duff,  R 
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TABLE  7  (4) 


BASIN 

INFILT. 

SLOPE  OF 

TYPE  OF  SURFACE 

CAP. 

PLAT  (°) 

Arizona 

KIRKLAND  (SIERRA  PRIETA)  QUAD,  (granite) 

2.3 

Pine  duff,  R 

2.1 

99  99  99 

New  Mexico 

BIG  BURRO  MTNS.  QUADRANGLE  (granite) 

4.5 

Leaves, R 

7.8 

Bare,  U 

4.0 

99  99 

3.7 

r>  v 

8.3+ 

99  r) 

8.3+ 

Pine  duff,  R 

FORT  BAYARD  (PINOS  ALTOS)  QUAD,  (granite) 

3.5 

Bare,  U 

1.8 

99  99 

1.5 

99  99 

TWIN  SISTERS  (PINOS  ALTOS  MTNS.)  QUAD. 

2.1 

Bare,  U 

(granite) 

3.5 

99  99 

